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ABSTRACT 

The  object  of  this  thesis  is  to  study  several  topics 
in  the  general  theory  of  relativity  which  involve  space- 
time  regions  of  intense  gravitational  field. 

Salient  features  of  the  Schwarzschild-Kruskal  mani¬ 
fold  are  reviewed.  A  system  of  null  coordinates  of 
Kruskal’s  type  for  a  collapsing  spherical  cloud  of  dust  is 
presented  which  gives  a  singularity-free  description  of 
the  complete  space-time  right  down  to  the  singular  event 
of  maximum  implosion.  Questions  regarding  the  qualitative 
effects  of  departures  from  sphericity  on  gravitational 
collapse  are  discussed  in  the  light  of  some  recent  advances. 
A  solution  of  Einstein’s  field  equations  is  derived  which 
represents  a  class  of  thin  spherical  shells  of  charged 
dust  collapsing  (and,  in  general,  bouncing)  in  the  field  of 
a  charged  spherical  body  placed  at  the  center.  A  method 
is  described  which  develops  the  interior  field  and  physical 
properties  of  a  slowly  spinning,  nearly  spherical  mass  shell 
as  a  power  series  in  the  angular  velocity  go  ,  on  the  pre¬ 
sumption  that  the  exterior  field  is  the  Kerr  metric;  results 
are  worked  out  explicitly,  correct  to  the  third  order  in 
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CHAPTER  I 

INTRODUCTION 

In  scope,  depth,  and  logical  simplicity,  the 
general  theory  of  relativity  ^^(GR)  today  still  stands  as 
the  best  (classical)  theory  of  space,  time,  and  gravitation. 
Created  in  1915  by  Einstein,  it  has  since  stood  its  ground 
against  the  tests  of  empirical  observations  (few  though 
there  are)  and  arguments  of  principle.  And  although  on 
empirical  grounds  there  might  appear  little  to  distinguish 
GR  from  some  other  theories  ’  (of  gravitation),  in  assess¬ 
ing  the  inherent  plausibility  of  these  later  theories  one 
must  bear  in  mind  that  they  were  formulated  with  full  know¬ 
ledge  of  the  important  effects  (such  as  the  deflection  of 
light  near  a  mass  center)  already  known  to  exist,  after 
having  been  first  predicted  by  GR  and  subsequently  con¬ 
firmed  . 


By  showing  how  gravitation  may  be  reduced  to 
being  an  aspect  of  space-time  geometry,  and  by  bringing 
down,  in  turn,  this  geometry  from  the  lofty  height  of  the 
a  priori,  Einstein  has  effected  a  profound  modification  on 
our  concept  of  the  external  world.  But  owing  partly  to  the 
intrinsic  mathematical  difficulties  of  the  theory,  as  mani¬ 
fested  most  characteristically  in  the  nonlinearity  of  its 
field  equations,  and  partly  to  the  extreme  sophistication 
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needed  to  experimentally  check  those  of  its  predictions  that 
can  he  worked  out,  subsequent  progress  of  the  theory  has  not 
been  as  great  as  the  initial  enthusiasm  which  it  had  genera¬ 
ted  might  have  led  one  to  expect.  Moreover,  since  under 
normal  situations  the  gravitational  field  is  incomparably 
weaker  than  the  other  so-called  fundamental  fields,  many 
physicists  are  inclined  to  believe  that  GR  would  never  be¬ 
come  important  on  a  microscopic  level.  Such  a  consideration, 
however,  in  view  of  our  present  state  of  understanding, 
would  appear  somewhat  premature;  in  particular,  the  essen¬ 
tial  non-linear  character  of  GR  may  yet  lead  to  unsuspected 
( 4 ) 

subtleties.  It  may  be  mentioned  that  one  of  the  reasons 

responsible  for  the  considerable  rise  in  renewed  interest 

in  GR  during  the  last  fifteen  years  or  so  is  precisely  the 

lack  of  essential  progress  in  other  branches  of  theoretical 

physics:  some  hope  that,  perhaps  by  suitably  amalgamating 

GR  and  quantum  theory,  a  satisfactory  theory  of  elementary 

(4) 

particles  will  emerge. 

Another  reason  has  to  do  with  recent  or  imminent 

technological  advances  which  will  permit  one  to  increase 

the  accuracy  of  earlier  experiments  and  also  to  test  other 

( 5 ) 

predictions  of  GR .  The  interaction  between  theories 

and  experiments  is  certainly  vital  to  the  growth  of  any 
empirical  science.  But  it  is  well  to  note  again  that, 
except  in  cosmology,  almost  all  of  the  physical  predictions 
of  GR,  as  normally  envisaged,  involve  geometries  that 
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deviate  extremely  little  from  flatness  -  a  condition 

■t 

sufficient  for  the  Newtonian  gravitational  theory  to  hold 
with  tremendous  accuracies  (provided  the  velocities  of 
material  particles  are  small  in  comparison  with  the  speed 
of  light  c).  Only  when  the  gravitational  field  is  intense, 
i.e.,  when  the  Newtonian  potential  <J)  becomes  comparable 
with  c2,  do  we  expect  significant  differences  between  the 
predictions  of  GR  on  the  one  hand,  and,  on  the  other  hand, 
those  of  Newtonian  and  indeed  other  gravitational  theories. 

Nor  a  collection  of  matter  of  mass  m  and  characteristic 
dimension  R, 

|  tp  |  / c  2  1  Gm/c  2R 

at  the  surface,  where  G  is  the  Newtonian  gravitational 
constant.  At  the  surface  of  the  sun  for  example,  we  have 
Gm/c2R  ~  10“6.  One  would  therefore  look  to  astrophysics 
for  possible  processes  where  general  relativistic  effects 
might  play  an  essential  role. 

Such  a  process  could  conceivably  be  taking  place 

* 

in  "Quasi-st ellar  objects"  -  astrophysical  objects 

discovered  since  I960  and  probably  belonging  to  a  class  of 

their  own  with  no  parallel  in  known  observations.  Though 

these  objects  are  now  believed  to  be  much  smaller  in  size 

4  4 

than  a  typical  galaxy,  their  radio  power  is  about  10  ergs/sec 


For  a  recent  monograph  on  Quasi-stellar  objects,  see 
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and  their  optical  power  about  10  ergs/sec,  which  is  100 

times  the  total  energy  output  of  a  giant  galaxy.  Other 

peculiarities  such  as  that  connected  with  their  light 

variations  have  also  been  observed,  but  the  central  problem 

of  the  moment  is  to  find  a  mechanism  to  account  for  the 

energy  output.  Reasons  can  be  adduced  to  practically  rule 

out  thermonuclear  reactions  as  being  responsible  for  the 

energy  output,  and  so  attention  has  turned  to  gravitation. 

The  reason  is  that  gravitational  energy  depends  on  Gm2/’R, 

i.e.,  on  the  square  of  the  mass,  so  that  it  can  play  a 

decisive  role  in  the  evolution  of  a  very  massive  body 

5  8 

(m  ~  10  -  10  solar  masses).  Thus  a  model  for  quasi- 

( 7  ) 

stellar  objects  has  been  proposed  which  consists  of  a 
very  massive  body  undergoing  violent  gravitational  contrac¬ 
tion  (collapse).  The  gravitational  collapse  of  this  super- 
star  could  perhaps  supply  the  necessary  energy  if  it  were 
to  shrink  to  dimensions  such  that  the  Newtonian  potential 

became  comparable  with  c2  -  thus  necessitating  the 

employment  of  the  full  machinery  of  GR . 


( 8 ) 

Relativistic  gravitational  collapse1  ;  by  itself 
forms  a  fascinating  study,  for  it  has  brought  to  light 
peculiar  features  wholly  unexpected  on  the  basis  of 
Newtonian  theory.  Chapter  III  gives  a  quantitative  summary 
of  the  collapse  of  a  model  which  is  particularly  simple 
but  generally  thought  to  be  representative:  the  collapse 
of  a  spherical  ball  of  dust .  In  that  connection  we  shall 
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introduce  a  coordinate  system  of  the  Kruskal^^  type,  so 
as  to  present  the  complete  history  of  the  collapse  in  a 
unified  picture.  The  chapter  ends  with  a  discussion  on 
whether  or  not  our  conventional  picture  of  collapse  can 
still  be  upheld  when  the  assumption  of  strict  spherical 

symmetry  -  which  is  of  course  never  realized  in  actual 

bodies  -  is  dropped. 

A  collapsing  ball  of  dust  eventually  ends  in  a 
state  of  singularity  characterized  by  infinite  curvature 
and  energy  density,  at  which  the  analysis  must  come  to  halt. 
This  state  of  affairs  can  be  avoided  in  the  case  of  charged 
bodies,  though  apparent  paradoxes  may  arise.  The  collapse 
of  charged  spherical  shells  is  dealt  with  in  Chapter  IV. 

Chapter  II  contains  a  summary  of  some  of  the  most 

important  aspects  of  the  Schwarzschild  (exterior)  solution  ^  ^ 

including  its  analytic  completion^'*' ^  This  remarkable 

solution  represents  the  vacuum  field  outside  any  spheri- 

(12) 

cally  symmetric  (uncharged)  distribution  of  matter  ,  and 
is  one  of  the  only  two  solutions  of  Einstein’s  gravita¬ 
tional  equations  which  have  any  experimental  basis  (the 

(13) 

other  being  the  Friedmann  solution  in  cosmology).  The 

trajectories  of  test  particles  and  light  rays  have  already 

( l4  ) 

been  worked  out  and  may  be  compared  with  the  corres¬ 

ponding  Newtonian  situation.  In  regions  of  intense  field, 
qualitative  differences  occur.  The  Schwarzschild  solution 
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is  also  interesting  in  that  it  may  be  used  to  illustrate 
such  basic  ideas  as  completeness  of  space-time  manifold 
and  singularities,  and  the  peculiar  difficulties  associated 

with  the  problem  of  coordinates  in  GR  -  all  of  which  are 

trivial  in  the  Newtonian  picture.  Finally,  the  solution 
is  of  great  importance  for  both  the  qualitative  and  the 
quantitative  analyses  of  relativistic  gravitational  collapse. 


While  the  Schwarzschild  metric  has  been  in  exist¬ 
ence  for  more  than  fifty  years  and  may  be  complemented  by 
some  interior  solution  to  form  the  complete  geometry  of 
a  space-time  manifold  due  to,  for  example,  a  static  sphere 
of  perfect  fluid,  no  exact  solution  representing  the  com¬ 
plete  field  of  some  isolated  rotating  object  of  a  realistic 
nature  has  yet  been  exhibited,  although  at  large  distances 

from  such  an  object,  where  the  field  is  assumed  weak,  the 

(15) 


form  of  the  metric  can  be  obtained 


from  the  linearized 


Einstein  field  equations.  However,  in  1963  was  discovered 
a  remarkable  two-parameter  vacuum  field  solution  (the  Kerr 


solution 


(16) 


)  which  turned  out  to  have  features  that  might 


be  expected  of  a  space-time  exterior  to  some  object 
possessing  mass  and  angular  momentum^  ?  the  exact  nature 
being  still  unknown.  In  Chapter  V,  a  method  is  described 
which  develops  the  interior  field  and  the  physical  proper¬ 
ties  of  a  slowly  spinning,  nearly  spherical  thin  mass  shell 
as  a  power  series  in  the  angular  velocity  w,  on  the  pre¬ 
sumption  that  the  exterior  field  is  the  Kerr  metric. 
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Results  are  worked  out  explicitly  to  the  third  order  in  w, 

and  we  extend  a  recent  strong-field  study  by  Brill  and 
( 19 ) 

Cohen  on  the  classical  problem  of  the  "relativity  of 

inertia"  in  the  sense  of  Mach^^  (see  Part  A,  Chapter  V), 

(  21  ) 

first  formulated  and  worked  out  by  H.  Thirring  in  1918 

within  the  framework  of  GR,  but  using  the  weak  field  approx¬ 
imations  . 
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CHAPTER  II 

SPHERICALLY  SYMMETRIC  VACUUM  GRAVITATIONAL  FIELD 

1 .  General  Relativity:  A  Summary. 

The  two  most  important  postulates  underlying 
Einstein’s  theory  of  gravitation  are,  firstly,  that  the 
geometry  of  space-time  is  Riemannian,  and,  secondly,  that 
of  all  mathematically  possible  Riemannian  spaces,  only 
those  that  are  related  to  matter  distributions  in  a  specific 
way  can  be  physically  realized.  A  few  standard  geometrical 
concepts  and  relations  will  be  needed  to  make  the  state¬ 
ment  more  precise. 

A  Riemannian  space  may  be  defined  as  a  geometri¬ 
cal  space  in  which  the  "distance”  or  "interval"  ds  between 
two  neighbouring  points  xa  and  xa+dxa  (referred  to  a 
suitable  coordinate  system)  is  given  by  a  (not  necessarily 
positive-definite)  quadratic  differential  form 

ds2  =  g  Q(x)  dxadx^  (1) 

a  B 

where  g  „  ,  the  metric  tensor,  behaves  like  a  covariant 
a  B 

tensor  of  the  second  rank  under  coordinate  transformations, 
and  is  assumed  non-singular  in  the  sense  of  matrix  theory: 

g  =  det(gaf3)  1  0  . 


(2) 
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Condition  (2)  enables  us  to  define  the  contra- 
variant  metric  tensor  ga^  by 


a3  xa 

s  sBy  ■  «Y  • 


(3) 


a  S 

By  convention,  g  0  and  g  are  used  to  lower 

j  °a3 

or  raise  tensorial  indices;  for  example,  with  a  covariant 

vector  A  one  can  associate  a  contravariant  vector 
a 

Aa  =  ga®Ag  . 


5a3 


The  following  entities  are  constructed  out  of 
a  6 

and  g  and  their  derivatives: 


Christoffel  symbols  of  the  first  and  second  kind  (res¬ 
pectively  )  : 


^a,3y  ^^y^a3  +  ^3^ay  ^a®3y^  * 


ra  _  aSr 

1 3y  s  o , 3y  5 


(4) 

(5) 


where  3  =  3/3xa  . 

a 


Riemann  tensor: 


=  3  (r“ 

3y6  y  3o 


3^(r“  )  +  r*  r?1  -  r“ 

6  3y  36  Ay  3y  A  6 


(6) 


Ricci  tensor 


R  „  =  Ry 


a3 


a3y 


=  3g3  3Q[ln(-g)]  -  ^7 Y  3  [ In ( -g)  ]  - 


a  3 


a3  y 


-  3  rYQ  +  rY  r5R  . 

y  a3  oa  y3 


(7) 
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Curvature  invariant: 


R  =  ga6R 


a8 


(8) 


Einstein’s  tensor: 


Go=R0-^g0R 

aB  aB  aB 


(9) 


For  vector  and  tensor  fields  defined  throughout 
a  domain  of  space-time,  the  covariant  derivatives  (indica¬ 
ted  by  a  vertical  stroke)  are  as  follows: 


Aa  =  3  Aa  +  ra  Ay 
|  B  B  By  ’ 


A 


=  -  rIoA„  , 


a  B  B  a 


aB  y 


(10) 


(11) 


=  3YAa3  +  r“6A3<S  +  f^6  Aa<s  ,  (etc.).  (12) 


We  now  give  a  brief  summary  of  general  rela¬ 
tivity  . 

Space-time  is  represented  by  a  four-dimensional 
normal  hyperbolic  Riemannian  space  of  signature  (+++-) . 
This  implies  that  at  any  given  event  E,  the  metric  tensor 
can  be  reduced  to  the  Minkowskian  form 


>aB 


Cgag)  =  diag  (1,1,1,-c2) 
E 


Here  and  in  the  following,  greek  indices  run  from  1  to  4. 
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by  a  suitable  coordinate  transformation  (c  is  the  funda¬ 
mental  velocity). 


The  physical  interpretation  of  the  metric  is 
given  as  follows.  Let  E  and  E'  be  neighbouring  events 
x ^  and  +  dx^1 .  An  observer  (with  four-velocity  ) 

momentarily  at  E  measures  the  spatial  separation  of  the 
two  events  as 

[(g  +  V  V  /c2)dx^dxv] 

y  v  y  v 

and  their  time-difference  as 


V  dxy  . 
y 


The  geometry  of  space-time  is  related  locally  to 
the  distribution  of  matter  (characterized  phenomenologically 
by  an  energy-momentum  tensor  T„  )  through  Einstein's 

gravitational  field  equations 


Gae  =  -K  Tag  ’  K  E  8*G/0*  >  (13) 

in  which  G  is  the  Newtonian  gravitational  constant.  The 
existence  of  four  differential  identities 


Ga^ | D  =  0  (Contracted  Bianchi  Identities) 

I  p 

means  that  a  consistent  choice  of  T^  must  satisfy  the 
"conservation  equations" 


0  . 


3 


(14) 
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Thus,  given  ,  eqs .  (13)  are  a  system  of  six  (non¬ 

linear)  partial  differential  equations  for  the  ten  independ¬ 
ent  components  of  g^  .  This  amount  of  indeterminacy  is, 
in  fact,  just  right,  on  account  of  the  arbitrariness  of  the 
coordinate  system. 

The  field  equations  (13)  are  the  analogue  of  the 
Poisson  equation  of  Newtonian  gravitational  theory, 

V2(f>  =  4t rGy 

(where  V2  is  the  3-dimensional  Laplacian  operator,  <J> 
the  gravitational  potential,  and  y  the  mass  density), 
which  satisfyingly  represents  a  limiting  case  of  (13)  for 
weak,  quasi-static  fields,  for  which  the  general  relativis¬ 
tic  line-element  reduces  to 

ds  2  =  ( l-2cj)/c  2  )  (dx2+dy  2+dz 2b ( l+2(j)/c  2  ) c  2dt 2  (15) 


in  appropriate  coordinates. 


equations 


For  vacuum  space-times  (T  =  0),  the 

ot  p 

reduce  to  Ga^=0  or,  equivalently. 


field 


(16) 


An  important  material  system  figuring  in  the 
study  of  general  relativity  is  the  perfect  fluid,  formally 
defined  by  Ta®  having  the  form 
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TaB  =  ( y+P/c 2 )uau^  +  P  ga3  ,  (17) 

ot  a 

in  which  u  =  dx  (x)/dx  (where  x  is  the  proper  time) 
is  the  four-velocity  of  macroscopic  fluid  particles,  P 
the  isotropic  pressure,  and  y  the  mass-energy  density 
as  measured  by  an  observer  moving  with  the  fluid  element. 
When  P  =  0  ,  we  speak  of  "dusts" . 

From  (17)  we  observe  that 

T  uD  =  -yc  u  (u  u  =  -c  )  .  (lo) 

3  a 

Thus  ua  is  the  time-like  normalized  (i.e.,  u  ua=-c2) 

*  a 

ot  6 

eigenvector  of  T  ,  with  corresponding  eigenvalue  -yc2. 

For  a  continuous  medium  having  an  arbitrary  energy-momentum 
ot  8 

tensor  T  we  shall,  following  Synge,  take  this  statement 
as  a  definition  of  its  velocity  field  ua  and  proper 
energy  density  yc2. 

Each  space-time  contains  a  class  of  geometrically 
distinguished  curves  which  form  the  closest  analogue  of 
straight  lines  in  ordinary  Euclidean  spaces.  They  are  the 
geodesics,  i.e.,  space-time  curves  whose  parametric  equa¬ 
tions  xy=xy(A)  satisfy 

d2xy/dA2  +  Ty  (dxa/dA) (dx3/dA)  =  0  .  (19) 

a3 

Here  A  is  known  as  an  affine  parameter,  and  is  deter¬ 
mined  up  to  linear  transformations.  The  geodesic 
equations  (19)  possess  the  first  integral 


' 


bnB 
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g  0 (dxa/dX ) (dx^/dA )  =  const 

ct  p 


(20) 


i .  e 


t 


a  , 

v  v  =  const  , 
a  9 

Ct  ct 

where  v  =  dx  (A)/dX  is  the  tangent  vector.  (Since  the 
metric  of  space-time  is  not  positive-definite,  the  norm 
A^A  of  a  vector  A^  can  be  positive,  zero,  or  negative; 

A^  is  then  designated  as  "space-like",  "null",  or  "time-like", 
respectively).  Thus  there  are  three  subclasses  of  geodesics: 
space-like,  null,  and  time-like  geodesics.  The  latter  two 
are  of  fundamental  importance,  as  we  shall  presently  see. 


Unlike  any  other  known  field  theory,  general 
relativity  is  one  in  which  the  equations  of  motion  for 
localized  sources  may  not  be  postulated  separately,  but 
rather  follow  from  its  field  equations.  In  particular,  we 
have  the  following  basic  result: 

In  vacuo,  the  world  line  of  a  (non-spinning  and 
uncharged)  test  particle  is  a  time-like  geodesic;  the  world 
line  of  a  "test  photon"  (or  neutrino)  is  a  null  geodesic. 

If  non-gravitational  forces  and  internal  interac¬ 
tions  can  safely  be  ignored,  a  proof  of  this  "geodesic 
hypothesis”  becomes  very  simple.  Consider  a  small  body  of 
dust  moving  in  a  vacuum  domain.  The  energy-momentum  tensor 
is  thus  Ta^  =  yuau^  ,  with  y  =  0  outside  the  world  tube 
representing  the  history  of  the  object.  Prom  the  conser- 


. 
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vation  equations  ( 1 4 )  one  immediately  obtains 


/  3  \  a  .  a  3  ~ 

(MU  ) | gu  +  u  | ^yuM  =  0  , 


which  yields,  upon  contracting  with  ua  , 


(yu6) i3  =  0  , 


since  u  ua  =  -c2  and  consequently  ua i nu  =  0  .  Thus, 
a  ^  J  6  a  5 


provided  y  ^  0  , 


n  a  6 
0  =  u  | ^u 


( ~  a  ,  Ys  3 

(3bu  fByu  )u 


d2xa/dx2  +  r“  (dx®/dx ) (dxY/dx )  , 


which  shows  that  the  dust  particles  follow  geodesics 


* 


The  remaining  part  of  this  chapter  is  devoted  to 
a  review  of  the  salient  features  of  the  Schwarzschild 
geometry:  trajectories  of  test  particles  and  photons 

(Sect.  2)  and  the  Kruskal  completion  of  the  Schwarzschild 
manifold  (Sect.  3) • 

We  shall  employ  "geometrized  units"  (c=G=l) 
throughout  the  remainder  of  the  thesis.  In  this  system, 
mass,  time,  and  length  all  have  the  same  dimension. 


' 
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2 .  The  Schwarzschild  Solution.  Motions  of  Test  Particles 
and  Photons. 

The  best  known  non-trivial  solution  of  the 
Einstein's  gravitational  field  equations  is  the  Schwarzschild 
solution.  Its  traditional  form  is^^ 

ds2  =  (1  -  2m/r)-1  dr2  +  r2(d02  +  sin20  dcj)2)  -  (1  -  2m/r)dt2. 

(21) 

The  line-element  is  static  in  the  sense  that  the  coeffi¬ 
cients  are  independent  of  time  t  and  it  is  invariant 
under  time  inversion  (t-*-t).  The  Schwarzschild  radial  co¬ 
ordinate  r  has  the  geometrically  invariant  meaning  that 
a  2-space  S:  r,t  =  const  has  an  area  4-rrr2  (independ¬ 
ently  of  t),  and  (0,<J>)  are  polar  coordinates  on  S.  This 
2-space  is  intrinsically  indistinguishable  from  an  ordinary 
spherical  surface  of  radius  r.  Thus  the  line-element  is 
also  spherically  symmetric.  At  large  spatial  distances 
(r  >>  2m),  it  is  approximately  Minkowskian,  and  by  requir¬ 
ing  that  Newton's  law  should  hold  at  such  a  region  of  weak 
field,  we  can  identify  the  constant  m  with  the  total  mass 
of  the  body  producing  the  Schwarzschild  field.  (Cf.  the 
linearized  line-element  written  down  previously,  eq.(15))- 
Accordingly  we  assume  m  >  0  . 

The  Schwarzschild  solution  satisfies  the  vacuum 
field  equations  (16).  It  owes  its  great  importance  to  a 


'>c  61 
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( 12 ) 

theorem  of  Birkhoff  ,  according  to  which  it  represents 
the  external  vacuum  field  of  any  spherically  symmetric 
uncharged  (static  or  non-static)  distribution  of  matter. 

On  this  solution  are  based  two  of  the  three  famous  experi¬ 
mental  tests  of  the  theory  of  general  relativity.  (However, 
these  do  not  check  the  full  range  of  validity  of  (21), 
since  the  ratio  2m/r  <<  1  outside  the  sun) .  Recently  it 
has  regained  attention  in  connection  with  the  study  of 
relativistic  gravitational  collapse. 


One  always  gets  a  firmer  grasp  of  the  geometry 
of  any  solution  if  one  knows  how  particles  and  light  rays 
behave  in  it.  For  the  Schwarzschild  geometry,  investiga¬ 
tions  on  the  motions  of  test  particles  and  photons  have 

( 14 ) 

been  thoroughly  carried  out  in  the  past  .  We  shall 
review  some  of  the  results  obtained,  as  they  not  only  il¬ 
lustrate  the  qualitative  differences  from  corresponding 
Newtonian  predictions  that  can  arise  when  the  geometry  of 
space-time  deviates  strongly  from  flatness,  but  are  also 
relevant  for  an  understanding  of  certain  features  of  gravi¬ 
tational  collapse. 


Within  the  framework  of  Newtonian  theory,  the 
motion  of  a  particle  in  the  gravitational  field  of  a  mass 
center  is  well-known.  If  the  angular  momentum  is  not  zero, 
the  path  of  the  particle  is  either  an  ellipse  or  a  hyper¬ 
bola,  depending  on  the  total  energy.  Gravitational  capture 


Ao- 


18 


Is  Impossible  except  for  radial  motions. 

Since  the  paths  of  test  particles  in  general 
relativity  are  time-like  geodesics,  they  are  in  principle 
completely  determined  once  the  metric  is  known.  For  the 
Schwarzschild  field,  the  geodesic  equations  (19)  possess 
the  following  first  integrals: 

(dr/dx)2  =  E2  -  U(r,£)  ,  U(r,£)  =  ( l-2m/r ) ( l+£2/r 2 ) ,  (22) 

r2(dtj)/dx)  =  £  =  const  ,  (23) 

and  without  loss  of  generality  0  can  be  set  equal  to 
tt/2  .  The  affine  parameter  x  is  the  proper  time  of  the 
particle,  (dx2  =  -ds2)  ;  the  constant  E  is  the  "total 

energy"  (including  the  rest  mass  energy  y  =  1),  and  the 
constant  £  the  "angular  momentum". 

Fig.  1  (p.26)  is  a  schematic  plot  of  U(r,£)  vs.  r  for 
three  typical  values  of  £  .  One  gets  a  qualitative  idea 
of  how  the  radius  of  a  moving  particle  varies  with  time  from 
an  examination  of  eq.(22)  and  fig.  1.  For  £  >  4m  (i.e., 
4myG/c),  in  which  case  U  has  a  maximum  value  U  >1, 
the  situation  is  as  follows. 

(1)  E  <  1.  The  particle  neither  falls  to  the  center 
r  =  0  nor  flies  off  to  spatial  infinity  ("finite  motion"). 

The  radius  of  closest  approach  is  r  .  >  4m.  Unlike  the 

Newtonian  case,  here  the  orbit  is  in  general  not  a  closed 
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curve.  (This  is  manifested  in  the  secular  shift  of  the 
perihelion  of  Mercury). 

(2)  1  <  E  <  U  .  The  particle,  coming  from  infinity, 

max 

goes  off  again  to  infinity  in  analogy  to  the  hyperbolic 

motion  in  the  Newtonian  case.  r  .  >  3m  . 

mm 

(3)  E  >  U  .  The  particle,  coming  from  infinity, 

max 

spirals  pass  r  =  2m  and  is  eventually  captured. 

For  &  <  4m  ,  all  particles  coming  from  infinity 
(E  >  1)  are  captured;  with  E  <  1  and  /3  •  2m  <  £  <  4m, 
the  particles  either  are  captured  or  execute  finite  motions 
^rmin  >  5  depending  on  E  and  £.  However,  when 

£  <  /3  *2m,  there  exists  no  finite  motion,  and  a  particle 
inevitably  ends  up  at  r  =  0. 

We  now  consider  light  rays.  According  to 
Newtonian  theory,  their  motions  are,  of  course,  uniform 
along  straight  lines.  In  general  relativity,  they  are 
represented  by  null  geodesics,  and  in  the  present  case  can 
easily  be  obtained  from  (22),  (23)  by  a  limiting  process. 

Each  light  ray  may  be  characterized  by  its  impact 
distance  b  at  infinity.  For  b  >  3/3  m,  the  light  ray 
is  deflected  from  its  linear  motion,  but  never  comes  closer 
than  r  =  3m,  and  eventually  goes  off  to  spatial  infinity. 
(For  a  ray  grazing  the  surface  of  the  sun,  the  deflection 
amounts  to  1.75".  This  prediction  of  Einstein  was,  as  is 
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universally  known,  brilliantly  confirmed  during  the  time 
of  total  solar  eclipse  of  1919) .  For  b  <  3/3  m,  all 
photons  are  gravitationally  captured.  A  photon  can  move 
in  a  circular  orbit  at  r  =  3m,  but  that  orbit  is  unstable. 

Finally,  we  mention  that  a  light  ray  emitted  by 
a  stationary  source  located  at  a  radius  r  =  const  >  2m 
can  escape  (and  hence  be  observed  by  some  observer  at 
spatial  infinity)  only  if  its  inclination  to  the  radial 
direction  is  less  than  a  critical  value  To(r)  determined 

by 


sin2  T0  =  27  m2(r-2m)/rJ  ,  To(2m)  =  0  .  (24) 

For  r  >  3m,  so  that  To  >  tt/2  ,  all  outward  going  rays 
escape  to  infinity.  As  r  decreases  from  3m,  the  cone 
with  half-angle  To  begins  to  narrow.  In  the  limit 
r  -*  2m,  only  rays  emitted  in  the  radial  directions  can 
escape;  all  other  emitted  rays,  though  some  of  them  may  go 
a  long  way  out,  are  ultimately  "pulled  back"  by  the  intense 
gravitational  field. 

3 •  The  Schwar zschild-Kruskal  Manifold. 

The  beautiful  simplicity  of  the  Schwarzschild 
metric  (21)  is  in  a  way  deceptive,  for  it  exhibits  patho¬ 
logical  features  that  have  troubled  physicists  since  the 
early  days  of  general  relativity.  Clearly,  (21)  breaks 
down  at  r  =  0  and  r  =  2m.  The  latter,  known  as  the 


ora 
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"Schwarzschild  surface",  is  of  great  interest.  Its  nature 
had  not  been  fully  understood  until  quite  recently.  When 
r  ->  2m,  the  metric  components  g  -*  °°  and  g  +0,  but 

11  4  4 

det(ga^)  remains  finite.  For  lower  values  of  r,  the 
curves  r  =  const  lose  their  time-like  character,  while 
the  curves  t  =  const  are  no  longer  space-like.  But  per¬ 
haps  the  Schwarzschild  surface  is  not  physically  realiza- 

(22) 

ble?  Thus  it  can  be  shown  ,  for  example,  that  no 
static  fluid  sphere  with  any  reasonable  equation  of  state 

Q 

can  have  a  radius  r  less  than  -g-(2m),  which  means  that 
for  this  broad  class  of  sources,  the  Schwarzschild  surface 
is  actually  buried  inside  the  matter  where,  of  course,  the 
vacuum  field  equations  no  longer  hold.  Nevertheless,  one 
may  wish  to  explore  the  Schwarzschild  solution  from  a 
geometrical  point  of  view  when  there  is  no  "real  mass". 

Nor  are  physical  motivations  in  fact  lacking.  For  there 
are  dynamical  objects  which,  starting  from  large  radii, 
continuously  contract  until  they  actually  cross  the 
Schwarzschild  surface  (gravitational  collapse) .  A  full 
study  of  the  dynamics  requires  an  understanding  of  the  geo¬ 
metry  on  and  inside  the  Schwarzschild  surface. 

Modern  advances  in  this  area  reveal  the 
Schwarzschild  surface  as  a  regular  surface,  though  possess¬ 
ing  the  most  remarkable  physical  properties.  The  apparently 
singular  behaviour  of  the  metric  at  r  =  2m  merely  re¬ 
flects  a  poor  choice  of  coordinates.  This  can  be  demonstra- 
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ted  by  embedding  the  entire  Schwarzschild  geometry  in 

(  23 ) 

higher  dimensional  pseudo-Euclidean  spaces  ,  or  by 
exhibiting  suitable  coordinate  transformations  ^"*".  On 
the  other  hand,  the  singularity  at  r  =  0  is  intrinsic 
(i.e.,  coordinate-independent),  as  can  be  seen,  for 
example,  from  the  expression  for  the  invariant  scalar 

R  q  x~Ra^y6  =  48m2/r6  .  (25) 

apyo 

A  maximal  analytic  extension  of  the 

(  9 ) 

Schwarzschild  space-time  has  been  discovered  by  Kruskal 
and  independently  by  Szekeres^^  .  In  terms  of  "null 
coordinates"  (u,v),  related  to  (r,t)  by 


uv  =  ( r/2m-l ) exp ( r/2m)  , 

u/v  =  exp(t/2m)  , 


(26) 


metric  (21)  takes  on  the  form 

ds2  =  (32m3/r) exp(-r/2m)dudv  +  r2dfi2  ,  (27) 

d$72  E  d02  +  sin20dcj)2  ,  (28) 


which  is  manifestly  regular  for  all  r  >  0,  i.e., 

for  °°  >  uv  >  -1.  Since  the  transformation  (26)  is  well- 

defined  for  r  >  2m,  and  G  n  is  a  tensor,  it  follows 

}  aB 

that  G  calculated  from  (27),  vanishes  when  r  >  2m. 

aB 

But  G  n  is  analytic  in  the  domain  of  regularity  of  (27). 
aB 

Hence  (27)  is  a  solution  of  the  vacuum  field  equations 


for  all  r  >  0  . 
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Figure  2  (p.26)  is  a  sketch  of  a  section 
(05cj)  =  const)  of  the  Schwarzschild-Kruskal  manifold.  Lines 
of  constant  r  are  hyperbolas  in  the  (u,v)  map,  and  lines 
of  constant  t  are  straight  lines  passing  through  the 
origin.  However,  the  Schwarzschild  time  t  runs  off  to 
infinity  on  the  coordinate  axes  (u  =  0  or  v  =  0)  and 
is  undefined  for  uv  <  0  (quadrants  II  and  IV,  fig.  2). 

The  pathology  of  the  Schwarzschild  metric  (21)  is  intimate¬ 
ly  connected  with  this  inadequacy  of  the  Schwarzschild 
time  coordinate.  The  extended  manifold  is  maximal  in  the 
sense  that  any  geodesic  of  (27)  is  extendable  to  infinite 
value  of  its  affine  parameter  unless  it  runs  into  the 
intrinsic  singularity  at  r  =  0. 

One  sees  at  once  from  (27)  that  lines  (u,0,4>)  =  const 
and  (v,0,cj))  =  const  are  null  (ds2  =  0).  Hence  the  Kruskal 
null  coordinates  have  the  important  property  that  the  co¬ 
ordinate  lines  represent  the  paths  of  radial  light  rays. 

The  lines  u  =  const  are  radial  incoming  light  rays  in  the 
sense  that  dr/dt  <  0  whenever  t  is  defined.  Similarly, 
v  =  const  represent  outgoing  radial  light  rays. 

The  extended  Schwarzschild  manifold  manifests  a 
curious  duplicity,  in  that  each  of  quadrants  I  and  III 
(fig.  2)  corresponds  to  the  entire  region  (2m  <  r  <  °°  , 

_ oo  <  t  <  °°)  where  the  Schwarzschild  metric  (21)  is 
regular.  Moreover,  the  metric  (27)  and  equations  (26)  are 


s 
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invariant  under  the  transformation  ("reflection") 

(u,v)  -*  (-u,-v).  Consequently  the  idea  of  identifying  the 

points  (u,v,0,<f>)  with  ( -u , -v  ,  0  ,  cf> )  have  been  proposed  and 
(24) 

examined  .  This  freedom  of  identifying  points  without 
disturbing  the  analyticity  of  the  metric  (aside  from  a 
sort  of  "conical"  singularity  at  the  origin  u,v  =  0)  is 
allowed  because  the  gravitational  field  equations  only 
determine  the  local  geometry.  The  resulting  topology  of 
the  manifold,  however,  is  such  that  it  is  impossible  to 
maintain  a  global  distinction  between  past  and  future  in 
the  region  r  <  2m.  Analyses  of  the  paths  of  radial  parti¬ 
cles  show  the  existence  of  self-intersecting  time-like 
lines,  thus  suggesting  the  annoying  possibility  of  causal 
violations . 


The  conventional  interpretation 


(25) 


of  the  top¬ 


ology  of  the  Schwarzschild-Kruskal  manifold,  in  contrast 
to  the  "elliptic  interpretation"  just  mentioned,  assumes 
different  points  of  the  Kruskal  map  to  be  physically 
distinct,  but  imposes  on  the  manifold  a  global  time-direc¬ 
tion.  This  is  achieved  by  introducing,  at  the  expense  of 
the  reflectional  symmetry  intrinsic  to  the  manifold,  a 
time  coordinate  T  and  an  orthogonal  space  coordinate  X 
(see  fig.  2),  where 

T  =  (u-v)/2  , 


X  =  (u+v)/2 


3 


(29) 


’  i 
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in  terms  of  which  metric  (27)  becomes 

ds2  =  ( 32m3/r ) exp ( -r/2m) ( dX2-dT2 )  +  r2dfi2  ,  (30) 

with 

( r/2m-l ) exp (r/2m)  =  X2  -  T2  .  (31) 

The  conventional  interpretation,  though  appar¬ 
ently  accepted  by  most  relativists,  is  not  free  of 

(26) 

objections.  For  a  critique,  see  . 

Many  of  the  features  of  the  Schwarzschild- 
Kruskal  manifold  manifest  themselves  in  the  phenomenon  of 
relativistic  gravitational  collapse,  and  will  be  further 
discussed  in  that  context  (Chapter  III). 


U(rsA) 
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ii 


r/’2m 


Pig.  1.  Effective  radial  potential  curves  [cf.  eq.(22)]for 
different  angular  momenta  i.  (t  <  2 /3m  <  £  <  4m  <  £  ) . 

&  3  2! 


Pig.  2  .  Two  dimensional  subspaces  0,(j)=const  of  uhe 


Schwarzs chi ld-Kruskal  manifold. 
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CHAPTER  III 

GRAVITATIONAL  COLLAPSE 


1 .  Introduction 

Any  collection  of  matter  ("star")  will  contract 
under  its  own  weight  unless  prevented  by  non-gravitational 
forces  or  rotation,.  The  question  is  whether  the  balance 
between  the  outward  push  of  pressure  and  the  inward  pull 
of  gravitation  may  not  be  eventually  upset  if  the  mass, 
and  hence  the  resulting  gravitational  force,  of  a  system 
is  increased  indefinitely .  In  other  words,  given  that  a 
collection  of  matter  of  a  specified  nature  is  in  a  state 
of  equilibrium,  can  we  infer  that  its  mass  must  therefore 
be  less  than  a  certain  finite  value? 

(27) 

According  to  Chandrasekhar  ,  who  uses  the 
Newtonian  equations  of  hydrostatic  equilibrium,  there  is 
indeed  an  upper  limit  to  the  mass  of  a  cold  star  (white 
dwarf)  maintained  in  equilibrium  by  the  pressure  of  a 
degenerate  electron  gas,  of  the  order  of  one  solar  mass, 
Within  the  framework  of  general  relativity,  Oppenheimer  and 
Volkoffv  ;  have  shown  that  for  systems  of  degenerate 
neutron  gas,  the  maximum  mass  for  equilibrium  is  also  of 
the  order  of  one  solar  mass.  Recently,  Wheeler  and  his 
coworkers^^  examined  exhaustively  the  case  of  a  static 
sphere  subject  to  an  equation  of  state  thought  to  repre- 
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sent  matter  at  the  end  of  thermonuclear  evolution.  Again 
there  exists  a  critical  mass  of  the  same  order,  above 
which  there  is  no  equilibrium  configuration. 

Thus  when  a  star  which  in  the  course  of  its 
evolution  has  sufficiently  cooled  down  possesses  a  mass 
appreciably  larger  than  the  critical  value  for  equilibrium, 
it  will  start  to  collapse.  The  study  of  relativistic 
gravitational  collapse  was  pioneered  by  Oppenheimer  and  his 
student  Snyder  in  1939V  ,  but  until  about  five  years  ago, 

little  further  progress  had  been  made.  The  recent  dis¬ 
covery  of  astronomical  objects  known  as  "quasi-st ellar" 
objects^)  has  stimulated  new  interest  in  the  problem,  and 
marked  the  first  time  that  general  relativistic  effects 
had  been  taken  seriously  by  astrophysicists.  However,  for 
many  relativists  the  fascination  of  the  collapse  problem 
lies  primarily  in  the  weird  consequences  which  general 
relativity  appears  to  predict. 

Consider  a  spherical  star  collapsing  adiabati- 
cally  (i.e.,  no  energy  transfer).  Its  interior  geometry 
is  described  by  a  certain  solution  to  the  Einstein's  field 
equations  depending,  of  course,  on  the  specific  energy- 
momentum  tensor  and  equation  of  state  assumed.  On  the 
other  hand,  by  Birkhoff  theorem,  its  exterior  geometry  is 
uniquely  described  by  the  Schwarzschild-Kruskal  line-element. 
Hence  the  history  of  (an  element  of)  the  star's  surface 
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must  be  a  time-like  world  line  of  the  Schwarzschild-Kruskal 
manifold.  For  a  star  which  collapses  to  a  radius  smaller 
than  its  Schwarzschild  radius  r  =  2m  (where  m  is  the 
mass  of  the  star) ,  such  a  world  line  may  be  qualitatively 
represented  by  the  curve  labelled  p  =  p^  in  fig.  4  (p.43c), 
with  the  portion  to  the  left  of  the  curve  representing  the 
exterior  region.  From  the  dispositions  of  the  future  null 
cones  shown  in  the  figure,  it  is  easy  to  see  that,  once 
the  star  has  collapsed  past  the  Schwarzschild  surface 
(v  =  0,  fig.  4),  curious  things  happen:  (1)  no  signals  or 
particles  emitted  on  the  star’s  surface  can  ever  escape  to 
the  external  world  (r  >  2m),  and  (2)  the  collapse  cannot 
be  stopped  but  must  rather  proceed  inexorably  until  the 
entire  star  gets  "swallowed  up”  at  the  singularity  r  =  0. 
Paradoxically,  to  an  external  observer  who  stays  out  of 
the  enigmatic  Schwarzschild  surface,  the  star  never  shrinks 
to  a  radius  less  than  2m  at  all.  This  we  conclude  from 
the  fact  that  the  world  line  of  the  star's  surface  cuts  the 
Schwarzschild  surface  at  time  t  =  +°° ,  where  t  is  the 
Schwarzschild  coordinate  time,  essentially  the  time  measured 
by  our  observer.  Thus  what  he  sees  is  that  the  star  begins 
to  slow  down  as  the  radius  r  =  2m  is  approached,  reach¬ 
ing  it  only  asymptotically. 

We  shall  return  to  a  more  detailed  discussion 
of  these  peculiar  features  of  collapse  and  the  important 
question  of  their  generality  in  Sects.  3  and  4,  turning  our 
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attention  now  to  the  interior  of  the  collapsing  star  and 
the  precise  world  lines  followed  by  its  particles.  The 
collapse  of  a  fluid  sphere  of  uniform  density  and  zero 
pressure  ("homogeneous  ball  of  dust"),  which  is  instruc¬ 
tive  because  of  its  simplicity,  has  been  thoroughly 

(  g  ) 

analyzed  .  Our  aim  here  is  to  render  the  results  more 
physically  transparent  by  introducing  a  unified  system  of 
null  coordinates  which  covers  both  the  interior  and  the 
exterior  domains  of  the  collapsing  object  (Sect. 2). 

2 .  Null  Coordinates  for  Spherically  Symmetric  Gravitational 

* 

Collapse . 


For  the  interior  domain  of  a  homogeneous  ball 
of  dust,  the  metric,  as  determined  by  the  gravitational 
field  equations  and  the  symmetry  of  the  problem,  turns  out 
to  be^^  identical  with  that  of  a  Friedmann  universe: 

(ds2)j-nt  =  S2 ( x ) { ( 1-kr 2 ) ~ 1  dr2  +  r2  d^2}  -  dx2,  ( k=0  ,  ±1 ) 

(1) 

dfi2  =  d62  +  sin20d(})2  ,  (2) 

where  S(x),  which  essentially  measures  the  radius  (=rS) 
of  the  ball,  is  determined  by 

(dS/dx)2  =  6/S  -  k  ,  ( 3 


*  Based  on  V.  de  la  Cruz,  reference  33- 
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3  being  a  positive  constant.  Here  t  measures  the  proper 
time  along  the  world  lines  of  the  dust  particles,  and 
r,  0,  <p  are  co-moving  coordinates,  i.e.,  the  particles  are 
at  rest  in  this  coordinate  system,  so  that  the  four-velocity 
of  the  dust  has  the  form  ua  =  6a  .  The  only  non- 

4 

vanishing  component  of  the  dust  energy-momentum  tensor 

Ot  06 

T  ^  =  yu  Ug  gives  the  energy  density 

-  T4  =  y 

4 

=  33/ (  8ttS  3 )  .  (4) 

From  eq.(3),  we  easily  see  that  the  constant  k  has  the 
following  significance,  as  well  as  the  usual  interpretation 
in  term  of  the  curvature  of  the  spatial  section  x  =  const, 
k  =  0  :  collapse  from  a  state  of  rest  at  infinity;  k  =  +1  : 
collapse  from  a  state  of  rest  at  a  finite  radius;  k  =  -1  : 
collapse  from  infinity  with  non-zero  initial  velocity. 

The  exterior  metric  of  the  ball  of  dust,  comple¬ 
menting  (1),  is  simply  the  Schwarzschild  metric,  eq.(II.21): 

(ds2)Ext  =  (l-2m/r)_1  dr2  +  rz  dft2  -  (l-2m/r)dtz  ,  (5) 

where  m  is  the  total  mass  of  the  dust. 

To  get  a  unified  picture,  the  interior  co¬ 
moving  coordinate  system  can  be  extended  to  the  exterior 


Greek  indices  take  values  1  -  4  ;  latin,  1-3- 
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region  ;  however,  it  retains  no  direct  geometrical 

significance  there,,  Extensions  of  Schwarzschild ' s  co- 

(  31 ) 

ordinates  to  the  interior  break  down  when  the  collapse 

has  proceeded  to  the  Schwarzschild  radius  r  =  2m. 

A  system  of  null  coordinates  of  Kruskal's  type 
covering  both  interior  and  exterior  of  a  ball  of  collap¬ 
sing  dust  will  now  be  given.  It  is  defined  down  to  the 
singular  event  of  maximum  implosion,  and  gives  a  direct 
picture  of  the  paths  of  radial  light  rays  and  the  disposi¬ 
tions  of  null  cones  throughout  the  manifold.  For  simpli¬ 
city,  attention  is  confined  to  the  special  solution  of 
Oppenheimer  and  Snyder  (k  =  0),  but  the  analysis  can  be 
readily  extended  to  the  other  cases  (k  =  ±  1),  with  quali¬ 
tatively  similar  results. 

For  k  =  0,  eq.  (3)  yields,  after  a  simple 
integration , 

S3/2  =  -  |  3"2  t  +  const  ,  (6) 

where  the  constant  of  integration  will  be  set  to  zero  by 
adjusting  the  zero  point  of  x.  Substituting  (6)  into 

(32) 

(1),  we  obtain  the  interior  metric  for  the  case  k  =  0: 

(ds2)int  =  t4/3  (dp2  +  p2dfi2)  -  dx2  (7) 

in  which 


P 


(jj  6)  '/3  ? 


(8) 
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is  also  a  co-moving  radial  coordinate.  Eqs .  (4)  and  (6) 
yield  the  expression  for  the  homogeneous  density  as  a 
function  of  time. 


y  =  1/(6ttt2)  . 


(9) 


Thus  (in  summary)  the  line-element  (7)  represents,  in 
terms  of  co-moving  coordinates,  the  field  in  the  interior 
of  a  homogeneous  ball  of  dust,  with  density  given  by  (9), 
collapsing  (for  x  <  0)  from  a  state  of  rest  at  infinity. 

We  will  first  seek  to  express  the  line-element 
(7)  in  terms  of  the  Schwarzschild  radial  coordinate 


r  =  px2/3 


(10) 


(the  right  hand  side  expression  being  the  square  root  of 
the  coefficient  of  dft2  in  eq . ( 7 ) )  ,  and  an  orthogonal 
time  coordinate  t  =  t(p,x).  Using  the  general  tensorial 
transformation  equations 


(ID 


we  find 


( 12A ) 
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where  use  has  been  made  of  (10)  and  (7), 


gtt  =  -  (t,x)2  +  T  4/3  (t,p)2  , 


(12B) 


and 


rt  n 
g  =  0  = 


-  |  px"1/3  t,  +  x  2/3  t, 

3  x  sp 


( 12C ) 


A  particular  solution  of  (12C)  is 


t  =  p2  +  \  x2/3  ,  (13) 

with  which  can  be  calculated  from  (12B).  Finally, 

the  line-element  (7)  takes  the  form 


(ds2)int  =  (l-2M(p)/r)  1  dr2  -  g  x 2/ 3 ( 1-2M( p ) /r )  1  dt 2  + 

+  r2  dft2  ,  (14) 

where 

M(p)  e  Trr  3y  ( x )  =  |  p3  (15) 

is  the  "mass”  of  a  sphere  interior  to  a  given  dust 
particle  with  co-moving  coordinate  p.  The  coefficients 
in  (14)  could  be  written  as  functions  of  the  new  coordi¬ 
nates  r  and  t  by  using  (10)  and  (13). 


Let 


i  .  e 


p  =  pp  =  const  , 


r  =  rp(x)  =  Ppx2/3 


(16) 

(17) 
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represent  the  boundary  of  the  ball  of  dust,  at  which 
(14)  must  join  continuously  to  the  exterior  metric  (5). 

The  condition  is  automatically  satisfied  by  the  coeffi¬ 
cients  of  dfi2.  For  those  of  dr2  to  match,  M(p^)  must  be 
identified  with  the  total  mass  of  the  dust  ball,  i.e., 

m  =  M(pb)  =  |  pb  .  (18) 


Finally,  one  requires 


|  t2/3  (l-2m/rb)  1  dt2  =  (l-2m/rb)dt2 


for  p  =  pb,  i.e. , 


dt/dt  =  i  (-t)1/3  (l-2m/rb)  1 


=  3  [|(t-pb2)]3/2[|(t-pb2) 


I  P  2]-i 
9  pb  J  5 


(19) 


the  last  equation  being  obtained  with  the  help  of  (13), 
(16),  and  (17).  Integration  of  (19)  yields 

t  “  -T3  -  3  Pb2T  +  g  Pb3  lnC(T  +  §  pb)/(T  '  1  pb)]  >(2Q) 


where 

T(t)  E  [  |  (t-pb2)  ]h  (21) 

and  the  constant  of  integration  has  been  set  to  zero. 


In  this  way,  Schwarzschild  coordinates  (r,0,<j>,t) 
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are  extended  continuously  into  the  interior  -  the 

resulting  metric  is  (14),  with  t  =  t(t)  defined  implicitly 
by  (20)  and  (21).  However,  (20)  shows  that  the  Schwarzschild 
time  coordinate  becomes  infinite  (and  useless)  when 

^|^-pb2 )  ±  ?  pb  =  0  , 


l .  e . 


p2  +  \  T2/3  = 


3  Pb2 


(22) 


For  the  boundary  p  =  pfe ,  this  happens  when 


Z  t2/3  =  2  d  2 
2  pb  * 


i  .  e  .  , 


rb  =  Pb1'7' 


4  3 

=  9  Pb 


=  2m  , 


(23) 


that  is,  when  the  radius  of  the  ball  is  equal  to  its 
Schwarzschild  radius,  as  should  be  expected.  Thus  some 
other  system  of  coordinates,  such  as  null  coordinates  of 
Kruskal’s  type,  will  be  needed  to  cover  the  entire  mani¬ 
fold. 


Recall  (Chap. II)  that  the  exterior  metric, 
maximally  extended,  can  be  expressed  in  terms  of  Kruskal 
null  coordinates  (u,v)  by 


, 

.•  \U:C 
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(ds2)^xt  =  (32m3/r)  exp  (~r/2m)dudv  +  r2dft2  ,  (24) 

(u,v)  being  related  to  (r,t)  by 

u2  =  (r/2m-l)  exp  (r+t)/2m  ,  (25) 

v2  =  (r/2m-l)  exp  (r-t)/2m  .  (26) 

The  lines  u  =  const  and  v  =  const  represent  the  paths 
of  incoming  and  outgoing  radial  light  rays,  respectively. 


we  set 


and 


To  find  the  radial  null  lines  in  the  interior, 
dfi2  =0  in  ( 7 )  and  solve  ds2  =  0  ,  obtaining 

1/3 

x  =  (3t  +  p-pb)/2pb  =  const  (incoming)  ,  (27) 


y  =  (3T1/3-  p+pb)/2pb  =  const  (outgoing)  .  (28) 


Thus  (x,y)  are  null  coordinates  for  the  interior,  analy¬ 
tically  related  to  p  and  t  everywhere  except  at  the 
geometrical  singularity  (instant  of  maximum  implosion) 
t  =  0.  Since  0  <  p  <  pfe  and  t  <  0  ,  we  have  x  <  0  and 
y  <  \  .  Along  the  curves  p  =  const,  y  extends  to  ^(l-p/p^) 
before  reaching  maximum  implosion.  See  fig.  3  (p.43a). 


Now  an  external,  say  incoming  radial  null  line 
u  =  u  =  const  (for  fixed  0 ,  c/> )  is  uniquely  determined  by 
the  event  (r  ,t  )  at  which  it  meets  the  boundary,  u  being 

i  i  i 

found  by  substituting  r  =  r  ,  t  =  t  into  (25).  In  the 

i  i 

interior,  we  now  seek  an  analytic,  monotonic  function 
u  =  u(x)  such  that  u(x)  =  u  represents  the  continuation 


'  • 
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of  the  external  null  line  u  =  u  ,  i.e.  such  that 

1 

u  =  const  labels  the  same  incoming  light  rays  throughout 
the  space-time  manifold. 


Let  the  event  (r  ,t  )  be  also  specified  in  terms 
of  the  co-moving  coordinates  by  (p^,T  ).  Then,  from  (17) 
and  ( 20 ) , 


and 


r 

l 


(29) 


t 

i 


*5 


Pi 


1/3 


4 

9 


In  [  (  t 


!/3  _ 


2 

3 


Pb)/(Ti 


1/3. 


'3pb 


)] 


(30) 


Substituting  these  into  (25)  and  making  use  of  (27),  we 
obtain 

u2  =  (r  /2m-l)  exp  [(r  +  t  )/2m] 
ii  ii 

=  ( x ^  — 1)  exp  (2X[+  x/|x  3)  ,  (31) 

where 


x,  =  x(pb>T1)  • 

Hence  the  desired  function  u  =  u(x)  is 

u  =  (x-1)  exp  (x+ix2+ix3)  ,  (x  <  0)  . 


(32) 


(33) 
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Similarly,  the  other  function  v  =  v(x)  can  be  found  to 
be 

v  =  vi  =  (y+i)  exp  (-y+^y2-^y3)  >  (y  <  o)  .  (34) 

Eqs.  (33)  and  (3*0  together  define  a  coordinate  chart 
(chart  I),  in  terms  of  which  the  interior  metric  becomes 

(ds2)I  =  m2 (x+y) 4 (xy)“3  dudv  x 

x  exp{-[i(x3~y 3)  +  i  (x2-y2)  +  (x-y)]}  +  r2dft2  .  (35) 

At  the  boundary  p  =  p^  ,  x  =  y  ,  this  joins  continuously 
to  the  exterior  Kruskal  solution  (24). 

Since  v(y)  ceases  to  be  monotonic  at 
y  =  0,  the  domain  of  validity  of  I  does  not  include 
the  whole  of  the  interior.  Another  chart  is  needed  to 
cover  the  domain  0<y<^,  i.e.,  the  region  AOB 

in  fig.  3-  Define  chart  II  by 

u  =  (x-1)  exp  (x+^x2+^x3)  ,  (x  <  0)  ,  (36) 

v  =  v2  =  (y+l)~ 1  exp  (y-^y2+^y3)  ,  (0  <  y  <  \)  . 

(37) 
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Then  the  interior  metric  in  terms  of  chart  II  is 
(dsI 2)n  =  -m2 ( 1+y ) 2 ( xy ) ~ 3 ( x+y ) 4  dudv  x 

x  exp {-[^ ( x  3+y  3 )  +  i(x2-y 2 )  +  (x+y)]}  +  r2dft2  . 

(38) 

The  coordinates  v  defined  by  charts  I  and  II  match 
continuously  (though  not  analytically)  at  y  =  0  ,  so  that 
no  concusion  could  arise  from  using  the  same  labels  for 

the  coordinates.  Moreover,  the  lines  v  =  const  for 

% 

0  <  y  <  h  (for  1  <  v  <  1.0013)  do  not  reach  the  bound¬ 
ary  p  =  ,  so  there  is  no  question  of  continuity  of  v 

across  this  boundary  to  the  exterior  v  of  eq.(24).  The 
charts  I  and  II,  which  are  adjacent  instead  of  overlapping, 
are  properly  connected  to  each  other  by  the  (x,y)  chart, 
which  overlaps  both. 

Figure  4  is  a  sketch  in  which  we  drew,  in  a 

(u,v)  diagram,  not  exactly  true  to  scale,  the  curves 

p  =  0  and  p  =  Pb  until  they  meet  the  curve  t  =  0 

(uv  =  -1).  The  exterior  curve  uv  =  -1  and  the  interior 
2 

curves  uv  =  -1  and  p  =  0  together  form  the  curve  r  =  0 
2 

of  vanishing  spherical  area  (shaded  along  its  side).  The 
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dotted  line  at  v  =  1  (y  =  0)  separates  the  region  of 
chart  II  from  the  region  of  chart  I.  The  size  of  the  region 
of  chart  II  is  exaggerated  in  the  figure.  Actually,  its 
three  corners  lie  at  (u  =  -0.8692,  v  =  1  ;  x  =  -1,  y  =  0  ; 
p  =  0  ,  t  =  -m/6),  at  (u  =  -0.9889,  v  =  1.0113  ;  -x  =  y  =  0.5; 
p  =  t  =  0)  ,  and  at  (-u=v=  1  ;  x=y=  0  ;  p=pb,  x  =  0). 

The  relation  between  the  x,y  scale  and  the  u,v  charts  are 
indicated  by  the  following  little  table. 
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-1.50 

0 

0 

• 

CM 

1 

-2.50 

u 

-1.0000 

-0.9889 

-0.8692 

-0.5578 

-0.2085 

-0.0358 

y 

-2 . 00 

-1.50 

O 

O 

1 — 1 

1 

-0.50 

0.00 

0.50 

V 

-785.77 

-21.2605 

0.00 

0.9739 

1.00 

1.0113 

The  curves  x  =  const  and  p  =  const  (with  the  exception 

of  x  =  0  ,  p  =  0  ,  and  p  =  p^)  are  not  shown  in  fig.  4. 

The  curves  p  =  const  rise  and  bend  over  till  they  are 

tangent  to  v  =  1  (y  =  0).  From  here,  we  continue  these 

curves  for  y  >  0  into  the  region  of  chart  II,  so  that  in 

the  u,v  plane  these  curves  show  an  inflection  point  at 

* 

v  =  1.  The  only  exception  is  the  curve  P  =  Pb  which  meets 
v  =  1  at  u  =  -1  under  45°  (dv/du  =  -1),  and  which  does 
not  continue  beyond  y  =  0.  The  apparent  light-like  direction 
which  all  the  other  curves  p  =  const  instantaneously  seem 


*  The  curve  p  =  p^  is  simultaneously  a  geodesic  of  the 


exterior  and  interior  geometries . 
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to  take  at  v  =  1  (y  =  0)  has  no  physical  reality  and 
may  be  regarded  as  a  local  breakdown  of  charts  I  and  II 
where  they  meet,  as  (fig.  3)  we  always  have  dy/dx  =  1 
along  p  =  const.  Similarly,  always  dy/dx  =  -1  along 
the  curves  t  =  const,  so  that  the  lines  u  =  const  and 
v  =  const  bisect  the  angles  made  by  the  curves  p  =  const 
and  x  =  const. 

We  have  sketched  in  a  few  curves  r  =  const  and 
t  =  const.  The  Schwarzschild  time  t  is  defined,  for  the 
exterior,  by  (25),  (26),  and,  for  the  interior,  by  (20), 
(21).  Again  the  lines  u  =  const  and  v  =  const  bisect 
the  angles  between  the  intersecting  curves  r  =  const  and 
t  =  const. 


As  noted  before,  the  Schwarzschild  chart  breaks 
down,  for  the  exterior,  when  u/v  <  0  ,  and,  for  the 
interior,  when  9t2/3  +  2p2  <  6  pfe2  ,  that  is,  when 
3(x-y)2  <  4(l+y)(l-x).  The  present  null  coordinates 
enable  us  to  give  a  unified  map  (fig.  4)  of  the  complete 
exterior  and  interior  manifolds. 
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Fig.  3-  The  interior  domain  (bounded  by  the  lines  p  =  0, 
p  =  Pb  j  and  t  =  0)  of  a  collapsing  ball  of  dust  in  terms 
of  null  coordinates  (x,y)  [cf.  eqs . ( 27 , 28 ) ] .  Also  cf.  fig.  4. 
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LEGEND  FOR  FIGURE  4 

Unified  map  of  the  complete  exterior  and  interior 
manifolds  of  a  collapsing  ball  of  dust.  The  interior  is 
bounded  by  the  outside  boundary  (the  curve  p  =  p^) ,  the 
history  of  the  center  of  the  sphere  (the  curve  p  =  0) ,  and 
the  instant  of  maximum  implosion  (the  curve  t  =  0) .  On 
t  =  0  as  well  as  on  p  =  0,  we  have  r  =  0  (zero  spherical 
area).  I  (II)  denotes  the  region  of  validity  of  chart  I 
(chart  II).  The  dotted  line  at  v  =  1  separates  the  two 
regions.  The  size  of  the  region  of  chart  II  is  exaggerated 
in  the  figure.  Actually,  its  three  corners  lie  at 
(u  =  -0.8692  ,  v  =  1),  at  (u  =  -O.9889,  v  =  1.0113),  and 
at  (u  =  -1,  v  =  1).  Typical  lines  of  constant  Schwarzschild 
radial  coordinate  r  and  time  coordinate  t  are  shown 
[cf . eq . ( 20 ) ] .  Wavy  lines  represent  radial  signals  emitted 
on  the  surface  of  the  dust  sphere  at  regular  proper  time 
intervals.  Also  shown  are  some  future  null  cones.  Parts 
of  the  figure  are  not  drawn  to  scale  in  order  to  make  some 
of  the  features  of  this  diagram  qualitatively  more  conspi¬ 
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3 •  Features  of  Spherically  Symmetric  Gravitational  Collapse. 

Before  going  on  to  the  question  of  the  generality 
of  the  features  of  collapse  outlined  in  the  introductory 
section,  we  pause  to  describe  these  features  in  greater 
detail,  basing  our  discussion  on  the  dust  model,  which  is 
historically  the  first  example  studied  of  an  object  that 
actually  collapses  past  its  Schwarzschild  radius. 

The  picture  as  seen  by  an  external  stationary 

observer  (world  line  r  =  r  =  const  >  2m) ,  whose  clock 

i 

P 

registers  time  according  to  t  =  ( l-2m/r  ) 2  t  ,  is  as 

i  i 

follows.  At  first,  when  the  stellar  radius  is  still  large, 
the  picture  is  qualitatively  Newtonian,  since  the  distor¬ 
tion  of  geometry  produced  by  the  star  is  then  negligible. 

As  the  star  shrinks  toward  the  Schwarzschild  surface 
r  =  2m,  its  motion,  instead  of  being  accelerated  as  one 
would  expect  from  the  Newtonian  description,  appears  to  an 
external  observer  to  slow  down,  reaching  the  radius  r  =  2m 
only  asymptotically.  At  the  same  time,  the  star  becomes 
redder  and  redder.  Fig.  4  gives  a  qualitative  explanation. 
The  wavy  lines  represent  radial  light  rays  sent  out  at 
regular  proper  time  intervals  by  an  audacious  astrophysicist 
riding  on  the  surface  of  the  star  and  falling  with  it 
("co-moving  observer").  By  this  means,  the  external  obser¬ 
ver  (r  =  5m/2  in  the  figure)  is  informed  of  the  progress 
of  the  collapse.  We  see  that,  as  the  collapse  proceeds 
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toward  the  Schwarzschild  surface,  the  light  rays  are 
received  at  increasingly  longer  intervals;  it  will  take  an 
infinite  time  to  receive  the  signal  sent  out  at  the  instant 
when  the  co-moving  observer  crosses  the  Schwarzschild 
surface,  while  all  later  signals  could  never  reach  the 
external  observer  at  all.  At  the  same  time,  the  signals 
will  also  be  progressively  red-shifted  (partly  because  of 
Doppler  effect,  partly  because  the  photons  will  have  to 
spend  more  and  more  energy  to  overcome  the  gravitational 

pull)  -  infinitely  so  for  the  last  signal  received.  In 

addition  (but  this  we  cannot  deduce  from  fig.  4),  more  and 
more  of  the  light  rays  emitted  at  the  stellar  surface  at 
oblique  angles  will  be  recaptured  by  the  growingly  Intense 
field,  in  a  way  similar  to  the  case  of  light  rays  emitted 
by  stationary  observers  (cf.  Chapter  II).  Thus  the  lumi¬ 
nosity  drops  to  zero  asymptotically. 

Nevertheless,  to  the  co-moving  observer  it  takes 
a  perfectly  finite  time  for  the  star  to  collapse  from  an 
arbitrary  radius  rQ  >  2m  to  the  Schwarzschild  surface  and 
thence  to  the  singular  state  r  =  0.  Indeed,  if  the 
density  corresponding  to  radius  rQ  is  yQ,  then  (9)  gives 
at  once 

Ax q  =  (  6tt y q )  ^  (i.e.,  =  (6TrGyQ)-^)  (39) 


for  the  time  of  collapse  to  the  singularity.  (By  coinci¬ 
dence,  this  is  identical  with  the  prediction  of  Newtonian 
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theory).  For  ^  1  g/cm3  ,  Atq  is  of  the  order  of 

10  sec.  Since  the  proper  time  x  has  an  invariant  mean¬ 
ing  in  general  relativity  (dx2  =  -ds2) ,  this  description 
by  the  co-moving  observer  is  essentially  coordinate-independ¬ 
ent  . 


Thus  the  star  does  cross  the  schwarzschild 
surface  in  a  finite  proper  time  from  any  arbitrary  radius, 

and  in  fact  that  event  will  not  appear  to  be  specially 

* 

remarkable  to  an  innocent  co-moving  observer.  However, 
when  this  has  happened,  gravitational  self-closure  takes 
place:  the  star  loses  all  possibilities  of  communication 

with  the  external  world.  In  other  words,  the  Schwarzschild 
surface  is  an  "event  horizon"  for  all  external  observers. 
Only  the  gravitational  field  of  the  star  will  be  felt.  The 
star  itself  reaches  a  state  of  zero  volume  and  infinite  den¬ 
sity  at  proper  time  x  =  0,  producing  a  singularity  of 
infinite  curvature.  Further  analysis  would  not  seem  to  be 
possible  or  meaningful  .  However,  if  our  ordinary  idea  of 
causality  is  to  be  preserved,  the  star  can  never  re-emerge 
to  the  same  external  world. 


*  The  physical  conditions  on  the  Schwarzschild  surface  are 
in  general  quite  normal,  except  for  the  intense  gravi¬ 
tational  field.  However,  according  to  the  Principle  of 
Equivalence,  this  would  not  give  rise  to  any  local  change 
of  physical  laws. 

**  The  entire  star  probably  gets  "crushed  out  of  existence" 
by  the  infinite  tidal  forces  at  the  singularity  3 u  .  For 
a  different  point  of  view,  cf.  Israel  54  . 
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Whether  these  peculiarities  are  qualitatively 
specific  to  the  present  much-too-idealized  model  is  a  ques¬ 
tion  that  naturally  comes  to  mind.  Certainly  they  appear 
to  be  more  general  than  the  model  itself,  as  the  discussion 
in  the  introductory  section  indicates.  For  a  sufficiently 
massive  star,  neither  the  density  nor  the  pressure  is  im¬ 
moderately  large  when  the  radius  is  of  the  order  r  =  2m, 
and  so  perhaps  the  idealization  of  zero  pressure  is  not 
important.  But  once  the  Schwarzschild  radius  is  passed, 
we  saw  that  the  collapse  must  proceed  irreversibly  to  zero 
volume,  provided  the  assumption  of  spherical  symmetry  is 
not  violated,  and  provided  no  energy-transfer  occurs  through¬ 
out.  That  pressure  forces,  which  are  expected  to  increase 
with  compression,  are  powerless  to  stop  the  collapse  at 
these  late  stages  from  progressing  to  completion  might  appear 
strange.  But  we  have  to  bear  in  mind  the  fact  that  for 
normal  matter  a  larger  pressure  is  always  accompanied  by  a 
larger  internal  energy,  which  contributes  to  the  gravita¬ 
tional  mass  and  hence  to  self-attraction. 

Such  considerations  have  been  confirmed  by  recent 

work.  For  spherically  symmetric  collapse  of  (non-rotating) 

fluid  spheres,  it  is  now  known  that  the  employment  of  more 

( 35 ) 

realistic  physical  conditions  (inhomogeneous  matter  , 

non-vanishing  pressure  and  pressure-gradients1  ,  or  even 

(37) 

the  emission  of  isotropic  radiation  )  generally  pre¬ 
serves  the  qualitative  features  presented  above. 
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The  situation  with  regard  to  non-spherical 
collapse,  with  or  without  rotation,  is  much  less  clear.  The 
problem  is  still  too  complicated  to  be  handled  by  the  analy¬ 
tic  or  numerical  techniques  now  at  our  disposal.  It  has 

/  o  O  \ 

been  claimed^  ;  that  small  departures  from  sphericity  or 

small  rotations  are  likewise  without  qualitative  effects. 

We  note  that  on  the  basis  of  Newtonian  theory,  even  a  small 

( 39 ) 

rotation  could  prevent  unlimited  contraction.  This 

corresponds  to  the  fact  that  a  test  particle  with  an  arbi¬ 
trarily  small  angular  momentum  can  (according  to  Newtonian 
theory)  move  along  a  stable  orbit  around  a  gravitational 
force  center  (of  mass  m,  say).  On  the  other  hand,  as  was 
noted  in  Chap. II,  general  relativity  predicts  that  for  an 

angular  momentum  less  than  2/3m,  the  particle  will  inevi- 

(  39 ) 

tably  be  captured.  In  other  words  ,  the  gravitational 
field  can  be  strong  enough  even  to  crush  rotation,  a  situa¬ 
tion  that  never  occurs  in  the  Newtonian  case.  Heuristic 
arguments  of  this  nature,  however,  though  they  can  be 

/  o  O  \ 

refined^  ,  cannot  be  trusted  completely  .  The  next  section 
deals  with  the  question  of  non-spherical  collapse,  citing 
a  number  of  interesting  works. 

In  Chap.  IV  we  shall  study  the  collapse  of 
charged  bodies. 
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4 .  On  Asymmetric  Gravitational  Collapse  and  Singularities. 

We  have  seen  that  for  spherical  collapse,  an 
event  horizon  generally  develops  which  seals  off  the  star 
from  the  external  world.  The  subsequent  fate  of  the  star 

can  be  followed  until  it  ends  in  a  singularity  characterized 

* 

by  infinite  curvature  and  density. 

The  occurence  of  singularities  is  a  prevalent 
feature  of  the  known  exact  solutions  of  Einstein's  gravita¬ 
tional  field  equations.  Some  of  the  best  known  of  these 
are  the  Schwarzschild  solution,  the  Reissner-Nords trom  solu¬ 
tion^"1")  ,  the  Kerr  solution^  ^  ,  and  the  Friedmann's  cosmo¬ 
logical  solutions.  The  prediction  of  physical  singularities 
by  a  theory  for  possibly  physically  realizable  systems 
(e.g.,  cosmological  models  and  models  for  gravitational 

•j* ** 

collapse)  is  a  matter  of  serious  concern  ,  perhaps  signalling 
the  breakdown  of  the  theory  itself.  However,  all  exact 


*  Presumably  quantum  mechanical  fluctuations  in  the  curva¬ 
ture  of  space-time  must  be  taken  into  account  when  the 
gravitational  field  becomes  extremely  intense.  However, 
these  probably  would  not  become  important  until  a  density 
of  the  order  of  1093  g/cm3  is  reached  40  . 

**  For  opposing  views  on  whether  singularities  can  occur  in 
realistic  cosmological  models,  see  Lifschitz  et  al  44 
and  Hawking47. 

t  "Physicists  abhor  singularities,  and  they  like  to  think 
that  Nature  does  too"  -  K.S.  Thorne. 
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solutions  involve  some  special  assumptions  concerning  the 
metrics  or  the  distributions  of  matter,  which  are  often 
motivated  by  mathematical  convenience  rather  than  physical 
necessity.  To  what  extent,  in  particular,  is  the  assumption 
of  strict  sphericity  responsible  for  the  appearance  of  the 
singularity  in  gravitational  collapse? 

Closely  connected  with  the  concept  of  singularity 

is  the  notion  of  "geodesic  incompleteness."  Although  no 

entirely  satisfactory  general  definition  of  a  singularity 

has  yet  been  given,  a  region  where  the  energy  density  or 

some  invariants  of  the  Riemann  curvature  tensor  (e.g.,  the 

(42) 

Petrov  scalars  )  become  infinite  must  presumably  be 

regarded  as  singular.  However,  there  is  some  difficulty 

in  defining  the  concept  -  needed  for  global  analyses  - 

of  a  singular  space-time  manifold  as  one  containing  singular 

points,  since  a  new  manifold  may  simply  be  defined  without 

these  points.  Nevertheless,  the  two  manifolds  would  share 

a  common  feature:  there  are  geodesics  which  cannot  be 

extended  to  arbitrary  values  of  their  affine  parameters. 

The  three  types  of  geodesic  incompleteness  -  that  of 

time-like,  null,  and  space-like  geodesics,  respectively  - 

(43) 

are  however  not  equivalent  .  There  are  space-times  which 
are,  for  example,  space-like  complete  but  neither  time-like 
nor  null  complete.  But  because  time-like  and  null  geodesics 
represent  world  lines  of  test  particles  and  light  rays, 
time-like  and  null  completeness  would  seem  to  be  a  minimum 
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condition  for  a  space-time  manifold  to  be  regarded  as 
"singularity-free" » 

In  what  follows  we  shall  quote  and  briefly  assess 
the  important  works  of  Lifschitz  and  Khalatnikov,  of 
Penrose,  and  of  Israel,  which  together  form  the  basis  on 
which  present  day  speculations  on  the  qualitative  effects 
of  asymmetry  on  collapse  are  made. 

( A)  Lif schit z-Khalatnikov  solution  near  a  singularity. 

( lili) 

In  1959  and  I960,  Lifschitz  and  Khalatnikov 
sought  to  solve  the  following  problem  which  they  posed: 
Assuming  a  singularity  to  exist,  it  is  required  to  find 
near  it  the  form  of  the  broadest  class  of  solutions  of  the 
gravitational  equations,  so  as  to  judge  from  the  number  of 
arbitrary  functions  (of  the  spatial  coordinates)  it  con¬ 
tains,  whether  this  solution  is  general.  We  remind  our¬ 
selves  that  among  the  arbitrary  functions  contained  In  a 
solution,  generally  there  are  some  whose  arbitrariness  is 
connected  with  the  arbitrariness  of  the  coordinate  system 
itself,  and  are  hence  not  directly  associated  with  the 
physics  of  the  problem  .  The  number  of  physically  arbitrary 
functions  in  the  general  solution  of  the  gravitational 
equations  for  a  perfect  fluid  is  8,  corresponding  to  the 
8  pieces  of  initial  data  which  must  be  specified  to  fix 
the  future  evolution  of  the  system.  (These  are  the  spatial 


distribution  of  the  fluid  density,  the  three  components 
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of  the  fluid  velocity,  and  the  amplitudes  and  rates  of 

(13) 

change  of  amplitudes  for  two  modes  of  gravitational  waves.) 


For  their  purpose,  Lifschitz  and  Khalatnikov  make 

use  of  synchronous  coordinate  systems,  characterized  by 

* 

the  properties  g  =  0  ,  g  =  -1  ,  and  assume  the 

4  i  4  4 

equation  of  state  o f  an  ultra-relativistic  gas  P  = 
near  the  singularity  (at  which  the  pressure  and  density 
supposedly  become  infinite).  Furthermore,  they  confine 
their  explicit  analysis  to  solutions  whose  singularities 
occur  on  a  space-like  hypersurface,  contending  that  this  is 
the  type  of  singularity  which  is  of  major  interest  here. 

By  a  suitable  choice  of  the  synchronous  frame,  one  can  re¬ 
duce  the  equation  of  the  hypersurface  to  the  form  t  =  0  . 
Lifschitz  and  Khalatnikov  found  that,  near  the  singularity, 
the  geometry  of  space-time  and  the  fluid  behave  in  the 
following  way: 


o  r  2Pi  c.*,9 

ds2  =  [x  l.,-  !.•  +  t 


2P?  2P  o  n  ,  i  ,  i  ,  2 

^  t  m.m.  +  x  3  n.n.Jdx  dxu  -  dx  , 
i  J  i  j  i  J  3 


(0)  -2(l-p3)  (0) 

U  =  U  X  ,  U.  =  U.  't 

J  J 


( 1-p  3 ) /2  (0) 

^  ,  U  =  U  T 


( 40A) 

( 3p3“l)/2 


4  4 


( 40B) 


*  Latin  indices  run  from  1  to  3;  greek,  1  to  4.  x4  =  x 
is  the  proper  time  of  the  matter. 


' 

- 
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Here  y  is  the  mass-energy  density  and  the  four- 

velocity  of  the  fluid.  The  quantities  p^,  1^,  mu  ,  n^ ,  y^), 

and  ua(0)  are  17  functions  of  the  space  coordinates 

i  * 

x  ,  which  are  connected  by  7  algebraic  relations  .  Lifschitz. 

and  Khalatnikov  further  asserted  that  "It  can  be  shown 

that  the  higher  terms  of  the  expansion  of  the  metric  con¬ 
tain  no  other  arbitrary  functions."  If  this  is  correct, 
then  the  Li fschit z-Khalatnikov  solution  contains  altogether 
ten  (17-7)  arbitrary  functions  of  the  space  coordinates, 
three  of  which  can  be  fixed  at  will  by  performing  a  co¬ 
ordinate  transformation  among  the  x1 ,  leaving  us  with  7 
physically  arbitrary  functions,  which  is  one  short  of  the 
number  required  of  the  general  solution. 

Thus  the  results  of  the  investigation  of 

*  * 

Lifschitz  and  Khalatnikov  show  that  space-like  singulari¬ 
ties  characterized  by  an  infinite  density  (such  as  that 
which  occurs  in  the  spherical  collapse  of  dust-c f. Sect . 2 ) 
cannot  be  created  by  realistic  gravitational  collapse;  in 
particular,  a  slight  perturbation  of  the  sphericity  will 
introduce  qualitative  changes. 


But  we  note  here  that  whether  the  employment  of  synchro- 
nons  frames  might  entail  some  loss  of  generality  is  a 
question  that  has  been  raised  3  *  . 

Two  of  these  relations  read:  Ep^  =  £p^2  =  1.  Also, 
by  convention,  -1/3  <  pi  <  0  <  p2  <  2/3  <  P3  <  1  . 
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(B )  Penrose  theorem. 

The  temptation  (which  one  might  feel)  to  genera¬ 
lize  the  results  of  Lifschitz  and  Khalatnikov  was  some¬ 
what  checked  by  the  more  recent  piece  of  work  of  Penrose*'^' 

and  its  further  ramifications,  mainly  in  the  hand  of 

(  47 )  / 

Hawking  .  Penrose  theorem  states  that  space-time  (con¬ 
taining  a  collapsing  system)  cannot  be  null-complete  if 

(1)  At  some  initial  instant  of  time,  the  universe 
has  an  infinite  volume; 

(2)  There  exists  a  Cauchy  hypersurface  H.  (A  Cauchy 
hypersurface  is  defined  as  a  complete  connected  space¬ 
like  3-surface  which  intersects  every  time-like  or  null 
line  once  and  only  once.  This  means  that  data  given  on  the 
hypersurface  are  sufficient  to  fix  the  entire  future  evolu¬ 
tion  of  the  universe).  Condition  (1)  implies  that  H  is 
open ; 

(3)  In  the  manifold  M+  which  represents  the  future 
time  development  of  H,  a  time  direction  can  be  globally 
assigned.  (Causality  condition); 

(4)  At  every  point  of  M+  ,  the  local  energy 

measured  by  any  observer  is  non-negative.  That  is,  if 

T  is  the  energy-moment  urn  tensor,  then  T  tMtV  >  0  for 

y  v  yv 

u 

all  timelike  vectors  tK  ; 


■ 
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and 


(5)  There  exists  in  M+  a  "trapped  surface"  T  - 

defined  generally  as  a  closed,  spacelike  2-surface  with  the 
property  that  the  null  geodesics  of  each  of  the  two  systems 
(of  outgoing  and  incoming  null  geodesics)  which  meet  T 
orthogonally  converge  locally  in  future  directions  at  T. 


For  spherically  symmetric  collapse,  trapped 

surfaces  exist  everywhere  inside  the  event  horizon  (cf. 

* 

fig.  4).  One  could  plausibly  argue  that  ,  provided  devia¬ 
tion  from  sphericity  is  not  too  large,  trapped  surfaces 
will  still  exist.  Hence,  if  conditions  (1)  to  (4)  are 
accepted,  it  appears  that  singularities  cannot  be  avoided 
by  destroying  the  spherical  symmetry. 


In  case  space-time  does  not  admit  a  Cauchy  sur¬ 
face  (as  happens  for  example  with  the  geometries  of  Kerr 
and  of  Reissner-Nordstrom)  ,  Hawking  and  Ellis  have 

given  a  modified  version  of  the  Penrose  theorem. 


The  rigorous  proof 0f  the  Penrose  theorem 

requires  advanced  techniques  from  differential  geometry  and 

topology,  which  lie  beyond  the  scope  of  the  present  work. 

(14) 

But  essentially,  the  proof  runs  as  follows  : 


# 


But  see  (C),  below. 


' 


L 

.  .  ' 
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Assume  that  all  the  hypotheses  of  the  theorem 
are  satisfied,  and  also  that  space-time  is_  null-geodesic 
complete.  Consider  the  null  geodesics  which  issue  perpen¬ 
dicularly  from  the  trapped  surface  T.  Since  T  is  trap¬ 
ped,  adjacent  null  geodesics  will  start  to  converge,  and 
by  extending  them  to  indefinite  values  of  their  affine 
parameters  (null-geodesic  completeness),  they  will  continue 
to  converge  (energy  condition  (4))  until  they  cross.  In 
this  way,  the  totality  of  such  null  geodesic  segments  forms 
a  closed  3-dimensional  null  hypersurface,  which  can  be 
approximated  arbitrarily  closely  by  a  closed  space-like 
3-surface  B.  Construct  all  time-like  geodesics  orthogonal 
to  B,  and  extend  them  Into  the  past  until  they  reach  the 
Cauchy  hypersurface  H  (condition  (2)).  This  family  of 
timelike  geodesics  provides  a  smooth,  nearly  one-to-one 
mapping  of  the  closed  hypersurface  B  onto  the  open  hyper¬ 
surface  H.  But  such  a  mapping  is  impossible  according  to 
a  fundamental  theorem  of  topology.  Hence  a  contradiction 
is  produced,  and  the  theorem  is  proved. 

The  Penrose-type  theorems  are  not  incompatible 
with  the  implication  of  the  works  of  LIfschitz  and 
Khalatnikov  as  we  stated  it,  for  these  theorems  say  nothing 
about  the  nature  of  the  singularities  that  may  occur.  The 
Schwar zschild  r  =  0  singularity,  which  is  hit  by  every 
time-like  and  null  geodesic  that  passes  inside  the 
Schwarzschild  surface,  is  actually  not  typical.  For  the 


. 
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Reissner-Nordstrom  solution,  for  instance,  only  radial  null 

or  space-like  geodesics  can  reach  the  singularity. 

(49) 

Grischuk  has  found  a  general  solution  to  the  gravita¬ 

tional  field  equations  for  dust  which  contains  a  time-like 
singularity.  All  these,  together  with  the  results  of 
Lifschitz  and  Khalatnikov,  indicate^^  that  the  singulari¬ 
ties,  if  they  do  occur  in  asymmetric  gravitational  collapse, 
should  consist  of  isolated  points  or  timelike  surfaces  which 
most  of  the  world  lines  managed  to  avoid. 

( C )  Israel’s  Theorem. 

The  significance  of  the  Penrose-type  theorems 
hinges  on  the  assumption  that  when  the  strict  sphericity 
is  only  very  slightly  perturbed,  the  Schwarzschild  surface 

will  -  though  distorted  -  still  retain  its  essential 

characteristics  as  an  event  horizon,  since  otherwise  the 
occurence  of  trapped  surfaces  for  the  general  case  cannot 
be  plausibly  argued  for.  For  our  purpose,  an  event  horizon 
(with  respect  to  observers  far  from  the  collapsing  system) 
may  be  defined  as  a  hypersurface  which  nowhere  extends  to 
spatial  infinity,  and  which  divides  space-time  events  into 
two  non-empty  sets,  one  consisting  of  events  that  in 
principle  could  be  observed  by  some  of  the  observers,  the 
other  consisting  of  events  that  forever  lie  outside  their 
power  of  observation.  In  the  special  case  of  static 
space-times  (i.e.,  those  independent  of  time  and  time- 


'  v  '  ' 


’  . 

I 
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' 
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reversal)  which  are  not  necessarily  spherically  symmetric, 
the  surfaces  g  =0  are  the  natural  analogue  of  the 

4  4 

Schwarzschild  surface. ^0) 

The  production  of  an  event  horizon  is  itself  a 

remarkable  feature  of  spherical  collapse,  and  perhaps  the 

assumption  of  its  continual  existence  for  more  general 

( 51 ) 

cases  should  not  be  lightly  made.  Indeed,  Israel  has 

very  recently  shown  that  any  static  perturbation  that  des¬ 
troys  the  spherical  symmetry  of  a  source  of  the  Schwarzschild 

field,  also  destroys  the  event  horizon.  More  precisely, 

( 51 ) 

Israel’s  theorem  states  that 

Among  all  static,  asymptotically  flat  vacuum 
fields  with  closed,  simply-connected  equipotential 
surfaces  g  =  const  ,  Schwarzschild ’ s  solution 

4  4 

is  the  only  one  which  has  a  non-singular  event 
horizon  g  =  0  . 

4  4 

What  is  the  bearing  of  Israel's  theorem  on  the 

(  52 ) 

question  of  asymmetric  collapse  ?  Suppose  that  an 
event  horizon  forms  during  the  collapse  of  an  asymmetric 
(non-rotating)  star.  Because  of  the  slowing  down  of  pro¬ 
cesses  for  the  external  observer,  it  seems  permissible  to 
regard  the  limiting  external  field  as  static.  Thus,  if 
the  regular  event  horizon  is  to  be  retained,  the  star  must, 
before  passing  through  it,  divest  itself  of  all  quadrupole 
and  higher  moments;  then  the  Penrose-type  theorems  may 


'  i 
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become  operable,  and  the  collapse  should  proceed  in 

qualitatively  much  the  same  way  as  previously  described, 

* 

ending  in  a  state  of  singularity  .  The  chief  theoretical 

difficulty  is  to  find  a  plausible  mechanism  to  accomplish 
( 52 ) 

this  task  .  But  unless  this  is  accomplished,  Israel's 

theorem  either  prohibits  (due  to  the  formation  of  singu- 

( 59 ) 

larities  on  the  event  horizon  )  or  renders  unlikely 
(due  to  the  absence  of  an  event  horizon)  the  development  of 
trapped  surfaces.  There  would  then  be  no  reasons  why  the 

4* 

star  should  not  eventually  bounce  out  or,  indeed,  pulsate  , 
as  happens  generally  in  the  corresponding  Newtonian  case. 


* 

Or  the  star  may  re-emerge  into  another  universe.  Cf. 
Chapter  IV. 

t  For  specific  considerations  in  favor  of  pulsating  stars 
as  a  model  for  the  baffling  quasi-stellar  sources,  see 
Israel  26  . 


. 
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CHAPTER  IV 

GRAVITATIONAL  BOUNCE* 


1 .  Introduction . 

We  shall  be  concerned  in  this  chapter  with  the 
relativistic  gravitational  collapse  of  a  charged  spherical 
shell  falling  in  a  spherically  symmetric  external  field. 

The  study  of  such  simple  artificial  problems,  while  of  no 
direct  relevance  to  astrophysics,  can  nonetheless  serve  a 
useful  purpose,  since  it  brings  into  relief  basic  issues 
of  principle  in  the  general  relativistic  theory  of  collapse 
which  are  still  far  from  understood. 

As  we  saw  in  the  last  chapter,  general  relativity 
leads  to  the  following  picture  for  the  evolution  of  a  con¬ 
tracting  spherical  body.  Once  the  compression  passes  a 
certain  critical  limit,  characterized  roughly  by  the 
Newtonian  potential  becoming  comparable  with  c2,  the  sub¬ 
sequent  history  is  one  of  continuing  collapse  which  cannot 
be  halted  by  pressure  forces.  The  irreversibility  of  this 
picture  is  surprising,  and  differs  radically  from  the 
corresponding  Newtonian  picture,  where  the  motion  is  in 
general  oscillatory.  If  one  examines  the  relativistic 

*  ( 55 ) 

Based  on  V.  de  la  Cruz  and  W.  Israel  . 
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derivation  to  see  how  the  element  of  irreversibility  enters, 
one  finds  that  it  stems  from  two  largely  unconnected  causes. 

(A)  External  irreversibility:  development  o f  an 
event  horizon  at  r  =  2m .  The  surface  of  the  contracting 
spherical  body  passes  (in  finite  proper  time)  within  the 
critical  Schwarzschild  sphere  r  =  2m.  To  an  external  ob¬ 
server,  light  emitted  from  this  sphere  suffers  infinite 
gravitational  and  Doppler  red-shift,  and  r  =  2m  therefore 
appears  as  an  event  horizon  which  the  contracting  body  seems 
to  be  approaching  asymptotically  as  t  ->  00  .  If  ordinary 
ideas  of  causality  are  to  be  maintained,  he  can  never  see 
the  body  re-emerge  from  this  sphere. 

(B)  Intrinsic  irreversibility:  spacelike  character  of 
the  curves  r  =  const  near  r  =  0  .  The  exterior 
(Schwarzschild)  field  of  the  body,  analytically  extended 

to  r  =  0  ,  has  the  property  that  the  curves  r  =  const  <  2m 
are  space-like.  The  history  of  a  particle  on  the  surface 
of  the  body  is  a  time-like  curve  of  the  exterior  manifold. 

It  is  easy  to  see  that  the  particle  can  reverse  its  inward 
motion  at  r  =  rQ  <  2m  only  if  (a)  its  world-line  is  mo¬ 
mentarily  space-'l'ike ,  and  (b)  it  subsequently  travels  into 
the  past.  Assuming  that  classical  general  relativity  re¬ 
mains  valid  even  under  the  extreme  conditions  prevailing 
near  r  =  0  ,  one  is  forced  to  the  conclusion  that  no  re¬ 
bound  is  possible  and  that  the  entire  mass  piles  up  irre¬ 
versibly  on  the  singular  curve  r  =  0. 


' 


rt 
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Of  these  two  arguments,  (A)  seems  (at  least  in 
so  far  as  the  condition  of  spherical  symmetry  is  strictly 
met  ;  cf.  Sect. 4,  Chap. Ill)  on  surer  ground,  since  it  does 
not  depend  on  extrapolation  to  extreme  conditions.  For 
masses  of  astrophysical  interest,  compression  to  r  =  2m 
does  not  produce  immoderate  densities  or  curvatures.  Modi¬ 
fication  due  to  quantized  gravitation,  possible  inapplica¬ 
bility  of  Einstein's  field  equations  in  regions  of  extremely 
intense  gravitational  fields,  or  other  new  physical  effects 
of  an  unanticipated  kind  might  profoundly  affect  the  situa¬ 
tion  near  r  =  0,  but  should  not  be  important  near  r  =  2m. 

The  exact  nature  of  such  modifications  is,  of  course,  unknown. 
One  could  try  to  take  their  effects  into  account  in  a  crude 
way  by  supposing  that  the  standard  Schwarzschild  metric  is 
modified  to 

ds2  =  ( l-2m/r+a/r2+  •  •  •  )“  \lr  2-(l-2m/r+b/r2  +  .  .  .  )dt 2  +  r2dft2  , 

dfi2  =  d02  +  sin2  0  d(f)2  . 

If  the  usual  astronomical  predictions  of  Einstein's  theory 
are  to  be  preserved,  the  constants  a  and  b  would  have 
to  be  small  compared  with  m2 .  In  that  case,  we  still  have 
an  event  horizon  (the  sphere  on  which  g  =0),  but  the 

4  4 

arguments  of  (B)  are  clearly  liable  to  break  down.  The 
possibility  cannot  be  ruled  out  the  a  collapsing  spheri¬ 
cal  body  reverses  its  motion  near  r  =  0  and  re-expands . 

The  intriguing  question  is  how  such  a  picture  can  be  recon- 


, 
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died,  with  the  apparent  irreversibility  of  the  collapse  as 
seen  by  an  external  observer. 


We  shall  illustrate  some  of  the  possibilities  by 

•4* 

focussing  attention  on  the  special  line-element 


ds2  =  f_i  dr2  -  fdt 2  +  r2dft2  , 
f(r)  =  1  -  2m/r  +  e2/r2  , 


(1) 


( 2|  1 ) 

which  is  formally  the  Reissner-Nordst rom  metric  for  the 

gravitational  field  of  a  charged  particle.  In  Sects.  2 
and  3  we  derive  the  equation  of  motion  of  a  thin  spherical 
shell  in  such  a  field.  It  is  found  (Sects. 4,7)  that  bounce 
can  occur  under  a  great  variety  of  conditions.  In  particu¬ 
lar,  the  shell  can  bounce  inside  an  event  horizon.  In  that 
case  the  manifold  represented  by  (1)  is  incomplete.  If  it 

is  extended  analytically  (Sect. 5)  in  the  manner  of  Graves 

( 4l ) 

and  Brill  ,  the  extended  space-time  appears  as  a  perio¬ 
dic  lattice  of  geometrically  similar  asymptotically  flat 
spaces,  joined  by  "tunnels"  in  which  space  is  closed.  Re- 
emergence  of  the  rebounding  shell  from  the  event  horizon 
then  takes  place  in  a  new,  distinct  space  (Sects. 6, 7). 


Our  results  are  similar  to  those  recently  obtained 
by  Novikov^*^  ,  who  has  considered  the  homologous  collapse 
of  a  uniformly  charged  ball  of  dust.  The  shell  model  has 


An  an  Illustration,  it  may  be  remarked  that  (1)  represents 
the  external  field  of  an  (uncharged)  spherical  body  in 
the  theory  of  F.  Hoyle  and  J.V.  Narlikar  3  . 
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the  advantage  that  complete  solution  can  be  given  in  a  sim¬ 
ple  explicit  form.  Novikov's  results  can  be  recovered  as 
a  special  case  (Sect. 7). 

Sect.  8  concludes  the  chapter  with  some  general 

remarks . 

2 .  Dynamics  of  a  Thin  Shell. 

( S  7 ) 

We  shall  adopt  the  approach  to  thin  shells 

in  general  relativity  in  which  their  histories  are  charac¬ 
terized  in  a  purely  geometrical  way  by  the  extrinsic 
curvatures  of  their  imbeddings  in  space-time. 

Let  V  be  a  four-dimensional  Riemannian  manifold 

of  metric  (+++-)  ,  with  coordinates  x  and  metric  tensor 

g^g  .  Let  Z  be  a  time-like  hypersurface  in  V  ,  i.e., 

a 

the  unit  normal  n  to  Z  is  space-like: 

n  na  =  +1  .  (2) 

a  v  ' 

In  terms  of  the  coordinates  xa  in  V  ,  the  parametric 
equations  of  Z  are  of  the  form 

xa  =  fa  USC2^3)  ,  (3) 

*  i 

where  6  a re  intrinsic  coordinates  of  Z  .  An  infi- 


Throughout  this  chapter,  Greek  indices  refer  to  4-dimen¬ 
sional,  Latin  indices  to  3-dimensional  quantities. 
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a 


nitesimal  displacement  dx  in  £  is  given  by 


dxa  =  Ofa/8e1)d^1  , 


(4) 


which  defines  the  triad  of  holonomic  basis  vectors  tangent 
to  £  and  associated  with  : 


e^)  =  3 f  /3€ 


(5) 


The  intrinsic  geometrical  properties  of  £  are  completely 


determined  by  its  metric  3-tensor  g. 


a 


e / . s  e / o \  ,  and 


Jij  (i )  ( j  )ot 

intrinsic  operations  such  as  covariant  differentiation  may 


a 


be  analogously  defined  in  £  .  Thus,  if  A  is  a  vector 
field  tangent  to  £  ,  we  can  associate  with  it  a  3-vector 
field  A^  in  £  given  by 


n  «  oi  « a  .  i  a 
A.  =  A  e /•  .  v  ,A  =  A  e  /  .  , 

i  a  (l)  5  (l)  5 


(6) 


and  define  the  covariant  derivative  of  A^  with  respect  to 
5^  to  be 


A.  .  =  3A./35J'  -  Ak  I\  .  . 

i;j  1 


(7) 


where 


rksij  “  *  (S  Slk/3CJ  +  3  g.k/3f-3  glj/3Ck) 


are  the  3-dimensional  analogue  of  the  4— dimensional 

Christoffel  symbols  T  0  . 

°  a ,  3y 


Extrinsic  properties  of  £  associated  with  its 
imbedding  in  V  are  measured  by  the  absolute  derivatives 
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ct  i 

<5n  /6£  of  the  unit  normal.  We  recall  that  the  absolute 
derivative  of  a  vector  function  Ba  defined  on  a  curve 
parametrized  by  t  is 


6Ba/St  =  8Ba/9t  +  BAf(f  dxVdt 

Ay 

Since  ,  from  ( 2 ) , 


(8) 


r.  a  ,  ~  _i 

i .  e  .  ,  on  /6  6  are 


na6na/6^1  =  0 

perpendicular  to 

r  a  .  c,  >.  i  „  j  a 

6  n  / 6 E,  =  K .  e /  „  > 
i  ( J  > 


(9) 

a 

n  ,  we  can  write 

(10) 


thus  defining  the  extrinsic  curvature  3-tensor 
Z  .  Explicitly, 


of 


K.  .  =  -  n  Se  “,/<5£i 

ij  a  (jr 

=  -  n  {d2fa/d^d^  + 
a 

+  OfVa^1)  OfVn1)  > 

Ay 


(ID 


If 


is  a  vector  field  tangent  to 


£ 


then1"  ^  ^ 


6Aa/6?1 


n  1  jr  a 

A  K.  .  n 
ij 


(12) 


( 57 ) 

In  general,  one  has  the  following  four  relations  between 

the  extrinsic  curvature  3-tensor  K. .  and  the  normal  com- 
ponents  of  the  Einstein  tensor  on  E  : 


.ab 


kab 


Kb  ,  -  9  K  = 

a  ;b  a 
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Tr  2  o  o  Ot  $ 

K  =  -  2  G  0n  n  , 
a8 

(13) 

^  a  8 

-  G  De ,  \  n 

a8  (a) 

(14) 

ij 


Ko  „  3  and  R  is  the  intrinsic  curvature 
i  J  5 


Here ,  K  =  g' 

invariant  of  Z  .  Each  term  of  (13)  and  (14)  is  a  4-scalar, 


i.e.  independent  of  the  coordinates  x 


a 


Now  let  the  time-like  hypersurface  Z  divide 

space-time  into  two  parts  V_  ,  V+  which  both  contain  Z 

as  part  of  their  boundaries  and  are  otherwise  disjoint. 

Let  n  (directed  from  V_  to  V  )  be  the  unit  space-like 

normal  to  Z  ,  and  let  K  ,  -  ,  K  ,  +  denote  the  extrinsic 

5  ab  *  ab 

curvatures  of  Z  associated  with  its  imbeddings  in 

( 57 ) 

V_ ,  V+  respectively.  Then  Z  is  the  history  of  a  thin 

shell  if 


^ab 


(15) 


is  nonvanishing.  The  surface  energy  tensor 

(57) 

shell  is  given  by  the  Lanczos  equations 


S  ,  of  the 
ab 


Y 


ab 


3ab 


Y  = 


=  -  8 


7T  S 


ab 


(Y  = 


ab 


yab> 


(16) 


the  analogue  of  Einstein’s  field  equations 


G 


a3 


8  7T 


(17) 


for  the  surrounding  continuous  medium. 
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eight  relations 


Corresponding  to  (13),  (14)  are  the  following 
t 


3R  -  K  ,  Kab  +  K2 I  1  ~ 


Lab 


b  i  + 

K  -  9  K  ~ 

a  ;b  a  1 


=  -  2  G  pnanB  | + 
a3  1 


a  i  ± 
a3  c(a)' 


-  G  „  e, _ ,n 


(18) 

(19) 


The  jump  of  (18)  across  Z  is 


2  [0  nV] 


( -K  Kab+K  2 ) 
ab 


-  ( -K  Kab+K2 ) 
ab 


(Kab+  -  Kab  )(Kab  +  Kab)  +  (K+-K  )(K++K  ) 


•Yab  (Ka+  +  Kab)  +  y(K++K  ) 


by  (15) 


/  \  /  „ab  .  ,,ab  \ 

-  (Yab-YSab)(K  +  +  K  -> 


8rr  S  (Kab  +  Kab) 
ab  + 


by  (16)  , 


or,  in  view  of  (17), 


2  [T  Qn°n3]  =  Sab  (K  *  +  K  ,  ). 

a  3  ab  ab 


(20) 


Similarly,  (19)  yields 


FT  e  ?  nB  1  = 
La3  (a)  n  J 


-  S 


a  ;b 


(21) 


t  Limits  of  a  field  quantity  1  as  an  event  P  on  Z  is 
approached  from  V_ ,  V+  respectively  are  denoted  by 
¥ |  ,  ¥ |  .  In  Sects.  2  and  3  square_brackets  denote 

jump  discontinuities:  [¥]  =  ¥ j  -  ¥  | 


69 


We  now  specialize  to  the  case  of  a  coherent 
shell  of  dust,  characterized  by  the  surface  energy  tensor 

o 3-d  a  b  , 

S  =  a  u  u  ,  ( 22 ) 

ga  q 

where  the  unit  time-like  vector  u  =  d£  /dx  tangent  to 

Z  represents  (t  being  the  proper  time)  the  4-velocity  of 

the  dust  particles,  and  a  is  the  sum  of  their  rest  masses 

per  unit  area.  The  dynamics  of  the  shell  in  vacuo  was 

(  57 ) 

considered  in  .  The  case  where  the  shell  falls  in  a 
continuous  medium  with  nonvanishing  energy  tensor  is  a 
straight  forward  generalization. 


From  (22)  , 


0  b  b  ,  ,  b  v 

S  .  =  u  ,  o  u  +(au)  .  u  , 

a  -b  a;b  •  b  p  a 


a  0  b  a  b  ,  ,  b ,  a 

u  sa  ;b  =  ua;bu  0  u  +  (o  u  ■>  ;buau 


/  b 

=  -  (0  u  )  b  , 


and  by  use  of  (21),  we  obtain 


/  b\,  ar  m  Ot  3  1  I'm  ■,•.0'  3  ~l 

(a  u  );b  =  u  [Tag  e(a)  n  ]  =  [T^  u  n  ] 


(23) 


and 


b 

a  u  ,  u  = 
a  ;b 


-  (6°  +  u  ub)  [T  e?  v  n13]  , 
a  a  a  3  ( c )  5 


(24) 


where 


u 


a  i  ± 


a  a  |  ± 
u  e ,  \ 

(a)  1 


rv  i  + 

=  dx  /dx 


. 
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represent  the  4-velocity  of  the  dust  particles  as  measured 
in  V+  ,  V_  .  The  4-acceleration  6u  / 6 t  | _  can  be  written, 
according  to  (12),  as 


6ua/6i|±  =  (6ua/6^1)  d^/dr 


a  b  ot  T.  a  b  oi, 

=  u  . ,  u  e  /  >.  -  K  ,  u  u  n 
p  b  ( a )  ab 


Prom  this  and  (20)  we  immediately  obtain 


O  0i  y  O  (  —  ,  £  Ot  /  rs  I  “  Opm  Ot  3  1 

an  6u  /6t  +  an  ou  /6t  =  -  2[T  nn  n 

a  1  a  1  a3 


(25) 


and 

n  6ua/6x|+  -  n  6ua/5i|  =  -  y  ,uaulD  =  4Tra  (26) 

a  1  a  1  1 ab 

with  the  aid  of  (16)  and  (22). 


3 .  Charged  Spherical  Shell  in  a  Spheri-symmetri c 
Electrovac  Field. 

We  consider  a  charged  spherical  shell  of  dust 
falling  in  the  electrovac  field  produced  by  a  spherically 
symmetric  concentration  of  mass  and  charge  at  its  center. 
For  such  a  spherically  symmetric  (static  or  nonstatic) 

/  r  O  \ 

universe,  an  extension  ^  ;  of  Birkhoff’s  theorem  shows  that 
the  line  element  is  reducible  to  the  standard  Reissner- 

Nordstrom  metric  (1)  -  with  appropriate  parameters  e,m 

-  in  any  region  free  of  matter. 

Let  r  =  R(t)  be  the  equation  of  £  ,  the  history 
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of  the  shell,  and 


(ds2)^  =  {R(x)}2  dft^  -  dx2 


(27) 


be  its  intrinsic  metric,  so  that  t  is  the  proper  time 
measured  along  the  streamlines  0,(|)  =  const.  The  interior 
and  exterior  line  elements  may  be  written 


(ds2)_  =  {f  (r)}-  dr2  +  r2d^2  -  f  (r)dt  2,  (r<R(x)), 

(28) 

( ds  2 )  =  {f  (r ) }_1  dr2  +  r2dfi2  -  fx(r)dt  2  ( r >R ( t ) ) , 

+  +  +  +  (29) 


where 


f  (r)  =  1  -  2m  /r  +  e  2/r2  , 

i  i 

(30) 

f, (r)  =  1  -  2m  /r  +  e  2/r^  . 

+  2  2 


Thus,  the  shell  has  charge  e  -e  and  gravitational  mass 

2  1 

m  -m 

2  i 


Both  (28)  and  (29)  must  induce  the  same  intrinsic 


metric,  namely  (27),  on  Z  ,  Comparison  of  the  coefficients 
of  d^2  confirms  that  the  interior  and  exterior  radial 


coordinates  agree  on  Z  .  Further, 


dx2  =  f  ( R ) dt  2  -  {f  (R) }-  dR2  =  fx(R)dt  2  -  {f  (R)}“dR2  . 

—  —  —  Hr  Hr  T 

(3D 


This  fixes  the  relation  between  t_  and  t+  on  Z  ,  and 
verifies,  as  expected,  that  the  simultaneous  imbedding  of 
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£  in  V  ,  V  is  possible. 

—  1 

We  proceed  to  write  out  explicitly  the  dynamical 
equations  (23)  to  (26).  Since  u  ,  n  are  orthogonal  unit 
vectors  in  the  2-space  of  r(=  x|)  ,  t + ( =  x*)  a  we  have 

u^  =  dx^/dr  =  (R,o,OjX_,  )  5  R  e  dR/’di  ,  (32) 

n*  =  (X+,o,o,-R)  ,  (33) 

where 


X  =  dt./dt 

.  + 

=  (f+(R)  +  R2>Vf+  by  (31),  (34) 

ot  — 

with  corresponding  expressions  for  u  ,  n&  „  By  intrinsic 
differentiation  of  u+u^  =  -1  and  use  of  (32)  we  find 

0  =  u  6ua/6 t | +  =  f"1  R  6u1/5t  -  f  X  6u'/6t  , 

and  this  may  be  used  to  eliminate  6u^/6t  from 

n  6ua/5 1  | +  =  X  6u1/6t  -  R  6uV6t|  + 
a 

yielding 

n  6ua/6ij+  =  (fX)-1  6u‘/5t 
a 

°*  .  j  + 

=  (fX)"1  {R  +  r}  uAuy} 

Ay 

=  (f+X+)_1  {R  +  %  df+(R)/dR}  ,  (35) 

oc  i _ 

with  a  corresponding  expression  for  n&  6u  / 6 t | 
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The  Reissner-Nordstrom  metric 
with  the  energy  tensor 


(1)  is  associated 


(59) 


-  T4 

4 

_  _  iji  1  _  iji  2 

1  2 

— 

T3  =  e2/87rr‘f 

3 

( oth 

er  components 

zero) ,  so  that 

in  our  case 

m3  a  ± 

T  u  nn 
a  3 

= 

0  , 

[T6 

a 

nanR]  =  [T1 
p  i 

n1 

n  +  T4  n4n  ] 

1  4  4 

=  [T1 

l 

n 

a 

a 

n  J 

=  (e 

l 

2 

e  2)/87rr4  . 

2 

From 

(23)  and  (37) 

9 

<  b  \ 
(au  ) ;b 

- 

0 

(36) 


(37) 


(38) 


(39) 


expressing  the  conservation  of  the  proper  mass  of  the 
shell.  For  the  co-moving  coordinates  employed  in  (27), 
u13  =  (0,0,1)  and  (39)  symplifies  to 


d  ln(R2a)/dx  =  0  , 


or 


47rR2a  =  k  =  const  .  (40) 

Equations  (24)  are  here  identically  satisfied. 
Equations  (25)  and  (26)  yield, when  (35),  (38),  and  (40) 
are  substituted  into  them. 


.  •  •  • 

1 

9 

. 

' 


74 


(e  2-e  2 ) /kR2  =  (f  X,)  1  {R  +  m  /R2  -  e  2/R3}  + 

2  1  +  +  2  2 

+  (f  X  )-J  {R  +  m  /R2  -  e  2/R3},  Mn  . 

-  -  i  i  ( 41 ) 

k/R2  =  (f. X,)-1  {R  +  m  /R2  -  e  2/R3}  - 

+  f  2  2 

-  (f  X  )_1  {R  +  m  /R2  -  e  2/R3}  . 

_  _  1  i 

(42) 


We  recall  from  (34)  that 
f,X,  =  {1  +  R2  -  2m  /R  +  e  2/R2}^  , 

-t-  -r  2  2 

f  X  =  {1  +  R2  -  2m  /R  +  e  2/R2}^  . 

-  i  l 

Setting  y  |  h(l  +  R2)  so  that  R  =  dy/dR,  and 
employing  x  =  1/R  as  independent  variable,  we  obtain  from 
(4l),  (42)  the  equivalent  set  of  equations 

(e  2-e  2)/k  +  k  =  -  2  ( 2y  -  2m  x  +  e  2x2)^  , 

2  1  dX  2  2 

(43) 

(e  2-e  2)/k  -  k  =  -  2  ( 2y  -  2m  x  +  e  2x2)^  . 

21  dx  ^  i  i 

(44) 


Either  of  these  may  be  integrated  at  once  and  the 

constant  of  integration  determined  from  the  other.  In  this 

* 

way  we  arrive  at  the  first  integral 


This  may  be  compared  with  the  corresponding  equation  of 
motion  in  Newtonian  theory : 

^(m  -m  )(dR/dt)2  +  ^[e  2+e  2-(m2+m2 ) ]/R  =  const  . 

2  1  12  12 


. 


.  1  'J  •  '  f/O  - 

■  !  • 
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1  +  (dR/dr ) 2  =  A  +  B/R  +  C/R2  ,  (45) 


where 


A  =  (m  -m  )2/k2  (46) 

2  1 

is  introduced  for  convenience,  and 


B  =  m  +  m  -  A(e  2-e  2)/(m  -m  )  ,  (47) 

12  2  1  2  1 

4C  =  A(e  2-e  2)2/(m  -m  )2  -  2(e  2+e  2)  +  A_1(m  -m  )2  . 

2  1  2  1  12  2  1 

(48) 


From  (40)  and  (46)  , 

4ttR2c  =  A  i(m  -m  )  , 

2  1 

which  enables  us  to  interpret  the  constant  A  in  terms 
of  the  binding  energy  W  ,  since 

-  W  =  (m  -m  )(1-A-^) 

2  1 


represents  the  difference  between  the  gravitational  mass 
m  -m  of  the  shell  (i.e.,  its  total  energy)  and  the  sum 

2  i 

of  the  rest  masses  of  its  constituent  particles.  It  thus 
represents  the  contribution  to  the  shell's  gravitational 
mass  due  to  its  kinetic  and  potential  energies. 


The  further  integration  of  (45)  would  be  elemen¬ 
tary.  However,  the  various  physical  possibilities  emerge 
more  clearly  from  a  qualitative  description  of  a  few 
representative  special  cases.  This  will  be  our  aim  in  the 


■.  ;  H\a  +  a  =  • 


b  f :  B 


'A  -I-  (s  e  +  :;  9)S  -  ■  (  m-  m)\ 

bn  6 

.. 

'  • 

■ 

-id 

■  ■ 
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next  few  sections. 

4 .  Charged  Shell  in  Vacuo. 

Suppose  that  no  mass  or  charge  is  present  apart 

from  the  shell  itself  (mass  m,  charge  e).  Then 

e  =m  =  0  ,  e  =  e  ,  m  =  m  and  the  equation  of  motion 
11  2  2 

(45)  reduces  to 

(1  +  (dR/dx)2}^  =  a  -  b/R  ,  (49) 

where  we  have  written  A  =  a  and 

b  =  ( a2e2-m2 )/2am.  (50) 

If  b  <  0  ,  we  can  imagine  the  shell  as  starting, 
either  from  infinity  with  initial  velocity  R  =  -  (a2-l)^ 
(for  a  >  1),  or  from  rest  at  a  finite  maximal  radius 
R  =  |b|/(l-a)  (for  0  <  a  <  1).  It  accelerates  as 
it  falls  inward  and,  upon  reaching  R  =  0  ,  produces  a 
singularity.  The  subsequent  history  is  therefore  a  matter 
of  conjecture  (the  possibility  of  a  rebound  is,  of  course, 
not  excluded ) . 

More  definite  conclusions  can  be  drawn  when 
b  >  0  (always  obtainable  for  a  shell  with  given  mass 
and  nonvanishing  charge  by  taking  a  sufficiently  large). 

In  this  case,  it  is  necessary  that  a  >  1  .  The  shell  is 
impelled  inwards  from  infinity  with  initial  velocity 


•  .V  •  -  :  '  '  ’’  ■  -  f 

&  "  "  £  :  "  •  :  TfT 

.n  :sf:\ 

n  ■' 

•  . .  .  - 
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/  2  \  ^2 

(a  -1)  .  It  is  decelerated,  and  comes  to  rest  at  a  finite 
radius  Rmln  =  b/(a-l)  ,  then  re-expands  symmetrically  to 
infinity.  As  measured  by  a  co-moving  observer,  who  uses 
the  proper  time  t  ,  the  time  required  to  implode  from 
any  given  radius  Rq  >  Rm^n  and  re-expand  to  this  radius 
is 


At  =  2 ( a2-l) ” 1  /  [ ( a+1 )RQ-b ] [ ( a-1 )R  -b ]  + 

+  4ab (a2-l)~3/2  ln{/  ( a-1) [ ( a+1 )RQ-b ]/2b  + 

+  /  ~ ( a+1 ) [ ( a-1 )RQ-b ]/2b  } 

(51) 


and  is  therefore  finite. 

Since  x  is  related  to  the  time  t  of 
stationary  observers  in  the  interior  flat  domain  by, 
according  to  (31), 

dt_2  =  ( 1+R2 ) dx  2  , 
we  obtain  from  (49) 

( dR/ dt_ ) 2  =  [ ( a-1) R-b ] [ ( a+1 )R-b ]/( aR-b ) 2  .  (52) 

The  denominator  does  not  vanish  for  R  >  R  .  .  The 

-  mm 

motion  as  seen  by  an  interior  observer  is  thus  qualitative¬ 
ly  similar  to  the  intrinsic  description  just  given. 


To  an  external  observer,  however,  the  sequence 


(. 
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of  events  may  appear  quite  different.  Prom  (49)  and  (31) 
we  obtain 


(dt+)2  =  _i  +  R2 _ 

dR  f2  f[(a-l)R-b][(a+l)R-b] 

f(R)  =  [(R-m)2  +  e2-m2]/R2  , 


(53) 


as  the  equation  of  motion  in  terms  of  the  exterior  time 
coordinate  t+  (essentially  the  proper  time  of  a  station¬ 
ary  observer  with  large  radial  coordinate).  For  a  shell 
with  e  >  m  ,  f  never  vanishes  and  the  coordinate 
ra  t+  cover  the  complete  exterior  manifold:  qualitative¬ 
ly  the  motion  seen  externally  is  as  previously  described. 
However,  for  m2  >  e2  ,  f  vanishes  at 

R  =  m  +  (m2-e2 ) ^ 
i 

before  Rm^n  =  m(a+l)/2a  -  a(m2-e2 ) /2m( a-1 )  <  m  is 

reached.  The  time  At  needed  to  implode  from  any  given 

radius  R  >  R  to  R  is 
o  1  i 

At  ~  lim  ^  R  2(m2-e2)-^  ln(l/e)  for  m2  >  e2  , 

+  1 

e+0 

(54) 

~  lim  (m2/e)  for  m2  =  e2  . 


Thus  we  reach  the  curious  conclusion  that  an  external 
observer  never  sees  the  re-expanding  shell  if  e2  <  m2  ; 


re-expansion  to  arbitrarily  large  radius  nevertheless 
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occurs  in  finite  time  according  to  an  interior  or  a 
co-moving  observer. 

5 •  Analytic  Completion  of  Reissner-Nordstrom  Manifold 

for  e2  <  m2  . 

To  resolve  the  apparent  paradox  of  the  previous 
Section,  we  require  a  picture  of  the  exterior  manifold 
when  e2  <  m2  .  The  coordinates  r,  0,  (p ,  t  then  no 
longer  furnish  a  complete  map.  The  problem  of  analyti¬ 
cally  completing  the  Reissner-Nordstrom  manifold  has  been 
dealt  with  by  Graves  and  Brill^1^  (for  e2  <  m2 )  and 
by  Carter  (for  e2  =  m2)  .  We  shall  present  a  some¬ 

what  simplified  review. 

(i)  The  case  e2  =  m2  .  In  this  case  the  function 
f(r)  in  the  Reissner-Nordstrom  metric  (1)  becomes 

f (r )  =  ( 1-m/r ) 2  , 

and  the  coordinate  t  remains  time-like  for  all  r. 
Introduce  an  angular  time-like  coordinate  0  ,  with  range 
_  oo  <  0  <  oo  }  such  that  t/2m  =  tg  0  for  -  tt/2  <  0  <  tt/2. 
The  (formally)  extended  line  element 

ds2  =  (1-m/r)-2  dr2  +  r2  d!22  -  4m2  ( 1-m/r ) 2  ( d  tg0)2 


(55) 
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represents  a  periodic  space-time  which  has  a  geometrical 
singularity  at  r  =  0  and  is  otherwise  free  of  singulari¬ 
ties.  The  r,  0  map  is  subject  to  local  breakdown  on  the 
lines  r  =  m  ,  0  =  (n+^)7T  .  That  r  =  m  is  actually  a 
regular  part  of  the  manifold  can  be  verified  by  express¬ 
ing  the  line  element  in  a  form  which  is  manifestly  regular 
for  r  >  0  : 


ds 2  =  2  dv’dr  -  (l-m/r)2dv'2  +  r2d^2  (56) 

where  the  advanced  time  parameter  vf  is  analytically 
related  to  r,  t  by 

dv’  =  (1-m/r)-2  dr  +  dt  (r  >  m)  .  (57) 

In  the  v' j  r  chart  one  can,  for  instance,  follow  any 
incoming  radial  null  geodesic  v*  =  const  originating 
in  a  region  with  r  >  m  (e.g.  region  la  of  fig. 5) 
down  to  r  =  0  (in  region  Ilia) .  This  chart  thus  provides 
a  regular  mapping  of  two  adjoining  regions  such  as  la  and 
Ilia.  By  analogous  use  of  a  retarded  time  parameter  we 
can  construct  a  chart  for  Ilia  and  lb .  An  infinite  chain 
of  such  overlapping  coordinate  patches  enables  us  to 
follow  any  null  or  time-like  geodesic  down  to  the  singu¬ 
larity  at  r  =  0  or  to  indefinite  values  of  its  affine 
parameter.  Use  of  the  r,  0  map  means  that  allowance 
must  be  made  for  local  breakdowns,  but  has  the  advantage 
of  providing  a  clearer  over-all  picture. 
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(ii)  The  case  e2  <  m2  .  In  this  case,  the  quadra¬ 
tic  coefficient  f(r)  in  the  Reissner-Nordstrom  metric 
(1)  has  real  unequal  factors 

f(r)  =  (r-r  )(r-r  )/r2  ,  (0  <  r  <  r  )  ,  (58) 

12  2  1 

r  =  m  ±  (m2-e2)^  .  (59) 

1  2 

Writing  the  metric  in  the  form 

ds2  =  f  ( f “ 1 dr+dt ) ( f “ 1 dr-dt )  +  r2dft2  , 
we  see  that  outgoing  and  incoming  radial  null  geodesics 


may  be  represented  by  equations  u  =  const 

respectively,  where 

and 

v  =  const 

2ku-1du  =  f-1dr  -  dt 

(60) 

2kv- 1 dv  =  f_1dr  +  dt 

(61) 

and  k  is  an  adjustable  constant.  In  the 

the  line  element  (1)  takes  the  form 

u,  V 

chart 

ds2  =  ( 4k2f/uv)dudv  +  r2dft2  . 

Integration  of  (60)  and  (6l)  yields 

2  2 
r  r 

(62) 

r  +  - - -  In  1 1-r/r  |  -  — - -  In  1 1-r/r  |  =  k 

r-r  1  r  -r 

12  12 

In  |  uv 

1  ,  (63) 

t  =  k  In  |v/u|  (r  >  r  or  r  < 

i 

r  )  . 

2 

(64) 
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The  constants  of  integration  have  been  set  to  zero  for 
convenience . 

Consider  now  the  chart  u  ,v  obtained  by 

1  1 

setting  k  =  k  =  r  2/(r  -r  )  .  We  find  from  (63) 

1  1  12 

-(r  /r  )2  r  -r 

u  v  =  -  l)(p  -  1)  21  exp[— - -  r]  (r  >  r  )  , 

12  r  2 

l 

(65) 

and  (62)  exhibits  no  singularity  at  r  =  r  .  The  chart 

u  j  v  in  fact  gives  a  regular  mapping  of  any  given 

subregion  of  the  manifold  which  has  r  >  r  .  A  (coordi- 

2 

nate)  singularity  does  develop  at  r  =  r  ,  however,  and 

2 

it  is  necessary  to  go  over  to  another  chart  before  that 
happens . 


Define  the  chart  u  ,v  by  setting 

2  2 

k  =  k  =  —  r  2/(r  -r  )  in  (63)  and  (64).  Then 
2  212 


u  v 
2  2 


-(r  /r  )2 
1  2 


r  -r 

exp[ - - - 2-  r] 

r  2 
2 


(r  <  r  )  , 
1 

(66) 


and  this  provides  a  regular  covering  for  any  subregion 


with  r  <  r 

1 


In  the  domain  of  overlap  r  <  r  <  r  we  have 

2  1 


r  2  In  I u  I  +  r  2  In  I u  I  =  -  (r  2  In  I v  I  +  r  2  In  I v  |)  . 
1122  1  1  2  2 


(67) 
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Since  u  =  const  represents  an  incoming  radial  geodesic 

and  therefore  must  correspond  to  a  constant  u  ,  we  have 

2 

u  =  u  (u  )  .  Similarly  v  =  v  (v  )  .  Hence  from  (67), 
112  112 


const 


e 


-a 


,  say , 


By  adjusting  the  scales  we  can  set  a  =  0  and  obtain 

2  22  2 

r  -r  r  -r 

121  2 


1  u 

= 

1  u 

$ 

1  v  1 

= 

1  v 

1 

2 

1 

2 

(68) 

The  complete  manifold  for  e2  <  m2  is  a 

periodic  lattice  of  alternating  regions  of  type  I 

(r  >  r  )  ,  type  II  (r  <  r  <  r  )  ,  and  type  III  (r  <  r  ). 
121  2 

Figure  6  (due  to  Carter is  a  schematic  over-all 
map  with  local  singularities  at  some  of  the  lattice  points . 
Figures  7(a)  and  (b )  are  Kruskal-type  diagrams  which  to¬ 
gether  give  a  faithful  map  of  any  subregion  covered  by  a 

pair  of  overlapping  charts  u  ,v  and  u  ,v 

11  22 

Because  of  the  cyclic  character  of  the  extended 
manifold,  it  is  natural  to  raise  the  question  of  possible 
topological  identifications.  For  instance,  in  Fig.  5  for 
e2  =  m2  ,  one  might  postulate  that  all  points  (r,0+2n7r), 
n  =  0  ,  ±  1,...,  represent  the  same  physical  event. 
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Such  "space-saving"  devices  are  tempting,  but  they  lead 
to  causal  paradoxes.  In  addition,  there  would  be  dyna¬ 
mical  difficulties  connected  with  gravitational  self¬ 
interaction,  since  a  world  tube  would  then  intersect  a 
space  t  =  const  more  than  once.  These  possibilities 
will  not  be  considered  further  here. 

6.  Charged  Shell  with  e2  <  m2  in  Vacuo. 

We  now  return  to  the  discussion,  begun  in  Sect.  4, 
of  the  charged  shell  in  empty  space,  and  proceed  to  con¬ 
sider  the  exterior  view  of  the  motion  for  b  >  0  , 
e2  <  m2  ,  when  an  event  horizon  exists. 

For  a  shell  with  e2  =  m2  ,  there  is  always  a 
special  solution  ( a  =  1  in  (49)  and  (50))  which  is  static. 
The  shell  is  then  at  rest  (in  neutral  equilibrium)  at 
any  radius  R.  The  world-line  ST  (Fig.  5)  represents 
the  history  of  such  a  shell  with  R  =  const  <  m  .  The 
extended  manifold  displays  an  infinite  sequence  of  r  =  0 
physical  singularities,  e.g.  for  <  0  <  3/2  tt  .  If 

we  wish,  we  can  remove  these  singularities  and  maintain 
strict  periodicity  by  introducing  an  endless  number  of 
"re-incarnations"  of  the  shell,  e.g.  at  S'T’.  Space-time 
is  then  flat  for  r  <  R  and  all  0  .  The  result  is  of 
some  interest  mathematically,  since  it  represents  a 
universe  containing  an  event  horizon  (r=m)  which  is 
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LEGENDS  FOR  FIGURES 

Figure  5.  Schematic  representation  of  the  extended 
Reissner-Nordstrom  manifold  for  e2=m2.  Shaded  sections 
of  the  map  are  not  part  of  the  manifold.  Dashed  lines 
represent  radial  null  geodesics;  the  apparent  constriction 
of  these  lines  at  r=m  is  due  to  local  defectiveness 
of  the  coordinates.  The  time-like  curve  KLM  represents 
the  history  of  a  thin  shell,  which  implodes  in  the  space 
la,  reverses  its  motion  at  L  after  passing  through  the 
event  horizon  r=m,  then  re-expands  in  the  space  lb . 

Figure  6.  Schematic  representation  of  the  extended 

Reissner-Nordstrom  manifold  for  e2<m2.  Null  lines  are 

inclined  at  45°.  FGHJM  is  the  history  of  a  shell  which 

collapses  from. infinity  in  the  asymptotically  flat  space 

la,  passes  through  the  event  horizon  r=r  ,  comes  to  rest 

at  J  with  a  minimal  radius  smaller  than  r  ,  then  re- 

2 

expands  into  the  asymptotically  flat  space  Ic.  ABCDE  is 
the  history  of  an  oscillatory  shell  or  uniformly  charged 
sphere.  Shading  on  the  curves  distinguishes  the  interior 
domain . 

Figure  7 •  Kruskal-type  diagrams  for  portions  of  the  over¬ 
all  map  of  figure  6,  showing  the  same  curves  ABCDE  and 
FGHJM.  Figures  7(a)  and  (b)  overlap  in  the  region  lib , 
and  may  be  regarded  as  linked  together  along  the  curve  r=rQ 
where  r  is  any  convenient  value  between  r  and  r 
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Figure  5 


Figure  6 


Figures  7 
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+ 

everywhere  free  of  singularity. 

The  history  of  a  shell  with  e2  =  m2  ,  a  >  1 
is  represented  by  the  time-like  curve  KLM  in  fig.  5.  To 
an  external  observer  in  the  asymptotically  flat  space  la 
the  shell  implodes,  then  appears  to  slow  down  as  it  ap¬ 
proaches  the  observer’s  event  horizon  r  =  m  ,  reaching  it 
only  asymptotically  as  t  °°  .  On  the  other  hand,  an 
observer  moving  with  the  shell  finds  that  it  passes  rapidly 
and  uneventfully  through  r  =  m  ,  contracts  to  a  nonzero 
minimal  radius  at  L  ,  then  re-expands  into  a  new  space 
lb  ,  identical  with  la  in  its  geometrical  properties,  but 
physically  distinct  from  it.  It  appears  that  we  are  forced 
to  accept  this  resolution  of  the  paradox  encountered  in 
Sect .  4 . 


The  path  FGHJM  (Figs..  6  and  7 )  of  a  bouncing 
shell  with  b  >  0  ,  e2  <  m2  has  a  similar  general  character 
The  bounce  carries  the  shell  into  a  different  space.  A  new 
and  peculiar  feature  is  the  appearance  of  a  time-like  sing¬ 
ular  curve  r  =  0  (  the  curve  XY)  in  the  region  outside 

the  shell.  It  has  to  be  interpreted  as  the  history  of  a 
particle  with  mass  m  and  charge  -e . 


This  does  not  contradict  the  Penrose  theorem  quoted  in 
Chapter  III,  since  two  of  the  hypotheses  of  that  theorem 
are  not  satisfied  here.  In  the  first  place,  the  manifold 
with  e2  =  m2  contains  no  "trapped  surface"  (even  though 
it  contains  an  event  horizon)  ,  since  outgoing  radial  null 
geodesics  have  dr/dt  =  (1-m/r)2  >  0  and  do  not  converge 
anywhere.  Secondly,  the  manifold  with  e2  <  m2  does  not 
admit  a  Cauchy  hypersurface. 
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7 •  Test  Shell;  Uniformly  Charged  Ball  of  Dust. 


We  now  turn  briefly  to  the  situation  where  the 
hollow  interior  of  the  shell  contains  nonvanishing  charge 
e  and  mass  m  .  We  shall  confine  our  discussion  to 

i  i 

the  ease  where  the  mass  y  =  m  -  m  and  charge 

2  1 

e  =  e  -  e  of  the  shell  itself  are  small  compared  with 
2  1 

m^  and  e^  ,  and  for  a  qualitative  description  it  will  be 
sufficient  to  consider  the  limit  of  a  "test  shell" 

(y+0,  £->0  with  e/y  finite).  In  this  limit  we  easily  obtain 
from  (45)  to  (48) 

(f)2  =  -  (1-  |2+  ^1)  +  A(1  -  gS)2  (69) 


where  we  have  written  e  =  e  ,  m  =  m  . 

i  i 

As  is  to  be  expected,  (69)  agrees  with  the 

equation  of  motion  of  a  radially  moving  test  particle  of 

_y 

mass  yA  2  and  charge  e  in  the  Reissner-Nordstrom 

field  (1).  The  latter  may  be  obtained  either^^)  from  the 

(61) 

equations 


2  v  ,  a  ,  3 

d  x  +  y  v  dx 

dx  dx  dx 


(e/yA  2)Fa  — 


a 


dx 


where  the  electromagnetic  field  tensor  F^  has  the  non¬ 
vanishing  components 


e/r  2 


3 
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or  from  the  Lagrangian 

L(R,dR/dt)  =  -  yA  ^  dr/dt  +  e 4>  dx^/dt 

=  -  yA-^  { f-f ~ 1  ( dR/dt )  2  +  ee/R 

(where  4>  =  (0,0,0,e/r)  is  the  electromagnetic  vector  po- 

tential)  by  forming  the  Hamiltonian  integral  H  =  y  . 

If  j e/m |  ,  | e/y |  and  A  are  each  less  than 

unity,  (69)  shows  that  the  shell's  radius  oscillates  between 

a  maximum  larger  than  r  =  m  +  /  m2-e2  and  a  minimum 

1 

smaller  than  r  =  m  -  /  m2-e2  .  The  history  of  the  shell 

2 

is  represented  by  the  curve  ABODE...  in  figs.  6  and  7. 

In  each  oscillation  the  shell  enters  a  new  space.  As  viewed 
by  a  co-moving  or  an  interior  observer  the  oscillation  is 
strictly  periodic;  however,  the  path  ABODE...  in  the 
exterior  space-time  is  not  cyclic,  but  subject  to  a  system¬ 
atic  time-shift.  If  a  given  maximum  occurs  for  t  =  t 
in  the  space  lb  (say),  then  succeeding  maxima  (in  Ic,  etc.) 
occur  for  t  =  t  +C,t  +20,  etc.  For  large  maximal 
radius,  the  constant  C  is  nearly  equal  to  the  proper 

period  of  pulsation,  and  both  agree  closely  with  the  corres- 

* 

ponding  period  calculated  from  Newtonian  theory  . 


* 

In  the  Newtonian  description  the  pulsating  shell  of  course 
always  remains  in  the  same  space. 
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The  occurrence  of  a  bounce  is  independent  of  the 
relative  sign  of  e  and  e  ,  so  it  clearly  has  little  to 
do  with  a  contest  between  gravitational  attraction  and 
electrostatic  repulsion.  For  a  neutral  shell  (e=0)  we  ob¬ 
tain  from  (69)  by  differentiation, 

d2R/dx 2  =  -  M(R)/R2  ,  M(R)  Em-  e2/R  . 

This  brings  out  clearly  the  physical  mechanism 
responsible  for  the  bounce.  Because  the  electrostatic  field 
energy  of  the  internal  charge  e  is  diffused  throughout 
space,  less  and  less  of  it  contributes  to  the  effective  in¬ 
terior  gravitational  mass  M(R)  as  the  shell  contracts. 
Ultimately  M(R)  becomes  negative  and  there  is  a  gravita¬ 
tional  repulsion. 

Finally,  let  us  note  another  interesting  special 
case.  If  we  set  e/y  =  e/m  ,  (69)  may  be  regarded  as  the 
equation  of  motion  of  a  particle  on  the  outer  surface 
r  =  R(t)  of  a  uniformly  charged  ball  of  dust  with  total 
charge  e  and  mass  m  ,  which  is  collapsing  homologously. 
For  e2  <  m2  ,  A  <  1  ,  the  motion  is  again  oscillatory, 
and  the  history  of  the  surface  is  given  qualitatively  by 
the  curve  ABCDEh.  .  .  in  figs.  6  and  7.  This  example  has 
been  discussed  by  Novikov^^. 
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8 .  Concluding  Remarks. 

The  collapse  of  a  spherically  symmetric  body  to 
an  event  horizon  appears  as  an  irreversible  process  to  an 
external  observer.  As  we  have  seen,  the  possibility  can¬ 
not  be  ruled  out  that  the  body  reverses  its  motion  within 
the  event  horizon  and  re-expands  symmetrically.  It  then 
appears  necessary  to  believe  in  the  existence  of  other 
asymptotically  flat  spaces  geometrically  similar  to  but 
distinct  from  ours,  which  will  accomodate  the  re-expansion. 
This  seems  at  least  as  fantastic  as  the  alternative  of 
irreversible  collapse  to  virtually  point-like  dimensions. 

In  assessing  the  possible  relevance  of  these  re¬ 
sults  to  realistic  gravitational  collapse,  it  is,  of  course, 
necessary  to  keep  in  mind  the  various  idealizations  and 
hypotheses  involved  (exact  spherical  symmetry,  asymptotic 
flatness,  analytic  cont inuability  of  the  manifold,  etc.; 
see  also  Chapter  III),  each  of  which  could  be  questioned. 

The  Israel  Theorem  (Chapter  III),  which  applies 
to  fields  due  to  electrically  neutral  objects,  can  be  ex- 
t  ended  to  read: 

Among  all  static,  asymptotically  flat 
electrovac  space-times  with  closed,  simply  connected 
equipotential  surfaces  g  =  const.,  the  only  ones 
which  have  regular  event  horizon  g 

4  *t 
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Reissner-Nordstrom  family  of  spherically 
symmetric  solutions  with  m  >  |e|  . 

(The  adjective  "electrovac"  means  "devoid  of  matter  but 
containing  electromagnetic  fields").  Thus  it  appears  that 
the  implication  of  the  Israel  theorem  for  the  question  of 
asymmetric  collapse  as  discussed  in  Sect.  4  of  the  last 
Chapter  will  still  be  retained  if  charges  are  added  to  the 
collapsing  objects. 

As  a  null  hypersurface  (i.e.,  a  characteristic 
hypersurface  of  the  field  equations),  an  event  horizon  is 
a  possible  locus  of  discontinuities  of  the  field.  It  is 
not  necessary,  and  perhaps  not  physically  justified,  to 
insist  on  analytic  continuation  of  a  manifold  through  an 

.  ,  .  (63) 

event  horizon  . 

As  to  the  idealization  of  asymptotic  flatness 
for  the  collapse  of  a  stellar  mass  in  our  expanding  universe, 
it  is  justified  in  our  present  epoch,  but  clearly  not  in 
the  remote  past.  It  will  not  always  be  justified  in  the 
future  if  the  universe  happens  to  be  oscillatory.  In  fact, 
a  lattice  structure  for  space-time  of  the  general  type  we 
have  been  considering  here  would  find  a  natural  interpre¬ 
tation  in  terms  of  an  oscillatory  universe. 
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CHAPTER  V 

MACH’S  PRINCIPLE,  SPINNING  SHELL 
AS  A  SOURCE  OF  THE  KERR  METRIC1'  . 

This  chapter  is  conveniently  divided  into  two 
parts.  Part  A  gives  an  introductory  account  of  Mach’s 
principle,  with  particular  reference  to  the  theory  of 
general  relativity,  which  aside  from  its  inherent  in¬ 
terest,  will  also  prove  useful  for  an  appreciation  of  the 
subsequent  work.  Part  R  opens  with  a  brief*  derivation  of 
the  Kerr  metric,  in  the  simple  and  natural  way  recently 
discovered  by  Ernst,  It  is  then  followed  by  the  main  body 
of  the  chapter,  dealing  with  the  interiors  of  rotating 
shells  and  their  physical  properties,  on  the  presumption 
that  the  exterior  geometry  is  described  by  the  Kerr  line- 
element  . 

PART  A.  MACH’S  PRINCIPLE, 

Mach’s  idea^1^  of  Inertia  as  something  cs.  usably 
determined  by  the  distribution  of  matter  of  the  universe 
was  one  of  the  great  currents  of  thought  that  guided 
Einstein  to  the  creation  of  his  general  theory  of  relati¬ 
vity,  ^5)  -j_t  is  based  on  the  kinematic  dictum,  already 
advocated  about  one  hundred  years  earlier  than  Mach  by  the 

^  Based  on  V,  de  la  Cruz  and  W,  Israel'^', 
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Irish  philosopher  Bishop  Berkeley ^ ^ ^ ^ ,  that  the  only 
meaningful  concept  of  motion  is  that  of  relative  motions 
between  masses.  In  direct  contrast,  Newton  had  postulated 
the  existence  of  an  absolute  space,  and  in  support  cited 
his  famous  bucket  experiment*"  : 

An  empty  bucket  is  hung  along  its  axis  of  symme¬ 
try  by  a  rope.  After  the  rope  is  twisted  a  number  of  times 
the  bucket  is  held  at  rest,  half-filled  with  water,  and 
then  released.  The  bucket,  which  is  set  into  rotational 
motion  by  the  tortion  of  the  rope,  in  turn  drags  along  the 
water.  Initially,  when  both  are  at  rest,  the  shape  of  the 
water  is  a  plane.  Thereafter,  it  changes  Into  that  of  a 
paraboloid,  reaching  a  maximum  curvature  when  the  water 
catches  up  with  the  bucket,  and  finally  becomes  plane  again 

Newton’s  contention  that,  since  the  relative  angu 
lar  velocity  of  the  water  and  the  bucket  does  not  determine 
the  shape  of  the  water  surface,  one  must  and  can  meaning¬ 
fully  ascribe  to  the  water  an  absolute  motion,  was  critici¬ 
zed  by  Mach,  who  pointed  out  that  the  distant  "fixed  stars" 
in  particular,  are  all  the  while  there.  According  to  Mach, 
the  centrifugal  force  acting  on  the  water  arises  from  the 
rotation  of  the  water. with  respect  to  the  fixed  stars  or. 
Indifferently,  the  rotation  of  the  fixed  stars  with  respect 
to  the  water-.  A  priori.,  it  appears  Impossible  to  decide 


who  is  more  nearly  correct.  But  In  so  far  as  Mach’s  Idea 


ari 

■ 
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represents  an  endeavor  to  explain  in  physical  terms  -  i.e., 

without  recourse  to  the  somewhat  metaphysical  concept  of 

absolute  space  -  why  certain  reference  frames  are  in  fact 

inertial,  it  is  to  be  preferred.  Moreover,  whereas  the 
observational  fact  that  the  fixed  stars  are  indeed  fixed 
(i.e.,  non-rotating)  with  respect  to  our  inertial  frames 
cannot  be  accounted  for  on  the  basis  of  Newtonian  theory  as 
otherwise  than  accidental,  it  can  be  taken  as  an  empirical 
support  for  the  Mach's  principle. 

After  the  advent  of  the  general  theory  of  relati¬ 
vity,  works  on  Mach's  principle  have  centered  on  the  ques¬ 
tion  of  whether,  or  to  what  extent,  Mach's  principle  has  been 
incorporated  into  general  relativity.  Of  particular  impor¬ 
tance  in  this  connection  is  the  quantitative  investigation 

(21) 

of  H.  Thirring  .  Using  the  weak  field  approximations  to 
Einstein's  field  equations,  he  studied  the  motions  of  test 
particles  located  near  the  center  of  a  slowly  rotating  mass 

shell  - his  version  of  the  Newtonian  rotating  bucket. 

The- result  shows  that  such  a  particle  will  be  subjected  to 
forces  which  are  completely  analogous  to  the  Coriolis  and 
centrifugal  forces  of  classical  mechanics,  but  differ  in 
magnitudes  from  these  by  factors  of  the  order  of  GM/c2R  , 
where  M  and  R  are  the  mass  and  radius  of  the  shell, 
respectively.  Because  of  the  approximations  used,  this 
dimensionless  quantity  has  to  be  much  less  than  unity. 
Nevertheless  the  Thirring  effect  must  be  regarded  as  a  strong 


' 

. 
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manifestation  of  Mach's  principle  in  general  relativity, 
and  it  led  to  the  speculation  that,  when  strong-field  cal¬ 
culations  were  done,  so  that  the  ratio  GM/c2R  could  be 
made  more  nearly  equal  to  unity,  the  interior  inertial 
frame  would  be  rigidly  dragged  around  by  the  rotating  shell, 
in  accord  with  the  expectation  of  Mach.  This  is  indeed 
the  case  for  a  slowly  rotating  shell,  to  the  first  order  in 

(19 ) 

its  angular  velocity,  as  Brill  and  Cohen  recently 

demonstrated.  In  Part  B,  we  shall  summarize  and  extend  the 
results  of  Brill  and  Cohen,  in  conjunction  with  our  study 
of  rotating  mass  shells  as  possible  sources  of  the  Kerr 
metric . 


While  emphasizing  the  significance  of  the  work 
of  Thirring,  Einstein  himself  was  the  first  to  recognize 
that  the  gravitational  field  equations  contain  solutions 
which,  from  a  Machian  point  of  view,  would  be  unacceptable. 
Thus-,  for  example.,  the  flat  space-time  £  =  0 )  ,  while 

consistent  with  the  vacuum  field  equations-,  would  endow 
a  test  particle  with  inertia,  which,  ipso  facto ,  is  not 
caused  by  other  masses.-.  But  since  the  gravitational  equa¬ 
tions,  like  those  of  Maxwell,  are  local  differential  equa¬ 
tions,  a  large  number  of  solutions  could  be  ruled  out  by 
admitting  only  certain  boundary  conditions-,  perhaps  differ¬ 
ing  according  to  the  matter  distributions ^  so  that  Mach's 
principle  would  act  as  a  supplementary  "selection  rule". 
Einstein,  however,  rejected  this  idea  as  ad_  hoc  and  artifi- 
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eial . 
tions 
to  be 
this 
doubt 


Instead,  he  argued  that  only  spatially  closed  solu- 

-  which  do  not  require  special  boundary  conditions 

imposed  -  are  to  be  accepted.  But  whether  even 

criterion  is  sufficiently  restrictive  is  opened  to 

(68) 


Results  of  this  nature  has  indicated  to  some^^ 
as  calling  for  modifications  of  the  general  theory  of 
relativity  along  lines  more  in  accord  with  the  spirit  of 
Mach.  Others  have  revived  the  suggestion  that  Mach’s 

principle  should  serve  as  a  selection  rule  supplementary 
to  the  field  equations,  but  in  a  more  general  sense  so  as 
to  include  spatially  closed  solutions.  In  this  connection 
we  may  perhaps  remark  that,  strictly  speaking,  Mach’s 
principle  has  no  relevance  outside  the  context  of  cosmo¬ 
logy.  If  we  are  prepared  to  believe  that  the  existence  of 
our  actual  world  is  not  accidental  but  should  rather  follow 
from  some  further  synthesis  of  physical  theories,  then  the 
role  played  by  Mach's  principle  would  appear  somewhat  ob¬ 
scure.  Finally,  there  are  some  who  prefer  to  regard  Mach’s 

principle  as  a  historical  relic  -  pregnant  in  its  own 

time  with  suggestive  ideas  -  which  is  now  properly  super¬ 

seded  by  general  relativity  and  can  be  forgotten  except 
for  heuristic  purposes. 
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PART  B.  SPINNING  SHELL  AS  A  SOURCE  OP  THE  KERR  METRIC 

1.  The  Kerr  Metric :  Derivation . 

( 1 6 ) 

As  we  noted  in  Chapter  I,  the  Kerr  solution^ 
is  of  considerable  interest  because  it  probably  represents 
the  geometry  exterior  to  some  finite  rotating  object. 

This  int epret at ion  was  first  deduced  from  the  form  of  the 
metric  itself,  since  its  original  derivation  was  very  for¬ 
mal,  being  discovered  by  Kerr  "accidentally"  in  the  course 

of  his  investigation  of  algebraically-special  vacuum  fields. 

( 71 ) 

But  recently  Ernst  was  able  to  present  an  alternative 

derivation  which  employs  only  elementary  techniques  of 
analysis,  and  has  the  additional  merit  that  it  proceeds 
from  the  assumption  of  axial  symmetry  rather  than  the  assump¬ 
tion  that  the  metric  is  algebraically-special.  We  will 
give  here  an  outline  of  Ernst’s  derivation,  which  will  also 
serve  to  introduce  his  elegant  and  promising  c  -function 
method . 


The  line-element  of  a  stationary,  axially- 

symmetric  vacuum  field  can  generally  be  cast  into  the 

( 72 ) 

"canonical  form"  of  Lewis  , 

ds2  =  P_1  [e/V(dp2+dz2  )  +  p2dcj>2]  -  F(dt-^d(j) ) 2  ,  (1) 

where  the  coefficients  are  independent  of  t  (stationary) 
and  of  <J>  (axial  symmetry).  Some  of  Einstein’s  vacuum 
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field  equations  =  0  had  been  employed  in  the  reduc¬ 

tion  of  the  metric  to  the  form  (1),  and  the  rest  are 
coupled  equations  for  the  three  functions  v,  F,  and  i p  . 
It  turns  out  that,  once  F  and  ip  are  known,  v  can  be 
obtained  in  a  relatively  straightforward  manner.  The 
equations  governing  F  and  ip  are 


FV2F  -  $F  •  VF  +  p —  2 F 4 Vx/j  •  $ip  =  0  , 
V  .  (p“2F2ViH  =  0  . 


(2) 

(3) 


Here,  the  three-dimensional  gradient  operator  is  to  be 
understood,  treating  (p,z,4>)  formally  as  the  usual  cylin¬ 
drical  coordinates  in  ordinary  Euclidean  3-space.  For 
v  ,  we  have  (subscripts  indicating  partial  differentiations) 


v 

P 


hp  1el*X(ipp2-\ pz2) 


(4) 


v 


z 


2pA  A 
K  p  z 


—  ^2  P 


-  1 


‘♦A 

e  ip  ip 
p  z 


(5) 


(where  e2v  e  F)  of  which  (2),  (3)  are  the  integrability 
conditions . 

A 

Let  denote  the  unit  normal  vector  in  the 

azimuthal  direction  and  Q  be  any  reasonable  function. 
Then  by  virtue  of  the  identity 

A 

V  •  (p_1  <j>  x  V  Q)  =  0 


(6) 
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eq .  (3)  will  be  automatically  satisfied  if  we  set 


p“  1 F 2 Vip  =  <p  x  VQ 


(7) 


Provided  ip  is  independent  of  <j>  (which  we  now  assume), 
(7)  is  equivalent  to 

P-2^Q  =  -  p_1  <J>  x  Vip  , 

and  hence  the  identity  (6)  implies  the  field  equation 


V  •  (F“2VQ)  =  0 


(8) 


for  the  new  potential  Q.  When  (2)  is  expressed  in  term 

of  Q  one  then  finds  that  the  complex  function  (" c -function") 


£,  =  F  +  i  Q 


(9) 


satisfies  the  simple  homogeneous  quadratic  differential 
equation 

(Re  t  )V2£  =  %t-  v£  .  (10) 

A  modification  of  the  ^ -equation  (10),  convenient  in 
certain  cases,  is  obtained  by  the  substitution 

£  =  u-n/u+n  ,  (in 

whereby  (10)  becomes 

(C£*-1)V2£  =  2£*VC  •  VC  .  (12) 

We  now  particularize  our  discussion  to  the 
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derivation  of  the  Kerr  metric.  For  this  purpose,  it  is 
advantageous  to  introduce  "prolate  spheroidal"  coordinates 
(x,y)  ,  where 

p  =  (x2-l)^  ( 1-y 2  )?s  ,  (13) 

z  =  xy  . 

In  terms  of  x,  y,  we  have  the  expressions 


V2A 


1 


(i-y2) 


[  (x2-l) 


3A 

9x  9x 


(l-y2) 


9A  3B, 

3y  3yJ  5 


A  , 

(14) 


for  arbitrary  functions  A,  B  independent  of  the  azimuth. 

With  the  help  of  (14),  it  is  easy  to  verify  that 
K  =  x  is  a  particular  solution  of  (12).  But  since  the 
operations  (14)  are  symmetric  in  x  and  y,  £  =  y  is  also 
a  solution.  If  we  now  look  for  a  linear  combination  of 
E,  =  x  and  £  =  y  which  will  again  satisfy  (12),  one 
arrives  at 

£  =  x  cos  a  +  iy  sin  a  ,  (15) 


This  solution  gives  rise  to  the  Schwarzschild  solution, 
e . q . (11.21) .  (F  =  (x-l)/(x+l),  Q  =  0,  e2v  =  (x2-l)/(x2-y 2 ) ; 
r  =  x+1  ,  cos  0  =  y  ;  lengths  measured  in  units  of  m)  . 
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where  a  is  a  real  parameter.  Solution  (15)  corresponds 

to  the  Kerr  solution.  To  obtain  the  metric  in  its  familiar 
(73) 

form  ,  introduce  a  =  tan  a  and  m  =  sec  a  ,  noting 

p 

that  lengths  are  measured  in  units  of  (m2-a2)2 

[since  sec2  a  -  tan2  a  =  1]  ,  and  identify  the  radial  and 

polar  coordinates  by 

r  =  x(m2-a2)^  +  m  ,  cos  0  =  y  .  (16) 

When  the  entire  metric  is  constructed  one  obtains 
ds2  =  (r 2  +  a2 cos  2 0 ) [dr 2 /( r 2 -2mr+a2 )  +  d02]  + 

+  [r2+a2+2mra2sin2 0/(r2+a2cos20 ) ] sin20  d<f> ^  + 

+  [4amr  sin2  0/(r  2+a2 cos  2  0  )  ]  dcj)dt  - 
-  [l-Pmr/(r2+a2cos2 0 ) ] dt 2  ,  (17) 

which  is  the  Kerr  line-element. 


2 .  Spinning  Shell  as  a  Source  of  the  Kerr  Metric : 
Introduction . 


The  Kerr  metric  (17)  is,  as  we  have  just  seen,  a 
2-parameter  particular  solution  of  Einstein's  vacuum  equa¬ 
tions  which  is  axially  symmetric  and  stationary  (but  not 
static  :  g  ^  ^  0  unless  a  =  0  ,  In  which  case  it  reduces 
to  the  Schwarzschild  line-element).  One  can  also  easily 
verify  that  it  is  asymptotically  MInkowskian. 
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In  a  way  reminiscent  of  the  cases  with  the 
Schwarzschild  and  the  Reissner-Nordstrom  solutions,  the 

Kerr  line-element  in  the  form  given  is  regular  only  for 

h 

r  >  m  +  (m2~a2) 2  when  |a|  <  m  .  Its  maximal  analytic 

(71) 

completion,  found  by  Boyer  and  Lindquist  ,  turns  out 
to  be  vastly  more  complicated  than  the  other  two  cases. 
However,  for  our  purpose  here,  only  the  regular  part  of 
(17)  will  be  needed. 

The  Kerr  metric  is  generally  believed  to  repre¬ 
sent  the  exterior  gravitational  field  of  an  isolated, 
uniformly  spinning,  spherical  (or  nearly  spherical)  body 
of  mass  m  and  angular  velocity  proportional  to  (-a). 
Until  quite  recently,  the  sole  basis  for  this  belief  was 
the  agreement  between  the  asymptotic  form  of  (17)  and  the 

weak  gravitational  field  of  such  a  body,  computed  from  the 

(  74 ) 

linearized  Einstein's  equations  : 

ds2  ~  (l+2m/r)dr2  +  r2d<72  -  (4J/r)sinz0  d(j)dt  - 

-  ( l-2m/r ) dt  2  , 

where  J  is  the  angular  momentum.  Despite  several  at- 

( 75 ) 

tempts  and  partial  results  ,  no  one  has  yet  been  able 
to  match  the  Kerr  metric  to  an  exact  interior  solution. 

However,  a  significant  advance  towards  this  goal 

( 19 ) 

was  achieved  recently  by  Brill  and  Cohen  .  These 
authors  derived  an  approximate  expression  for  the  field 


. 
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of  a  slowly  spinning  thin  spherical  shell  by  considering 
it  as  a  linearized  perturbation  of  the  static  spherically 
symmetric  solution  in  which  the  angular  velocity  oo  is 

treated  as  a  small  parameter,  but  without  restrictions  on 

(17) 

the  mass  or  radius  of  the  shell.  Cohen  was  then  able 

to  show  that  the  exterior  field  agrees  with  the  Kerr  metric 
to  first  order  in  a.  Thus,  if  a  possible  source  of  Kerr's 
metric  is  taken  to  be  a  thin  shell  held  together  by  a  sur¬ 
face  pressure  and  hoop  stresses,  then  it  follows  from  the 
results  of  Brill  and  Cohen  that  (to  first  order  in  a) 

(i)  the  shell  is  spherical  and  of  uniform  density,  (ii)  it 
spins  rigidly  (independent  of  latitude),  and  (iii)  the 
interior  space-time  is  flat.  The  last  result  means  that 
inertial  frames  can  be  defined  globally  in  the  interior. 

In  the  limit  when  the  radius  of  the  shell  approaches  its 
Schwarzschild  radius,  Brill  and  Cohen  found  that  the  inter¬ 
nal  inertial  frames  are  dragged  around  rigidly  by  the  shell. 
This  conclusion  appears  to  support  the  Machian  philosophy 
underlying  the  investigations  of  Thirring ,  as  discussed  in 
Part  A. 


It  is  of  considerable  interest  to  find  out  whether 
these  elegant  results  continue  to  be  valid  beyond  the  first- 
order  approximation  of  Brill  and  Cohen.  In  the  present 
paper  we  shall  describe  a  systematic  iterative  procedure  for 
developing  the  interior  field  of  a  spinning  shell  as  a 
power  series  in  m  on  the  assumption  that  the  exterior  field 
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is  the  Kerr  metric.  Our  method  -  described  in  more  detail 

in  Sect.  3  -  is  to  connect  the  Kerr  metric  continuously 

across  a  nearly  spherical  boundary  £  (whose  shape  can  be 
chosen  arbitrarily  to  within  quantities  of  order  a2)  to  an 
interior  axisymmetric  vacuum  solution  which  is  developed  by 
successive  approximations  in  a  power  series  in  a”1  „  The 
physical  characteristics  of  the  shell  are  then  interred  a 
posteriori  from  the  discontinuity  in  the  normal  derivative 
of  the  metric  tensor  across  £  .  In  lowest  approximation 
(Sect.  4)  we  merely  recover  the  results  of  Brill  and  Cohen 
in  a  new  and  very  simple  way.  The  next  approximation  (Sect. 
5)  gives  results  correct  to  order  aJ  and  again  requires 
only  simple  calculations,  because  departures  from  flatness 
in  the  interior  are  linear  to  this  order.  Details  of  calcu¬ 
lations  can  be  found  in  the  Appendix, 

We  shall  find  that  the  elegant  or  .  r ties  (i), 

(ii)  and  (iii)  all  break  down  in  the  next  approximation. 

But  in  the  (strong  field)  limit  of  compression  to  the 
Schwarzschild  radius,  a  "Machian"  feature  survives.  The 
rotation  tends  to  a  rigid  motion  In  this  limit ,  and  -  pro¬ 

vided  the  shape  of  £  is  suitably  chosen  — —  the  interior 
metric  tends  to  flatness,  with  the  internal  inertial  frames 
rigidly  attached  to  the  spinning  shell. 
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3.  Outline  of  Method. 


We  recall  (see  Chap.  IV)  that  thin  shells  in 
general  relativity  may  be  characterized  in  a  geometrical 


way  as  follows. 


(57) 


Consider  two  distinct  space-time 


r+ 


manifolds  V  and  V  which  are  partially  bounded  by 
time-like,  cylindrical  hypersurf aces  E-  and  E**  respec¬ 
tively.  We  could  describe  V-  and  V+  as  "compatible" 


if 


and 


.+ 


can  be  put  into  isometric  correspondence. 


In  that  case,  by  making  the  identification  E  =  E  (  =  E), 

—  -I- 

we  can  consider  V  and  V  together  as  comprising  a 

single  space-time  which  contains  a  (generally  singular) 

hypersurface  E  .  This  hypersurface  has  to  be  interpreted 

as  the  history  of  a  thin  shell  or  surface  layer  if  the 

extrinsic  curvatures,  K-,  and  K+,  ,  of  E~  in  V-  and 

5  ab  ab  } 

+  +  * 

E  in  V  are  unequal  .  The  surface  energy  tensor  S  ^ 

of  the  layer  is  given  by  the  Lanczos  jump  conditions1  7^', 

(  57 ) 

which  can  be  written  in  the  form 


-  8ttS  ,=Y,-g,Y 
ab  1  ab  toab  ’ 


(18; 


where  y  ,  =  K  *  -  K 

clD  clu 


ab  5 


Y  =  g 


ab 


ab 


We  shall  take  V  to  be  Kerr’s  vacuum  manifold 
with  metric  (17).  For  the  hypersurface  2^  we  could 
choose  quite  generally  the  equation  r  =  ( arb „ const . )  +  a‘x 
(arb.  function  of  6  with  equatorial  symmetry).  For  our 


Latin  indices  range  from  2  to  4  and  distinguish  3-tensors 
defined  on  E  .  Greek  indices  range  from  1  to  4. 


, 
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present  purposes  (.explicit  calculations  not  carried  beyond 
a3)  it  will  be  sufficient  and  convenient  to  assume  more 
specifically 

Z  :  r  =  R  (1+ek  V2  cos2  0)  (19) 

where  R  and  k  are  arbitrary  constants,  and 

V2  =  1  —  2m/R  +  a2/R2  ,  e  =  a2/'R2  .  (20) 

(71)  + 

If  V  =  0  ,  it  can  be  shown  that  E  is  a  regular  null 

+ 

hypersurface  of  V  (real  for  |a|  <  m) ,  the  most  natural 

t 

analogue  of  the  critical  hypercylinder  r  =  2m  in 
Schwarzschild ’ s  space-time,  to  which  it  in  fact  reduces  when 
a  =  0 . 


For  the  interior  space  V  ,  we  write  its  metric 
in  the  canonical  form  for  axi-symmetric  stationary  vacuum 
fields ,  eq . ( 1 )  : 

( ds 2 ) _  =  e2 ( V_A  }  (dp2  +  dz2 )  +  p£e~"Ad$2  - 

-  e2A(dT-T|jd<!)  2  ,  (21) 

in  which  the  coefficients  are  functions  of  p  and  z. 

The  vacuum  field  equations  for  this  line-element,  essential¬ 
ly  (2)-(5)a  are  (subscripts  indicate  partial  differentiation) 


t 


But  a  generalized  ^Schwarzschild  surface"  can  also  be 
defined  in  other  ways.  Compare  C.  V.  Vishveshwara ,  ref.  50, 
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■  p  (x;- 

Az>  -  %  P  !e',X(4'p2~'h2)  ’ 

(22a) 

=  2p  A  X 
p  z 

-  h  P~  e‘tX\ppilJz  , 

(22b) 

+  p”  A 

P 

+  xzz  =  -  h  p“2e4X(i| l>p2+tyz)  , 

(23a) 

a 

-3- 

1 

a 

l 

+  ^zz  =  "  4(Vp  +  Vz}  * 

(23b) 

Equations  (23a, b)  are  the  integrability  conditions  of 
(22a, b) . 


Adopting  51  =  (0  ,<j>,t)  as  intrinsic  co-ordinates 
of  E  ,  we  postulate  that  the  parametric  equations 
xa  =  xa  (  0  ,  cf> , t )  of  E  have  the  form 


E  :  p  =  f  (9)  ,  z  =  f  (0)  ,  (24) 

i  2 

$  =  <f>  +  aftt  ,  T  =  At  ,  (25) 

where  A  and  ft  are  constants. 


Since  we  expect  the  complete  solution  to  mirror 
the  symmetries  of  the  exterior  Kerr  metric  -  i.e.  invar¬ 

iance  under  reflection  in  the  equatorial  plane  ( 0  tt  -  0 ) 
and  under  reversal  of  the  sense  of  rotation  ($->--$,  a-*-a) — — 
it  follows  that  all  quantities  entering  (21),  (24),  (25) 

have  equatorial  symmetry  (except,  obviously,  for  f  (0)) 

2 

and  are  even  functions  of  a  (except  for  \p  ,  which  is  odd). 

We  have  now  to  impose  the  compatibility  require¬ 
ment  that  the  intrinsic  metric  of  E  ,  computed  from 
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(21),  (24)  and  (25)  should  agree  with  that  of  Z+  , 

computed  from  (17)  and  (19).  Equating  the  coefficients  of 

d02,d02,dt2  and  d0dt  furnishes  four  conditions  for  four 

functions  of  0  ,  namely  f  (0),  f  (0)  and  the  values  of 

1  2 

X  and  i/j  on  Z  ,  and  also  fixes  the  constants  A  and  Q 
Equations  (23a, b)  then  in  principle  determine  X  and  ip 
throughout  V  as  solutions  of  an  elliptic  boundary-value 
problem.  The  function  v  is  obtainable  by  a  quadrature 
from  (22*b)  together  with  the  requirement  of  "elementary 
flatness"^^  (v=0  on  the  axis). 

In  practice,  the  solution  of  (23a, b)  of  course 
presents  difficulties.  It  can  be  reduced  to  a  linear 
boundary-value  problem  of  standard  type  if  we  proceed  by 
successive  approximations,  treating  the  nonlinear  right 
sides  of  (22)  and  (23)  as  perturbations  whose  values  are 
assumed  known  from  the  previous  approximation.  Since  all 
quantities  in  (22)  and  (23)  are  small  for  small  a  (a=0 
corresponds  to  a  static  spherical  shell  with  a  flat  inter¬ 
ior)  ,  our  procedure  may  be  regarded  as  an  expansion  in 
powers  of  e  =  a2/R2  . 

As  soon  as  an  expression  for  the  interior  metric 
is  known,  it  is  straightforward  to  compute  the  extrinsic 
curvatures  from  the  defining  relations  (see  eq.(ll)  of 
Chap.  IV) 

K  (T  =  -  n  '  {  3  2xa^  ^  5b  +  r“  (3  xA/3  £a)  (3  Xy/^  £b)  } 

a.D  Oi  A  |J 


(26) 


■ 


' 
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in  which  n  is  the  outward  unit  normal  to  £  .  The 
a 

surface  energy  tensor  of  the  shell  is  then  read  off  from 

( 1 8 ) .  The  velocity  ua  and  the  proper  surface  density 

+ 

a ( 9 )  are  defined  by  the  eigenvalue  equation 

Sa  ub  =  -  oua  ,  uaua  =  -  1  .  (27) 

It  can  be  shown  easily  by  symmetry  argument  that  u  =  0, 
and  we  have  u^  =  um^  ,  where  oj(0)  =  d<j>/dt  is  the 
angular  velocity  of  a  shell  zone  as  measured  by  a  station¬ 
ary  observer  at  infinity.  According  to  ( 1 8 )  and  (27)  , 
a)  is  determinable  from 

"  =  y  t  /  ("Yj+Yt-T*)  •  (28) 

4 .  First  Approximation. 

In  lowest  approximation,  we  neglect  a2  .  To 
this  order,  Kerr's  metric  (17)  induces  the  intrinsic 
metric 

(ds2)^  =  R2  (d0  2+sin2  0d<J) 2 )  +  (  4ma/R)  sin20d<J)dt  - 

-  ( l-2m/R) dt 2  (29) 

on  the  shell  r  =  R+0(e).  One  verifies  immediately  that 
this  agrees  with  the  metric  induced  on  the  hypersurface 


t 


Cf.  eq.  (18)  of  Chap.  II. 


>U  J  0 

■ 
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£  :  p  =  R  sln8  +  0  ( e )  ,  z  =  R  cos9  +  0(e)  (30) 

by  a  flat  interior  line-element  dp2  +  dz2  +  p2d$2  -  dT2 
provided  we  choose  the  values 

n  =  2m/R3  +  0(e)  ,  A  =  (l-2m/R)^  +  0(e)  (3D 

for  the  constants  in  (25).  Thus,  in  lowest  order 5  the 
boundary  conditions  on  £  are  satisfied  by  the  trivial 
null  solution  of  (22)  and  (23);  we  have 

A,v  =  0(e)  ,  ip  =  a0(e)  .  (32) 

Computation  of  the  extrinsic  curvatures  and 
substitution  in  (18)  yields  the  following  nonvanishing 
components  of  the  surface  energy  tensor  (neglecting  e  ) : 

S®  =  sj  =  [R(l-V)-m]/8TrR2V  ,  (33) 

s£  =  -  (  1-V)/4ttR  ,  (34) 

S^  =  ma(  1+2V) /8ttR4 V  ,  S^  =  -  3ma  sin2  6/87TR2V  (35) 

The  derived  surface  density  and  angular  velocity 


a  =  (1-V)/4ttR  , 

(36) 

2am  ( 1+2V ) 

'  (37) 

R3  ( 1-V) ( 1+3V) 

are  independent  of  0  in  lowest  approximation.  According 


See  Sect.  4  of  Appendix  for  details. 


; 

. 


Ill 


to  (33),  the  shell  is  held  in  equilibrium  by  a  uniform 
surface  pressure.  (Effects  of  the  hoop  stresses  and 
departures  from  sphericity  depend  on  the  square  of  oj  and 
do  not  enter  until  the  next  approximation.)  Equations  (36) 
and  (37)  agree  with  results  obtained,  using  a  different 

(IQ  17) 

approach,  by  Brill  and  Cohen.  9 

The  inertial  frames  in  the  flat  interior  rotate 
with  angular  speed  (from  (25),  (3D) 

(d<j>/dt).  ,  =  -  afi  =  -  2am/R 3  ,  (38) 

r  $=const.  9 

as  measured  by  a  stationary  observer  at  infinity.  This  is 
numerically  smaller  than  the  angular  velocity  oo  of  the 
shell  itself,  as  measured  by  the  same  observer,  but  becomes 
equal  to  it  in  the  limit  R  -*  2m  (i.e.  V-*0). 

From  the  point  of  view  of  an  inertial  observer 
inside  the  shell  the  inertial  frames  at  infinity  spin  at 
the  faster  rate 

(d*/dTh=const.  “  an/V  ,  (39) 

which  tends  to  infinity  when  R  -*  2m.  The  asymmetry 
between  (38)  and  (39)  is,  of  course,  due  to  the  effects  of 
the  gravitational  redshift.  The  angular  velocity  of  the 
shell  relative  to  this  interior  observer  is 

a)'  5  (d$/dTD  ,  =  VDw+afi) 

<f)-oot  =  const . 


112 


6ma  V 

-  • 

R3  ( 1-V) ( 1+3V) 

Prom  (40)  we  see  that  w  ’  -*  0  when  R  2m  , 
that  a  rigid  Machian  dragging  of  the  interior 
frame  occurs  in  this  limit. 


(40) 

which  means 
inertial 


5 .  Second  Approximation. 


For  the  next  approximation,  which  neglects  terms 

beyond  order  a3  ,  we  see  from  (32)  that  the  right  sides  of 

(22)  and  (23)  can  be  replaced  by  zero.  Thus,  v  has  to  be 

a  constant  -  in  fact  zero,  to  ensure  "elementary  flat¬ 

ly) 

ness"  on  the  axis  of  symmetry  -  and  the  general  solu¬ 

tions  for  X  and  ip  in  the  interior  have  the  form 


X  =  a* 1 2  [c  +c  ( 2z2-p 2 )  +  ...]  , 

1  2 

ip  =  a3  [c  ?+c  ’  p2  +  c  '  (p4-4p2z2  )  +  ...]. 

1  2  3 


(41) 


The  line-element  (21)  is  invariant  in  form  under 
a  constant  scale  transformation  of  p  and  z  ,  and  we  can 
use  this  freedom  to  arrange  c  =  0  .  Similarly,  the  trans- 

i 

formation 


$  -*  $  +  a3cTT  ,  \p  ip  -  a3c'p2 
2  2 

leaves  (21)  invariant  to  order  a3  ,  and  effectively  makes 

c’  =  0  . 

2 


' 
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The  remaining  constants  in  ( 4 1 )  and  the  functions 
of  order  e  involved  in  (30)  and  (31)  have  now  to  be  deter¬ 
mined  from  the  isometry  of  E~  and  Z+  to  order^"  a3  . 
Equating  coefficients  of  dt2  in  the  two  line-elements  and 
employing  (30)  and  (31)  where  this  entails  only  errors  of 
order  e2  ,  we  find 


X  = 

3  (ma2/R5) (1+k) (2z2-p2 )  , 

(42) 

II 

CM 

< 

V2  [l-e(l+y)  +  3  ey(l+k)]  , 

(43) 

where  y  =  2m/R. 

From  the  coefficients  of  d<j)2 

and 

d02  we  derive  in  turn 

2  :  p  =  R  sin  0[l+^e  +  ^ey ( 2-k ) +ek ( 1-hu ) cos 2 6 ]  ,  (44) 

z  =  R  cos  0[1-  ^  ey ( 1+k )  +  ek ( 1-^y ) cos  2 0 ]  .  (45) 

Finally,  equating  coefficients  of  d(J)dt  yields 

ip  =  |  (ma3/R7)(l+3k)V(p4-4p2z2)  ,  (46) 

=  ( 2m/R 3 ) [ 1-  |  e ( 2+k ) +  i  ey(3k-4)]  .  (47) 

b  5 

The  interior  solution  is  thus  determined  in  terms 

of  four  disposable  parameters  m,  a,  R  and  k.  The  surface 

energy  tensor  and  physical  properties  of  the  shell  can  now 

+ 

be  read  off  from  the  extrinsic  curvatures  K  ,  obtained 

ab 

from  (26)  and  listed  in  the  Appendix.  Since  the  general 


t 


See  Sect.  1  of  Appendix  for  details. 


-  .  -P  :  5  ■ 
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formulae  are  somewhat  unwieldy  ,  we  shall  here  quote 
explicit  results  only  for  the  two  extreme  limiting  cases. 

In  the  weak-field  limit  m/R  ->  0  we  find 
4fTR2a/m  +  1  +(l/12)e(8k-13)  -  e ( 4k+13/4 ) cos 2 6  ,  (48) 

w  -*■  -  (  3a/2R2 )  [  1+  £  (  3-l6k )  +  e  ( 5/12-2k)  cos 2  6  ]  ,  (49) 

o  u 

showing  that  the  distribution  of  proper  mass  and  the  angu¬ 
lar  velocity  are,  in  general,  latitude-dependent. 

At  the  opposite  extreme,  letting  R  approach  the 
"gravitational  radius"  m  +  (m2-a2)2  (i.e.  letting  V  ->  0) 
we  find 

4ttRct  =  1  +  i  e  ( 2-k)  +  he  (k-5)  cos 2  0  +  0(V2)  ,  (50) 

m  =  -  a/(R2+a2 )  +  0(V2)  .  (51) 

The  rotation  is  rigid  in  this  limit.  (The  result  (51) 
actually  holds  exactly  to  all  orders  in  a.  This  is  easily 
proved  by  noting  that  the  time-like  vector  u  tangent 
to  Z  must  become  parallel  to  the  generators  of  the  null 
hypersurface  V  =  0  of  the  Kerr  manifold,  since  every 
other  direction  in  this  hypersurface  is  space-like.)  We 
observe  also  that  the  surface  density  a  remains  finite 
and,  in  general,  is  latitude-dependent.  However,  in  this 
limit,  it  is  not  the  density  that  plays  the  dominant  role, 
but  the  surface  pressure  and  the  hoop  stresses,  which 


■ 


< 
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become  infinite  (see  (33)  and  the  formulae  of  the  Appendix), 

and  which  contribute  all  of  the  gravitational  mass  m.  This 

* 

is  immediately  evident  from  the  formula 

m  =  -  /  (-g44)"^  (S4-S* 2-S3)dE  ,  (52) 

^  4  2  3  2 

2 

which  holds  generally  for  the  stationary  field  of  any  shell 
in  a  condition  of  steady  motion  in  vacuo. 

6 .  Discussion . 

We  have  not  entered  into  the  question  of  the 
convergence  of  our  approximation  procedure.  But  if  conver¬ 
gence  is  taken  for  granted,  we  may  infer  the  existence  of 
an  infinite  variety  of  exact  vacuum  solutions,  which  repre¬ 
sent  the  interior  fields  of  spinning  shells  of  differing 
degrees  of  ellipticity,  and  which  all  join  continuously  to 
the  Kerr  metric.  At  the  same  time  our  results  indicate 
that,  if  any  such  exact  solution  is  found  in  closed  form, 
it  is  unlikely  to  be  particularly  simple  or  attractive, 


Intrinsic  coordinates  £2 ,  £3 ,  £4  =  t  are  presupposed 

in  which  the  time-like  killing  vector  of  the  field 

(which  is  also  an  intrinsic  killing  vector  of  E  )  has 

components  n  =  (0,  0,  1).  The  integration  is  over  a 

two-dimensional  slice  t  =  const  of  E  (dE  is  an  invar- 

2 

iant  element  of  area).  For  a  general  proof  of  (52)  and 
of  a  corresponding  formula  for  the  angular  momentum,  see 
V.  De  la  Cruz  and  W.  Israel,  ref.  64,  Appendix  2. 
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since  for  no  choice  of  the  ellipticity  constant  k  is  the 
shell  both  rigid  and  uniform  (except  in  the  limit  V  •*  0). 

In  this  sense,  a  spinning  shell  does  not  appear  to  be  a 
"natural"  source  for  the  Kerr  metric. 

We  note  from  (42)  that  the  special  choice  k  =  -  I 
makes  A  =  0  .  It  follows  that  the  Interior  space-time 
tends  to  flatness  when  V  0  ,  since  ip  goes  to  zero  with 
V  according  to  (46).  From  (47)  and  (51)  we  find  again 
that  oo  ’  -*  0  (cf.(40)).  Thus  the  "Machian"  effect  obtained 
in  Sect. 4  persists  at  least  to  third  order  in  the  angular 
velocity.  But  it  should  be  remarked  that  the  credibility 
of  our  limiting  result  as  an  armchair  version  of  Newton’s 
bucket  experiment  is  somewhat  impaired  by  the  circumstance 
already  noted  that,  as  V  ->  0  ,  the  surface  pressure  be¬ 
comes  infinite  and  enormously  exceeds  the  density.  Of  course, 
a  model  for  Newton’s  bucket  less  extremely  idealized  than 
an  infinitely  thin  shell  should  at  least  ameliorate  this 
difficulty . 


We  have  throughout  confined  attention  to  shells 
whose  exterior  field  is  of  the  Kerr  type.  The  Kerr  solu¬ 
tions  actually  form  a  subset  of  measure  zero  in  the  class 
of  all  stationary,  asymptotically  flat  vacuum  fields.  But 
they  are  probably  the  only  fields  in  this  class  which 


' 
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*4* 

possess  non-singular  stationary  null  hypersurfaces 
(g*4")-'  =  o  .  If,  therefore,  one  Is  primarily  interested 
in  the  limit  of  compression  to  the  gravitational  radius, 
then  it  is  worth  noting  that  our  discussion  has  probably 
been  exhaustive  in  this  respect. 


Prom  the  formulas  for  small  stationary  vacuum  perturba¬ 
tions  of  the  Schwarzschild  manifold  given  by  Ernst 
(ref.  71) 3  it  can  be  seen  that  the  only  such  perturba¬ 
tions  which  preserve  asymptotic  flatness  and  a  regular 
event  horizon  are  members  of  the  Kerr  family „  See  also 
W.  Israel,  ref.  62. 
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APPENDIX 

This  appendix  pertains  to  Part  B  of  Chapter  V, 

The  following  quantities  and  expressions  (correct  to  0(  a") 
unless  otherwise  stated)  are  here  derived  in  some  detail: 
interior  metric  of  the  spinning  shell  (Sect,  1);  extrinsic 

it  + 

curvatures  K  ^  of  the  shell  with  respect  to  V  (Sect,  2) 
and  V  (Sect.  3);  shell  stress-energy  tensor  S  ,  proper 

cl  D 

energy  density  a  ,  and  angular  velocity  to  ,  correct  to 
order  a  (Sect.  4);  and  shell  energy  density  and  angular 
velocity  to  order  a3  (Sect.  5). 

We  introduce  the  following  notations: 

a  =  a/R  , 


0  e  0(an)  . 
n 


1 .  Interior  Line-element  of  the  Spinning  Shell,  Correct 
to  0 


v  = 


The  interior  line-element  is  given  by  (V.21)  with 
0  and  X  ,  ifj  of  the  forms  (V.4l).  Thus, 


(ds2)  =  ( 1-2X  )  ( dp  z +dz  2 +  p  dO^  )  - 


-  ( 1+2  X )  dT2  +  2i jj  d<MT  , 


(1A) 


Unless  otherwise  stated,  all  equations  in  the  Appendix 
are  correct  to  0 


3 


' 
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A  =  a2  [c  +c  (2z2-p2 )  +  • - • ]  (IB) 

1  2 

ip  =  a 3  [c  '  +c  ’  p 2  +c  '  (  p  ‘ z2  )+  •  •  *  ]  .  (1C) 

12  3 

Prom  (V.24,25)  and  the  first  order  results 
(V. 30,31) ,  we  see  that  the  parametric  equations  of  Z  in 
V  have  the  forms 

Z  :  p  =  R  sin  0  [l+a2f(0)]  , 

z  =  R  cos  0  [l+a2g(0)]  , 

(2) 

$=({>  +  afit  , 

T  =  At  , 

with 

n  =  2m/R 3  +  a2£/R2  ,  (3A) 

A2  =  V2 ( l+a2h/R2 )  ,  (3B) 

for  some  functions  f  ( 0 )  ,  g(,0)  and  constants  i  ,  h. 

Using  (2),  we  find. that  the  interior  metric 
(1A)  induces  the  following  intrinsic  metric  tensor  on  Z  : 

3g QQ  =  Rz tl^2 A+2a2 ( f  cos  2  0  +  f ’  sin0  cosO  + 

y  y 

+  g  sin20  -  g’  sin0  cos0)]  ,  (4A) 

3g~^  =  R2  sin2  0(l+2a2f)  -  2aR^  sinz  0  ,  (4b) 
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3S,j>t  =  afi  +  ^ 

3g~t  =  -  [  A  2  ( 1  +  2  A  )  -  a2R2  sin2  0  ( y/R2 ) * ]  ,  (4D) 

(other  components  vanish)  where  prime  stands  for  d/d0 
and  a  ,  ijj  are  evaluated  on  the  shell. 

Prom  (V.19)  , 

dr2  =  0  d02  on  E  » 

k 

The  intrinsic  metric  induced  by  the  exterior  (Kerr)  line- 
element  (V.17)  is  (|j  e  2m/R) 

'Sqq  =  R2  [1+a2 cos  2 0 ( 2kV2  +  1 ) ]  ,  (5A) 

3gJ(J)  =  R2  sin2  0  [l+a2(l+2kV2cos20+y  sin2  0)]  ,  (5B) 

=  R  ya  sin2 0 [l-a2cos' 0 (k  V2+lj]  *  (5C) 

3g^t  =  -  [l-y+a2y ( 1+kV2 ) cos2 0 ]  ,  (5D) 

(other  components  zero).  For  later  use,  we  list  here  the 
corresponding  expressions  for  the  non-vanishing  contra- 


variant  components: 

3g®^  =  R-2  [1-a2 ( l+2kV2 ) cos 2 0 ]  ,  (6A) 

3g^  =  (RV  sin  0 ) “2  [1-y  +  a2 cos  2 0 (y-2kV4 ) ]  ,  (6B) 

3gJt  =  ay(RV2 )_1  [l-a2cos20(l  +  3k-2ky)]  ,  (6C) 

3g^t'  =  -  V2  [1+a2  ( 1+y  sin20-ky  cos20)]  .  (6D) 


: 
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We  now  impose  the  condition  3g^,  =  3g  K  .  For 

cLD  cl  D 

(a,b)  =  (4,4),  eqs.(4D)  and  (5D)  yield,  with  the  help  of 
(3B)  , 


X  =  -  h  a2 [ ( 1+h-ky ) +y ( 1+k ) sin2 0 ]  on  £  , 

or,  defining  a  constant  3  by 

-y  3  =  1  +  h  -  ky  ,  (7) 

X  =  h  a2 y [3- ( 1+k) sin^ 0 ]  on  £  ,  (8) 

For  (a,b)  =  (3,3),  we  find,  when  (8)  is  taken  into  account, 

f  (  0  )  =  R~2  [  1+y  (  3-k) +k  (  2-y  )  cos 2  0  ]  +  0  .  (9) 

2 

Equating  the  ( 3 , 4 ) -component s  yields 

ijj  =  a3y  VR  sin2  0  [  6 -3h+ ( l+3k)  sin'1  0  ]  on  £  ,  (10) 

where  the  constant  6  is  defined  by 

Z  =  -  y  R~2 ( 2+V2 6 )  ,  (11) 


Z  being  the  unspecified  constant  in  (3A). 


Finally,  the  condition  3gQ0  =  3g00 
differential  equation  for  g(0): 


gives  rise  to  a 


0  =  y ( l+k-3 ) sin2 0-3k ( 2-y ) sin2 0  cos^0  + 

2 

+  2R2g  sin^0  -  2R2g'  sin0  cos0  , 
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whose  solution  is  easily  found  to  be 

2R2g(0)  =  k(2-y)cos20-y(l+k-0)  +  0  ,  (12) 

2 

in  which  the  constant  of  integration  has  been  set  to  zero. 


We  now  determine  the  constant  coefficients  in 

(1B,C)  for  X  and  ifj  .  Prom  these  equations  and  the 

substitutions  p  =  R  sin0  +0  ,  z  =  R  cos0  +0  ,  we  have 

2  2 

that ,  on  £  , 

X  =  a2  [c  +2R2c  -3R2c  sin20+***]  ,  (13) 

1  2  2 

i \>  =  a3[cf  +  R2(c’  -  4c'  R2)sin20  +  5c'  R4  sin4  0+  •  •  •  ]  . 

12  3  3 

(14) 

Comparison  of  (13)  with  (8)  yields 

c  =  ( 1/6 ) y R-2 [ 36-2 ( 1+k) ]  , 

i 

c  =  (l/6)uR--(l+k)  ,  (15) 

2 

c  =  c  =  •  •  •  =  0  . 

3  4 

Similarly,  from  (14)  and  (10)  we  obtain 

c’  =  0  , 
i 

c '  =  yVR-4 [6+(4-3k)/5]  > 

(16) 

c'  =  R"6yV(l+3k)/5  , 

3 

c '  =c!  =*..=0. 

4  5 
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By  choosing 


3  =  2(l+k)/3 
6  =  -  (4-3k)/5 


(17) 


we  make  both  c 
obtain  from  (7) 

i 


and  c’  vanish.  Correspondingly,  we 
2 

(li) , 


-h  =  1  +  u(2-k)/3  , 

-  |iR-2  [10+(3k-lt)V2)/5  . 


(18) 


To  summarize:  The  parametric  equations  of  the 
shell  in  V-  in  terms  of  intrinsic  coordinates  (0,<j),t) 

are 

p  =  R  sin0{l  +  \  a2  [1  +  y(2-k)/3  +  k( 2-y ) cos 2 0 ] }  , 

(19) 

z  =  R  cos0{l  +  ^  a 2  [k( 2-y ) cos 2 0  -  y(l+k)/3])  ,  (20) 

$  =  (}>  +  afit 

=  (J)  +  ayR"1  (1-a2  [2+(3k-4)V2/5]  )t  ,  (21) 

T  =  At 

=  V{1  -  h  a2  [1+y (2-k)/3] )t  .  (22) 

The  interior  metric  is  given  by  (1A),  with 

A(p,z)  =  (1/6)  a2y  R-2  ( 1+k )  (  2z  2  -p  2  )  ,  (23) 


'  • 
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\p(p  ,z)  =(l/5)a3yVR"3(l+3k)p2(p2-4z2)  .  (24) 


2.  Exterior  Curvature  3-tensor  K+,  of  Z  With 

- =i -  ab  —  - 

Respect  to  V+  . 

The  coordinates  in  the  exterior  manifold  V+ 

Oi 

are  x  =  (r,0,6,t)  ;  the  intrinsic  coordinates  of  Z  are 
=  (0  j(j)jt)  .  The  equation  of  Z  is  given  by  (V.19): 

r  =  r(0)  a  R( l+a2kV2 cos 2 0  )  .  (25) 


The  exterior  metric  is  given  by  (V.17)  or,  dis¬ 
carding  terms  of  0  , 

( ds  2 ) +  =  ga0dxadx3 

=  ( r  2  +  a2  cos  2  0  )  ( r  2  -2mr+a2  )  ~ 1  dr  2  +  (rz+az cos2 0 )d0z  + 

+  (r  2  +  a2  +  2ma2  sin2  0/r )  sin2  0d(J)  2  + 

+  [4  ma  sin2  0  ( 1-a2  cos 2  0/r 2 ) /r  ]  dcj)dt  - 

-  (l-2m/r  +  2ma2 cos 2 0/r 3 ) dt z  „  (26) 

The  unit  outward  normal  n  to  Z  is  deter- 

a 


mined  by 


n  =  b9  [r-r  (  0  )  ]  /'9xa 
a 

=  b ( 1 , 2a 2RkV2 sin0  cos0,  0,  0)  ,  (27) 


where  b 


is  a  normalizing  factor. 


' 
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1  =  n  na  =  g1]-(n  )2  +  g22(n  )2 
a  l  2 


y 


1.6.3 


b  =  (g11) 


1  1  \  — 2 


=  V  1  [1+^a2 ( 1-ky ) cos 2 0 ]  . 


Prom  the  defining  equation  (V.26), 


K 


ab 


-  na[82xa/Ha3?b  +  r^Y(3x3/9^a)OxY/95b) 


and  (25),  (26),  we  find 


K 


+ 


2  2 


n  r '  '  ( 0 )  -  nd[r  (r’)2+2r  r'+T  ] 
1  1 1  ,a  1 2  }a  2  2 }a 


n  r 
1 


;  1 


n1  [3g(r  ? )  2g  +  r  1  g  -kg 


11,1 


112 


n2[-^(r’)2g  +r’g  +  k  g  ]  . 

112  221  222 

5  5  5 


Now 


n  =  V-1  [1  +  k  a2 ( 1-ky ) cos  2 0 ]  , 
n1  =  V[1  -  k  a2 ( 1-ky ) cos 2 0 ]  , 

n2  =  2a2kV  sin0  cos0/R  , 

r’  =  -  2a2kV2R  sin0  cos0  , 

( r  '  )  2  =  0 


(28) 

]  (29) 

a=l ,  2 
]  - 

2  1 


"•  ;  . .  •  .  ■  i : 

CS)  ((d36\)ti’H*3«\R*«)t  ;.I  +  d3Ss36\tx*8!  i  -  -  *r 

..  S  I 
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ri;  =  -  2a2kV2R(cos2 0-sin20 )  , 


g  =  -  2a2  sin0  cos0/V2  on  E  , 
112 


g  =  2R(l+a2kV2  cos20)  on  E  , 

2  2  1 

3 


and 


g  =  -  2a2  sin0  cos0  on  E  . 

2  2  2 

3 


Hence 


K  +  =  VR  {1  -  3g  a2  [4k+ ( l+ky-10k) cos  2 0 ]  }  .  (30A) 

2  2 

Similarly , 


K 2  3  ^  n  (9x  /9e ^ (gab , 3+ga3,b  sb3,a^ 


a  3b=l 3  2 


=  0  , 


(30B) 


K  +  =  0  , 

2  4 


(  3  o  C  ) 


K  +  =  k  n!g  +  %  n2g 

33  331  332 

3  3 


=  VR  sin20{l  -  \  a2  [  ( l-6k+ky ) cos 2 0+y  sin20]} 


(  3  0  D ) 


K  +  =  -  ^  ay  V  sin20{l  -  \  a  2 ( 7  +  12k-5ky ) cos  2 0 }  , 

3  4 

(30E) 

K  +  =  -  h  (yV/R)  {1  -  3*5  a2(7  +  ^k-5ky)cos20}  .  (30F) 

4  4 


} 
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The  trace  Is 

K+  E  (K2  )+  +  (K3  +  K1* ) + 

2  3  4 

=  (V/R) {1  -  \  a2  [3 ( l-2k-ky ) cos2 0  +  4k]}  + 

+  (  2RV)  ~ 1  {1+V2  -  ^  a2  [2+ ( 2-y  )  cos  2  6  + 

+  k( -4+2y+3y 2 ) cos2 0 ] }  .  (31) 

3-  Extrinsic  Curvature  3-tensor  K~  of  Z  With  Respect 

- e -  ab  —  - - - 

to  V-  . 

The  coordinates  for  the  interior  region  V~  are 

OL 

x  =  (p,z,$,T).  Representing  the  shell  Z  by  an  equa¬ 
tion  of  the  form  p  =  p(z),  we  may  use 

na  =  (z  ,T)  (32) 

as  intrinsic  coordinates.  We  first  calculate  K  ,  K  .  , 
etc.,  and  then  make  an  intrinsic-coordinate  transformation 
(z,§,T)  ->  (0,cf),t)  to  obtain  KqQ  ,  etc. 

The  equation  of  Z  ,  p  =  p(z)  ,  is  derived  from 
(19),  (20)  by  eliminating  0  .  The  result  is 

p2  =  R2  +  a2R2  [l  +  y(2-k)/3]  - 

-  [  1  +  a2 ( l+y+yk-2k ) ] z  2  .  (33) 

The  unit  outward  normal  n  to  Z  is  given  by 

0 L 


r 

■  . 


ct 


•’  —  '  o  ;  3‘j. 

' 
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=  b3  [p-p (z ) ] /3xa 

=  b(l,-p'(z),0,0)  ,  (34) 

and 

1  =  n  na  =  b2gx 1  [l+( p ' ) 2 ]  (35) 

vX 

where  use  has  been  made  of  the  form  of  the  interior 
metric  (1A),  and  g11  is  to  be  evaluated  on  £  . 

From  the  defining  equation  for  [V. (26)] 

and  the  expressions  for  the  4-dimensional  Christoffel 
symbols  calculated  from  (1A),  we  obtain 

K~z  =  -  p  "  {g1  1  [l+(p  '  )2  ]  >_ls  +  (g1  1  U  +  (p’ )2]  F(P'AZ-Ap)  , 

K-,  =  K__  =  0  , 

z$  zT  s 

K”  =  tl+(p'  )2]  }“^(l-2A-pA  +pp’A  }  , 

<3?  11  p  A 

K~T  =  h  (gx x [i+(p' )2] }~^(^p-p’^z)  , 

KTT  =  (p’VAp)  > 

where  g1 1  ,  X,  Ap  =  3A/3p  ,  etc.,  are  to  be  evaluated  on 
E  .  With  the  help  of  (23),  (24),  (33),  and 

ii  =  (g  )-i  =  i  +  2X  , 

i  i 

it  is  a  straightforward  matter  to  calculate  K  etc., 

Zl  Z I  y 

as  functions  of  z  ,  or  of  0  through  (20).  In  this  way. 


we  find 
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K  =  (R  sinz@)  1  (1  +  3g  a2  [  ( l+4y+2k+2ky )  sin2  0  - 

Zj  Zj 


-  6k  -  y(5-4k)/3]  }  ,  (36) 

K“$  =  R  sin2 0 (1  +  h  a2  [l+y(5+2k)/3  - 

-  ( l+4y-6k+6ky ) cos  2 0 ] }  ,  (37) 

K“t  =  2a3yV(l+3k) ( 1/5-cos  2 0 ) sin2 0  ,  (38) 

K“t  =  a2R"1y(l+k)(sin20-2/3)  .  (39) 


We  now  make  the  intrinsic-coordinate  transforma¬ 
tion  na  =  (Zj$,T)  ->  E,a  =  (0,<|>, t),  defined  by  (20,  (21) 
and  (22),  and  use  the  rule  of  tensor  transformation 

Kab  =  K',b,(3na'/35a)(3rib'/3<f ) 
which  now  reads : 


K-0  =  (3z/30)2K"z  , 

(40) 

o 

II 

-P 

1  CD 
« 

II 

-e- 

1  CD 

w 

(41) 

K.  .  =  K  .  , 
cj xp  $$ 

(42) 

K,  ,  =  a^K . ,  +  AK,m  , 

cj) t  $$  $T  5 

(43) 

K"  =  a2ft2K~  +  A2K~m  . 
tt  $$  TT 

(44) 

Since 


y 


(3z/30) 


-  R  sin0{l  +  %  a2  [ 3W( 2-y ) cos 2 0-y ( 1+k) /3] } 
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n  =  uR“2{l-a2 [2+(3k-4)V2/5] }  , 


A  =  V  +  0  , 

2 

eqs . ( 40 ) - ( 44 )  yield,  when  (  36)  —  (  39)  are  substituted  into 
them,  the  following  non-vanishing  components  of  K&b  , 
referred  to  the  intrinsic  coordinates  (0,4>,t): 

K  =  R{1  +  %  a2  [  ( l+4y-10k+8ky ) sin2 0  +  6k  - 
2  2 

-  y  ( 7+l6k)/3] )  ,  (45A) 

K~  =  R  sin2  0 { 1  +  \  a2 [l+y( 5  +  2k)/3  - 

3  3 

-  ( l+4y-6k+6ky ) cos 2 0 ] }  ,  (45B) 


K  =  ay  sin20{l  -  \  a2 [ 3( l-2k) /5+ ( lly+8ky ) /15  + 

3  4 

+  ( 5+6k-6ky )cos20 ] }  ,  (450) 


K  =  a2yR_1 [ ( 1+k-ky ) sin2 0-2( 1-y ) ( l+k)/31  •  (45D) 

4  4 

Using  (6)  for  the  contravariant  components  of 
the  intrinsic  metric  gal3  (=  3ga^)  ,  we  obtain  the  trace 


K  =  (Kp 

d 

=  |  {1  -  h  a2 [-y(l+4k)/3  +  2k  + 

xl 


( 46 ) 


+  2( l+p-2k+ky ) cos2 6  ]  }  . 


cl 

4 .  Shell  Energy-Stress  Tensor  S  ,  Energy  Density  a 


u 

and  Angular  Velocity  co  :  Correct  to 

o  . 

From  Chap.  V,  we  have 

i 

8ttS.  .  =  g.  .  y  -  Y  -  •  j  y  •  -  =  .  - 

ij  &ij  '  ij  5  ij  ij 

k7  .  , 

ij 

(47) 

S .  .  u^  =  -  o  u .  , 

(48) 

with 

u2  =  0  , 

(49) 

u3  =  CO  u4  . 

(50) 

Eqs.  (47),  (48)  imply 

Yjj  =  (y  +  87ra)ui  , 

*\ 

CD 

•H 

Y  co  +  y  =  (Y  +  8Tra)(g  co+g  ) 

33  34  3334 

3 

(51) 

y  co  +  y  =  (Y+87Ta)(g  co+g  ) 

34  44  3444 

$ 

(52) 

so  that,  upon  dividing. 

0  =  co2  [g  y  “g  Y  ]  +  oo  [g  y 

3433  3334  4433 

g  Y  ] 

3  3  4  4 

+ 

+  [g  Y  -  g  Y  ]  a 

44  34  34  44 


which  is  equivalent  to 


oj  =  Y3/(<jOY£t+Y‘i-Yd)  • 

4  3  4  3 


(53) 
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All  the  above  equations  are  exact.  In  the  lowest 

approximation,  where  0  -  terms  are  neglected,  the  intrin- 

2 

sic  metric  of  Z  ,  referred  to  coordinates  0,<j),t  is  given 

a  + 

by  (V.29),  and  y  ,  y  ^  >  calculated  from  (neglecting 

0  -  terms)  of  Sects.  2  and  3,  are  very  simple.  Thus, 

2 


Y2  =  Y3  =  ( 1/R) ( V-l)  , 

2  3 

y3  =  -  h  (ay/R2V) [V(2+V)+y] 

4 

y4  =  3  ay  sin20/2V  , 

3 

y4  =  y/2RV  , 

4 


(other  components  zero)  and 


(54) 


y  =  ( 1/2RV) ( 1+3V2-4V)  .  (55) 


From  (47)  we  now  get  the  following  non-vanishing  components 
of  the  shell  energy-stress  tensor  (to  0  )  : 


8tt  S2  =  8 it  S3 

2  3 

=  ( 1/2RV ) ( 1-V ) 2  ,  (56A) 

8tt  S3  =  ^  (ay/R2V)  [V(2+V)+y]  ,  (56B) 

4 

8 7T  S4  =  -  3  ay  sin20/2V  ,  (56C) 

3 

8tt  S4  =  -  2  ( 1-V )  /R  .  (56D) 

4 


From  (53)  and  (54),  we  obtain  for  the  angular 
velocity  of  the  shell  (to  0  )  , 


- 


' 
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go  =  -  ay ( 1+2V) /R( 1+2V-3V2 )  .  (57) 

The  proper  energy  density  of  the  shell  is  most 
easily  calculated  from  (52),  with  the  result  (to  0  ) 

i 


a  =  (1-V)/4ttR  . 


(58) 


5*  Shell  Energy  Density  a  and  Angular  Velocity  to  , 
Correct  to  0 


Prom  (48)  and  (47),  one  gets 
8ttc 

au 

a  b 


o  o  a  b 
ott  S  ,  u  u 
ab 


Y  -  Yabu  u 


,  ir- \  /  Tr+  a  b ,  ,  T,—  a  b 

(R  -K  )  (Kabu  u  )  +  Kabu  u 


(59) 


a 

u  u 
a 


We  have  from  (49),  (50),  and  the  condition 

1  , 


-rr  a  b 

K  ,u  u 
ab 


4  N  2 


(u*) 


(w2K  +2ooK  +K  ) 

(u4)2  , 

(60) 

33  34  44 

-  (g  +2mg  +0JZg 

r1  . 

(61) 

4  4 


3  4  3  3 


Now  we  note  that,  to  evaluate  (60)  to  0  ,  one 

3 

needs  w  only  to  0  ,  which  is  already  given  in  (57). 

i 

With  the  help  of  (5)  for  the  intrinsic  metric  of  Z  refer¬ 
red  to  6,(J),t,  and  of  the  results  of  Sects.  2,3  for 
+ 


K 


ab 


we  find 


' 


- 
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(u4)z  =  V“z (l+az [1+y  sinz0-yk  cos 2 0  +  9y 2 sin2 0  V2/D2]}  , 

(62) 

032K+  +  2oa  K+  +  K+ 

33  34  44 

=  3§  (  yV/’R)  {-1  +  a2  [2y  sin2  0  (  2+9  V2 /D+9  V  VD2  )  + 

+  h.  (7+4k-5ky)cos20] }  ,  (63) 

oj2K“  +  2 id  K"  +  K~ 

33  34  44 

=  a2uV2R-1 {(l+9yV2D-2)sin26  +  [k  sin2 0-2 ( X+k) /3 ] }  , 

(64) 


where 


D  e  1  +  2V  -  3V2  .  (65) 

Using  these  and  (60)a  we  obtain 

K^buaub  =  -  (y/2VR){l  -  3^  a 2 ( 5+4k-3ky ) cos  2 0  - 

-  2a2 ( 1+5V+6V2+6V3 )/D( 1+3V) }  ,  (66) 

K“buaub  =  a2yR_1 (9yV2sin20/D2+(l+k) (1-3  cos20)/3> 

(6?) 

Substituting  (66),  (67),  (31),  (46)  into  (59),  one  then 
finds  for  the  energy  density  of  the  shell  (to  0  ) 


' 
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87ra  =  2(1-V)/R  -  \  a2VR  1  [  (  3y-4  +  6ky  +  8k)  cos  2  0  + 

+  9y2(l+2V-3V2)-1  sin2 0  -  (2+3y+4k)]  - 

-  a ZR~ 1  [ ( 2+3y-4k+3ky ) cos  2 0  +  2k  -  u(2+5k)/3]. 

:  .  (68) 

We  now  use  (51)  to  calculate  the  shell  angular 
velocity  03  to  0  ,  writing 

3 

oj  =  ao3  +03  (69A) 

1  3 

where  from  (57) 

03  =  -  y ( 1+2V ) R" 1  ( 1+2V-3V2  ) “ 1  ,  (69B) 

i 

and  03  is  a  higher  order  term  to  be  determined.  Eqs . 

3 

(51)  and  (69A)  yield 

03  [y  -(87Ta+y)g  ]  =  (87ra+y)(ao3  g  +  g  )  - 

3  33  33  1  3  3  34 

-  ( ao3  Y  +  y  )  ,  (70) 

1  3  3  3  4 

in  which  g  ^  are  given  by  (5),  o  by  (68),  and  Yab  and 
Y  can  be  obtained  from  the  results  of  Sects.  2,  3  for  the 
extrinsic  curvatures.  After  some  calculations,  one 
arrives  at 

tf 

2 

2R ( 1+2V-3V2  )  03 

3 

=  a3y( 1+2V-3V2 ) +  3V2 ( 3+ky ) cos  2 0  + 


,  co  + 


-  i. 
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+  4V[ V5+6y/5-3kV2/5+2V2(l+3k)cos20] ) 

+  3a3yV2 {2+6yV2-l8y2V2 (1+2V-3V2 )_1  + 

+  ( 2-lly-12k+l4ky+6y 2-3ky 2 ) cos  2 0  + 

+  l8y  2 V2 ( 1+2V-3V2 ) ~ 1  cos  2  0  + 

+  V[2+2y(7+4k)/3-2(l+6y-6k+8ky)cos20]  } 
Eqs  .  (69A,B,C)  determine  a)  to  0 

3 


(69C) 
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